TOP QUARK PHYSICS 



Conveners: M. Beneke, I. Efthymiopoulos, M.L. Mangano, J. Womersley 

Contributing authors: A. Ahmadov, G. Azuelos, U. Baur, A. Belyaev, E.L. Berger, W. Bernreuther, 
E.E. Boos, M. Bosnian, A. Brandenburg, R. Brock, M. Buice, N. Cartiglia, F. Cerutti, A. Cheplakov, 
L. Chikovani, M. Cobal-Grassmann, G. Corcella, F. del Aguila, T. Djobava, J. Dodd, V. Drollinger, 
A. Dubak, S. Frixione, D. Froidevaux, B. Gonzalez Pineiro, Y.P. Gouz, D. Green, P. Grenier, 
S. Heinemeyer, W. Hollik, V. Ilyin, C. Kao, A. Kharchilava, R. Kinnunen, V.V. Kukhtin, S. Kunori, 
L. La Rotonda, A. Lagatta, M. Lefebvre, K. Maeshima, G. Mahlon, S. Mc Grath, G. Medin, 
R. Mehdiyev, B. Mele, Z. Metreveli, D. O'Neil, L.H. Orr, D. Pallin, S. Parke, J. Parsons, D. Popovic, 
L. Reina, E. Richter-Was, T.G. Rizzo, D. Salihagic, M. Sapinski, M.H. Seymour, V. Simak, L. Simic, 
G. Skoro, S.R. Slabospitsky, J. Smolik, L. Sonnenschein, T. Stelzer, N. Stepanov, Z. Sullivan, T. Tait, 
I. Vichou, R. Vidal, D. Wackeroth, G. Weiglein, S. Willenbrock, W. Wu 



1. INTRODUCTION 

The top quark, when it was finally discovered at Fermilab in 1995 [[TJ, ^, ||], completed the three- 
generation structure of the Standard Model (SM) and opened up the new field of top quark physics. 
Viewed as just another SM quark, the top quark appears to be a rather uninteresting species. Produced 
predominantly, in hadron-hadron collisions, through strong interactions, it decays rapidly without form- 
ing hadrons, and almost exclusively through the single mode t— >Wb. The relevant CKM coupling V t b is 
already determined by the (three-generation) unitarity of the CKM matrix. Rare decays and CP violation 
are unmeasurably small in the SM. 
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Q_i' Yet the top quark is distinguished by its large mass, about 35 times larger than the mass of the 

next heavy quark, and intriguingly close to the scale of electroweak (EW) symmetry breaking. This 
unique property raises a number of interesting questions. Is the top quark mass generated by the Higgs 
mechanism as the SM predicts and is its mass related to the top-Higgs- Yukawa coupling? Or does it play 
an even more fundamental role in the EW symmetry breaking mechanism? If there are new particles 
lighter than the top quark, does the top quark decay into them? Could non-SM physics first manifest 
itself in non-standard couplings of the top quark which show up as anomalies in top quark production 
and decays? Top quark physics tries to answer these questions. 

Several properties of the top quark have already been examined at the Tevatron. These include 
studies of the kinematical properties of top production [Q], the measurements of the top mass [J5[ g], of the 
top production cross-section [80, the reconstruction of ti pairs in the fully hadronic final states H [l0[ ], 



the study of r decays of the top quark [11], the reconstruction of hadronic decays of the W boson from 
top decays [12], the search for flavour changing neutral current decays [O], the measurement of the 
W helicity in top decays [[14]], and bounds on ti spin correlations JTB]]. Most of these measurements are 
limited by the small sample of top quarks collected at the Tevatron up to now. The LHC is, in comparison, 
atop factory, producing about 8 million ti pairs per experiment per year at low luminosity (10fb -1 /year), 
and another few million (anti-)tops in EW single (anti-)top quark production. We therefore expect that 
top quark properties can be examined with significant precision at the LHC. Entirely new measurements 
can be contemplated on the basis of the large available statistics. 

In this chapter we summarize the top physics potential of the LHC experiments. An important 
aspect of this chapter is to document SM model properties of the top quark against which anomalous 
behaviour has to be compared. In each section (with the exception of the one devoted to anomalous 
couplings) we begin by summarizing SM expectations and review the current theoretical status on a 
particular topic. This is followed by a detailed description of experimental analysis strategies in the 

*To appear in the Report of the "1999 CERN Workshop on SM physics (and more) at the LHC". 



context of the ATLAS and CMS experiments. Particular emphasis is given to new simulations carried 
out in the course of this workshop. In detail, the outline of this chapter is as follows: 

In Section ^] we summarize SM precision calculations of the top quark mass relations and of the 
total top quark width. We then recall the importance of the top quark mass in EW precision measure- 
ments. We discuss, in particular, the role of EW precision measurements under the assumption that a SM 
Higgs boson has been discovered. 

Section |] deals with the tt production process: expectations for and measurements of the total 
cross section, the transverse momentum and tt invariant mass disuibution are discussed. A separate 
subsection is devoted to EW radiative collections to tt production, and to radiative corrections in the 
Minimal Supersymmetric SM (MSSM). 

The prospects for an accurate top quark mass measurement are detailed in Section |] Next to 
"standard" measurements in the lepton+jets and di-lepton channels, two mass measurements are dis- 
cussed that make use of the large number of top quarks available at the LHC: the selection of top quarks 
with large transverse momentum in the lepton+jets channel and the measurement of iJ/ip correlations 
in t—>£J/ifiX decays. This decay mode appears to be particularly promising and the systematic uncer- 
tainties are analyzed in considerable detail. 

Single top quark production through EW interactions provides the only known way to directly 
measure the CKM matrix element V t b at hadron colliders. It also probes the nature of the top quark 
charged current. In Section |5] the SM expectations for the three basic single top production mechanisms 
and their detection are documented, including the possibility to measure the high degree of polarisation 
in the SM. 

The issue of top quark spin is pursued in Section |6] Here we summarize expectations on spin 
correlations in ti production, the construction of observables sensitive to such correlations and the results 
of a simulation study of di-lepton angular correlations sensitive to spin correlations. Possible non-SM 
CP violating couplings of the top quark can be revealed through anomalous spin-momentum correlations 
and are also discussed here. 

As mentioned above, the search for anomalous (i.e. non-SM) interactions is one of the main moti- 
vations for top quark physics. In Section 7. the sensitivity of the LHC experiments to the following cou- 
plings is investigated: gti couplings and anomalous Wtb couplings in top production, flavour-changing 
neutral currents (FCNCs) in top production and decay. 

Section |T] is devoted to rare top decays. The SM expectations for radiative top decays and FCNC 
decays are documented. Decay rates large enough to be of interest require physics beyond the SM. 
The two Higgs Doublet Models, the MSSM and generic anomalous couplings are considered explicitly 
followed by ATLAS and CMS studies on the expected sensitivity in particular decay channels. 

Finally, the measurement of the top quark Yukawa coupling in tiH production is considered (Sec- 
tion [)]). The SM cross sections are tabulated in the various production channels at the LHC. For the case 
of a low mass Higgs boson, the results of a realistic study using a simulation of the ATLAS detector are 
discussed. 

The following topics are collected in the appendices: 6-quark tagging and the calibration of the 
jet energy scale in top events; the direct measurement of the top quark spin (as opposed to that of a top 
squark) and and of top quark electric charge; the total cross section for production of a fourth genera- 
tion heavy quark; a compendium of Monte Carlo event generators available for top production and its 
backgrounds. 

The internal ATLAS and CMS notes quoted in the bibliography can be obtained from the collab- 



orations' web pages [[If], |17j. Updated versions of this document, as well as a list of addenda and errata, 



will be available on the web page of the LHC Workshop top working group [18]. 



2 



2. TOP QUARK PROPERTIES AND ELECTROWEAK PRECISION MEASUREMENTS^] 

The top quark is, according to the Standard Model (SM), a spin- 1/2 and charge-2/3 fermion, transform- 
ing as a colour triplet under the group SU(3) of the strong interactions and as the weak-isospin partner 
of the bottom quark. None of these quantum numbers has been directly measured so far, although a 
large amount of indirect evidence supports these assignments. The analysis of EW observables in Z° 



decays [19] requires the existence of a T3 = 1/2, charge-2/3 fermion, with a mass in the range of 
170 GeV, consistent with the direct Tevatron measurements. The measurement of the total cross section 
at the Tevatron, and its comparison with the theoretical estimates, are consistent with the production of 
a spin- 1/2 and colour-triplet particle. The LHC should provide a direct measurement of the top quantum 
numbers. We present the results of some studies in this direction in Appendix ||. 

2.1 Top quark mass and width 

In addition to its quantum numbers, the two most fundamental properties of the top quark are its mass m f 
and width T t , defined through the position of the single particle pole m\ = 774 — iFt /2 in the perturbative 
top quark propagator. In the SM m t is related to the top Yukawa coupling: 



lit 



( f i)=2 3 / 4 G 1 J 2 m t (l + 5 t (ri), (1) 



where 5t(n) accounts for radiative corrections. Besides the top quark pole mass, the top quark MS mass 
rrit(fJi) is often used. The definition of m"t(/i) including EW corrections is subtle (see the discussion in 
[|20|]). As usually done in the literature, we define the MS mass by including only pure QCD corrections: 

mt(fi) = m t (1 + (5q C d(/u)) _1 • (2) 

The conversion factor <5qcd(/^) is very well known [jH]]. Defining m t = m t {mt) and a s = ot?f (yntj/n, 
we have 

5qcd (mt) = ~ a s + 8.2366 a 2 s + 73.638 a* + .. . 

= (4.63 + 0.99 + 0.31 + 0.1lljj;£)% = (6.05±g;^)%. (3) 

This assumes five massless flavours besides the top quark and we use a s = 0.03475 which corresponds 
to a^ s (mz) = 0.119 and m t = 165 GeV. The error estimate translates into an absolute uncertainty 
of ±180 MeV in m t — m t and uses an estimate of the four-loop contribution. Note that the difference 
between the two mass definitions, mt — m t , is about 10 GeV. This means that any observable that is 
supposed to measure a top quark mass with an accuracy of 1-2 GeV and which is known only at leading 
order (LO) must come with an explanation for why higher order corrections are small when the observ- 
able is expressed in terms of that top quark mass definition that it is supposed to determine accurately. 
We will return to this point in Section |] 

The on-shell decay width Tt is less well known, but the theoretical accuracy (< 1%) is more than 
sufficient compared to the accuracy of foreseeable measurements. The decay through t—>bW is by far 
dominant and we restrict the discussion to this decay mode. It is useful to quantify the decay width in 
units of the lowest order decay width with Mw and m^ set to zero and \ V t b\ set to 1: 



_ G F ml 

Incorporating Mw the leading order result reads 



1.76 GeV. (4) 



r LO (t^bW)/\V tb \ 2 = T 1 - 3-f + 2^ = 0.885 T = 1.56GeV. (5) 

\ mt / 



'Section coordinators: M. Beneke, G. Weiglein. 
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Table 1: Corrections to the top quark width To (Mw = 0, lowest order) in units of To. The best estimate of T(t^bW) /\ Vtb\ 



is obtained by adding all corrections together. Parameters: a s = 0.03475, My 



Mw 7^ correction at lowest order, see 
a s correction, Mw = 
a s correction, Mw ^ correction 



at correction, M w = [22, E3fl 



a 2 s correction, Mw / correction |p3l] 



EW correction [ 24 ] 



1.4 GeV and m t = 175 GeV. 



-11.5% 
-9.5% 
+1.8% 
-2.0% 
+0.1% 
+1.7% 



The correction for non-vanishing bottom quark mass is about —0.2% in units of Tq. Likewise corrections 
to treating the W boson as a stable particle are negligible. Radiative corrections are known to second 
order in QCD and to first order in the EW theory. Table [j] summarises the known corrections to the 
limiting case @. Putting all effects together we obtain: 



T(t^bW)/\V tb \ 2 » 0.807r = 1.42 GeV. 



(6) 



The top quark lifetime is small compared to the time scale for hadronisation [|25|]. For this reason, top- 
hadron spectroscopy is not expected to be the subject of LHC measurements. 



2.2 Role of m t in EW precision physics 

The EW precision observables serve as an important tool for testing the theory, as they provide an impor- 
tant consistency test for every model under consideration. By comparing the EW precision data with the 
predictions (incorporating quantum corrections) within the SM or its extensions, most notably the mini- 
mal supersymmetric extension of the Standard Model (MSSM) [Eq], it is in principle possible to derive 
indirect constraints on all parameters of the model. The information obtained in this way, for instance, on 
the mass of the Higgs boson in the SM or on the masses of supersymmetric particles is complementary 
to the information gained from the direct production of these particles. 

In order to derive precise theoretical predictions, two kinds of theoretical uncertainties have to 
be kept under control: the uncertainties from unknown higher-order collections, as the predictions are 
derived only up to a finite order in perturbation theory, and the parametric uncertainties caused by the 
experimental errors of the input parameters. The top quark mass enters the EW precision observables as 
an input parameter via quantum effects, i.e. loop corrections. As a distinctive feature, the large numerical 
value of m t gives rise to sizable corrections that behave as powers of mj. This is in contrast to the 
corrections associated with all other particles of the SM. In particular, the dependence on the mass of the 
Higgs boson is only logarithmic in leading order and therefore much weaker than the dependence on mt. 
In the MSSM large corrections from SUSY particles are only possible for large splittings in the SUSY 
spectrum, while the SUSY particles in general decouple for large masses. 

The most important m t -dependent contribution to the EW precision observables in the SM and 



the MSSM enters via the universal parameter Ap which is proportional to mf J27[], 
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where the limit m;,— »0 has been taken, sw (c\y) is t ne sin (cos) of the weak mixing angle, and S z (0) 
and S (0) indicate the transverse parts of the gauge-boson self-energies at zero momentum transfer. 

The theoretical prediction for Mw is obtained from the relation between the vector-boson masses 
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and the Fermi constant, 



M » I 1 "if J = ^ (1+Ar) ' <8) 



where the quantity Ar [ 28 ] is derived from muon decay and contains the radiative corrections. At one- 

c 2 

loop order, Ar can be written as Ar = Aa — ^ Ap + (Ar) n i, where Aa contains the large logarithmic 

s w 

contributions from the light fermions, and the non-leading terms are collected in (Ar) n i. 

The leptonic effective weak mixing angle is determined from the effective couplings of the neu- 
tral current at the Z-boson resonance to charged leptons, jj^ c = (^/2GfM^ ^ [gvln — dAl^ls], 
according to 



efr 4 V Re (g A ) J 

In sin 2 9^ the leading m t -dependent contributions enter via <5sin 2 9^ = — (c w s w )/(c w — s w )Ap. 

The precision observables M\y and sin 2 JP* are currently known with experimental accuracies of 
0.05% and 0.07%, respectively [ 19]. The accuracy in Mw will be further improved at the LHC by about 
a factor of three (see the EW chapter of this Yellow Report). Besides the universal collection Ap, there is 
also a non-universal correction proportional to m 2 in the Zbb coupling, which however is less accurately 
measured experimentally compared to Mw and sin 2 O 1 ^. The strong dependence of the SM radiative 
corrections to the precision observables on the input value of nit made it possible to predict the value of 
nit from the precision measurements prior to its actual experimental discovery, and the predicted value 
turned out to be in remarkable agreement with the experimental result []|, g]. 

Within the MSSM, the mass of the lightest CP-even Higgs boson, m^, is a further observable 
whose theoretical prediction strongly depends on m t . While in the SM the Higgs-boson mass is a free 
parameter, nih is calculable from the other SUSY parameters in the MSSM and is bounded to be lighter 
than Mz at the tree level. The dominant one-loop corrections arise from the top and scalar-top sector via 
terms of the form Gprnf ln(mf mf 2 /m|) [p9Q. As a rule of thumb, a variation of m t by 1 GeV, keeping 
all other parameters fixed, roughly translates into a shift of the predicted value of by 1 GeV. If the 
lightest CP-even Higgs boson of the MSSM will be detected at the LHC, its mass will be measurable 
with an accuracy of about Am^ = 0.2 GeV [30]. 

Due to the sensitive dependence of the EW precision observables on the numerical value of nit, 
a high accuracy in the input value of m t is very important for stringent consistency tests of a model, 
for constraints on the model's parameters (e.g. the Higgs boson mass within the SM), and for a high 
sensitivity to possible effects of new physics. It should be noted that this calls not only for a high 
precision in the experimental measurement of the top quark mass, but also for a detailed investigation 
of how the quantity that is actually determined experimentally is related to the parameter nit used as 
input in higher-order calculations. While these quantities are the same in the simplest approximation, 
their relation is non-trivial in general due to higher-order contributions and hadronisation effects. A 
further discussion of this problem, which can be regarded as a systematic uncertainty in the experimental 
determination of m t , is given in Section [TJ 



2.3 Physics gain from improving Am t from Am t = 2 GeV to Am t = 1 GeV 

During this workshop the question was investigated of how much information one could gain from the 
EW precision observables by improving the experimental precision in m t from Am t = 2 GeV, reachable 
within the first year of LHC running (see Section fO|), to Anit = 1 GeV, possibly attainable on a longer 
time scale (see Section |4.6| ). 

In order to analyse this question quantitatively, we have considered the case of the SM and the 
MSSM and assumed that the Higgs boson has been found at the LHC. For the uncertainty in Aahad (the 
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Table 2: Comparison of the current theoretical uncertainty from unknown higher-order corrections (Atheo) in Mw and 
sin 2 O^g* with the parametric uncertainties from the error in Aai la d and mt. 





Atheo 


S{Aa bad ) = 0.00016 


Am t = 2 GeV 


Am t = 1 GeV 


AM w /MeV 


6 


3.0 


12 


6.1 


Asin 2 ^ x 10 5 


4 


5.6 


6.1 


3.1 



hadronic contribution to the electromagnetic coupling at the scale Mz) we have adopted J(Aahad) = 
0.00016, which corresponds to the "theory driven" analyses of [31]. 

Concerning the current theoretical prediction for Mw and sin 2 9 l ^E t in the SM, the theoretical un- 
certainty from unknown higher-order corrections has been estimated to be about AMw = Q MeV and 
A sin 2 ^gff 5 * = 4 x 10 5 |(32]. In Table || the theoretical uncertainties for Mw and sin 2 9^ t from un- 
known higher-order corrections are compared with the parametric uncertainty from the input parameters 
Aahad and mt for Am t = 2 GeV as well as Amt = 1 GeV. The parametric uncertainties from the other 
parameters, supposing that the SM Higgs boson has been found at the LHC in the currently preferred 
range, are negligible compared to the uncertainties from Aa n ad and mt. The resulting uncertainties in 
M\y and sin 2 9^ have been obtained using the parameterisation of the results for these quantities given 
in [33]. As can be seen in the table, for Am t = 2 GeV the parametric uncertainty in m t gives rise to the 
largest theoretical uncertainty in both precision observables. While for sin 2 9 l ^ t the uncertainty induced 
from the error in mt is comparable to the one from the error in Aa na d> f° r Mw the uncertainty from 
the error in m t is twice as big as the one from unknown higher-order corrections and four times as big 
as the one from the error in Acthad- A reduction of the error from Amt = 2 GeV to Amt = 1 GeV 
will thus mainly improve the precision in the prediction for Mw- The uncertainty induced in Mw by 
Amt = 1 GeV is about the same as the current uncertainty from unknown higher-order corrections. 
The latter uncertainty can of course be improved by going beyond the present level in the perturbative 
evaluation of Ar. 

In Fig. [I] the theoretical predictions for Mw and sin 2 B 1 ^ 1 (see [ 34 ] and references therein) are 
compared with the expected accuracies for these observables at LEP2/Tevatron and at the LHC (for the 
central values, the current experimental values are taken). The parametric uncertainties corresponding to 
5(Aahad) = 0.00016 and Am t = 2 GeV, Am t = 1 GeV are shown for two values of the Higgs boson 
mass, mn = 120 GeV and mjj = 200 GeV, and the present theoretical uncertainty is also indicated 
(here mjj is varied within 100 GeV< ma < 400 GeV and Am t = 5.1 GeV). The figure shows that, 
assuming that the Higgs boson will be discovered at the LHC, the improved accuracy in m t and My/ at 
the LHC will allow a stringent consistency test of the theory. A reduction of the experimental error in mt 
from Amt = 2 GeV to Am t = 1 GeV leads to a sizable improvement in the accuracy of the theoretical 
prediction. In view of the precision tests of the theory a further reduction of the experimental error in 
My/ and sin 2 9 1 ^ would clearly be very desirable. 

While within the MSSM the improved accuracy in m t and Mw at the LHC will have a similar 
impact on the analysis of the precision observables as in the SM, the detection of the mass of the lightest 
CP-even Higgs boson will provide a further stringent test of the model. The prediction for m^ within 
the MSSM is particularly sensitive to the parameters in the t-t sector, while in the region of large Ma, 
and large tan f3 (giving rise to Higgs masses beyond the reach of LEP2) the dependence on the latter two 
parameters is relatively mild. A precise measurement of m/j can thus be used to constrain the parameters 
in the t-t sector of the MSSM. 

In Fig. |H it is assumed that the mass of the lightest scalar top quark, , is known with high 
precision, while the mass of the heavier scalar top quark, mj 2 , and the mixing angle 9^ are treated as free 
parameters. The Higgs boson mass is assumed to be known with an experimental precision of ±0.5 GeV 
and the impact of Am t = 2 GeV and Am t = 1 GeV is shown (the theoretical uncertainty in the Higgs- 
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M w [GeV] 

Fig. 1: The SM prediction in the Mw-sin 2 9^g* plane is compared with the expected experimental accuracy at LEP2/Tevatron 
(AM W = 30 MeV, sin 2 6^ c ff pt = 1.7 x 10~ 4 ) and at the LHC (AM W = 15 MeV, sin 2 6^ c ff pt = 1.7 x 1(T 4 ). The theoretical 
uncertainties induced by <5(Aahad) = 0.00016 and Amt = 2 GeV (full line) as well as Amt = 1 GeV (dashed line) are 
shown for two values of the Higgs boson mass mu ■ 
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Fig. 2: Indirect constraints on the parameters of the scalar top sector of the MSSM from the measurement of mi at the LHC. 
The effect of the experimental error in m t is shown for Am t = 2 GeV and Amt = 1 GeV. 
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mass prediction from unknown higher-order contributions and the parametric uncertainties besides the 
ones induced by m 4 ~ 2 , 9^ and m t have been neglected here). The two bands represent the values of m^ 2 , % 
which are compatible with a Higgs-mass prediction of rrih = 120.5 ±0.5 GeV, where the two-loop result 



of J35| ] has been used (the bands corresponding to smaller and larger values of are related to smaller 
and larger values of the off-diagonal entry in the scalar top mixing matrix, respectively). Combining the 
constraints on the parameters in the scalar top sector obtained in this way with the results of the direct 
search for the scalar top quarks will allow a sensitive test of the MSSM. As can be seen in the figure, a 
reduction of Am t from Amt = 2 GeV to Am t = 1 GeV will lead to a considerable reduction of the 
allowed parameter space in the m^-O^ plane. 

3. tt PRODUCTION AT THE LHCf] 

The determination of the top production characteristics will be one of the first measurements to be carried 
out with the large statistics available at the LHC. The large top quark mass ensures that top production 
is a short-distance process, and that the perturbative expansion, given by a series in powers of the small 
parameter a s {m t ) ~ 0.1, converges rapidly. Because of the large statistics (of the order of 10 7 top quark 
pairs produced per year), the measurements and their interpretation will be dominated by experimen- 
tal and theoretical systematic errors. Statistical uncertainties will be below the percent level for most 
observables. It will therefore be a severe challenge to reduce experimental and theoretical systematic 
uncertainties to a comparable level. In addition to providing interesting tests of QCD, accurate studies of 
the top production and decay mechanisms will be the basis for the evaluation of the intrinsic properties 
of the top quark and of its EW interactions. An accurate determination of the production cross section, 
for example, provides an independent indirect determination of mt- Asymmetries in the rapidity distri- 



butions of top and antitop quarks [36] are sensitive to the light-quark partem distribution functions of the 
proton. Anomalies in the total ti rate would indicate the presence of non-QCD production channels, to be 
confirmed by precise studies of the top quark distributions (e.g. p T and ti invariant mass spectra). These 
would be distorted by the presence of anomalous couplings or s-channel resonances expected in several 
beyond-the-SM (BSM) scenarios. Parity-violating asymmetries (for example in the rapidity distributions 
of right and left handed top quarks) are sensitive to the top EW couplings, and can be affected by the 
presence of BSM processes, such as the exchange of supersymmetric particles. As already observed at 
the Tevatron ^J, the structure of the ti final state affects the direct determination of m t . Initial and 
final-state gluon radiation do in fact contribute to the amount of energy carried by the jets produced in 
the decay of top quarks, and therefore need to be taken into proper account when jets are combined to 
extract mt. The details of the structure of these jets (e.g. their fragmentation function and their shapes), 
will also influence the experimental determination of the jet energy scales (important for the extraction 
of mt), as well as the determination of the efficiency with which 6-jets will be tagged (important for the 
measurement of the production cross section). 

It is therefore clear that an accurate understanding of the QCD dynamics is required to make full 
use of the rich statistics of ti final states in the study of the SM properties of top quarks, as well as to 
explore the presence of possible deviations from the SM. In this section we review the current state of the 
art in predicting the production properties for top quark pairs (for a more detailed review of the theory of 
heavy quark production, see [^7|]). The study of single top production will be presented in Section |] 



3.1 Tools for QCD calculations 

Full next-to-leading-order (NLO, 0(a|)) calculations are available for the following quantities: 



1. Total cross sections [38] 



2. Single-inclusive px and y spectra [ f3"9| ] 

3. Double-differential spectra (m t i, azimuthal correlations A<I>, etc.) 



2 Section coordinators: M.L. Mangano, D. Wackeroth, M. Cobal (ATLAS), J. Parsons (ATLAS). 
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All of the above calculations are available in the form of Fortran programs [ |40| , |41| ], so that kinematical 
distributions can be evaluated at NLO [42] even in the presence of analysis cuts. 

Theoretical progress over the last few years has led to the resummation of Sudakov-type loga- 
rithms [ |43| ] which appear at all orders in the perturbative expansion for the total cross sections [|4|, fj5fl . 
More recently, the accuracy of these resummations has been extended to the next-to-leading logarithmic 
(NLL) level [46, 47]. For a review of the theoretical aspects of Sudakov resummation, see the QCD 
chapter of this report. As will be shown later, while the inclusion of these higher-order terms does not 
affect significantly the total production rate, it stabilises the theoretical predictions under changes in the 
renormalisation and factorisation scales, hence improving the predictive power. 

Unfortunately, the results of these resummed calculations are not available in a form suitable to 
implement selection cuts, as they only provide results for total cross-sections, fully integrated over all 
of phase space. The formalism has been generalised to the case of one-particle inclusive distributions 
in p8|], although no complete numerical analyses have been performed yet. 

The corrections of C(a|) to the full production and decay should include the effect of gluon 
radiation off the quarks produced in the top decay. Interference effects are expected to take place between 
soft gluons emitted before and after the decay, at least for gluon energies not much larger than the top 
decay width. While these correlations are not expected to affect the measurement of generic distributions, 
even small soft-gluon corrections can have an impact on the determination of the top mass. Matrix 
elements for hard-gluon emission in ti production and decay (j>p—>W + bW~bg, with t and i intermediate 
states) are implemented in a parton-level generator tf49|[. The one-gluon emission off the light quarks 
from the W decays was implemented, in the soft-gluon approximation, in the parton-level calculation 
of @]. 

The above results refer to the production of top quarks treated as free, stable partons. Parton- 
shower Monte Carlo programs are available (HERWIG [Q, PYTHIA [|52|], ISAJET []53|]) for a complete 
description of the final state, including the full development of the perturbative gluon shower from both 
initial and final states, the decay of the top quarks, and the hadronisation of the final-state partons. These 
will be reviewed in Appendix |d[ Recently, 0(a s ) matrix element corrections to the decay of the top 
quark (t— >Wbg) have been included in the HERWIG Monte Carlo [54]. The impact of these corrections 



will be reviewed in Sections 3.3 and 4.62 



3.2 Total tt production rates 

In this section we collect the current theoretical predictions for cross sections and distributions, pro- 
viding our best estimates of the systematic uncertainties. The theoretical uncertainties we shall consider 
include renormalisation (fi R ) and factorisation (/jl f ) scale variations, and the choice of parton distribution 
functions (PDF's); 

We shall explore the first two by varying the scales over the range /io/2 < /** < 2^0, where 
/i = [i R = fi F and 

• fio = m t for the total cross sections 



Mo = v m t + Pt f° r s i n gl e inclusive distributions 



Mo = J m-t + (Pt t + Pt f)/^ f° r double inclusive distributions 



In the case of PDF's, we shall consider the latest fits of the CTEQ [55] and of the MRST [56, 57] groups: 



• MRST (a s (M z ) = 0.1175, (k T ) = 0.4 GeV) (default) 

• MRST(<7 |) (a s = 0.1175, (k T ) = 0.64 GeV) 

• MRST( 5 T) («s = 0.1175, (for) = 0) 

• MRST(a s ||) (a s = 0.1125, (for) = 0.4 GeV) 
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: tt production rates. Left: scale dependence at fixed order (NLO, dashed lines in the lower inset), and at NLO+NLL 
lines). Right: PDF dependence. See the text for details. 



• MRST(a s tt) (a s = 0.1225, (k T ) = 0.4 GeV) 

• CTEQ5M (a s = 0.118) 

• CTEQ5HJ (a s = 0.118, enhanced weight for Tevatron high-£r jets) 

• CTEQ5HQ (a s = 0.118, using the ACOT heavy flavour scheme g].) 

All our numerical results relative to the MRST sets refer to the updated fits provided in []57|]. These give 
total rates which are on average 5% larger than the fits in [{Sfifl. The total it production cross section is 
given in Fig. ||, as a function of the top mass. As a reference set of parameters, we adopt /zq = m t and 
MRST. Full NLO+NLL corrections are included. The upper inset shows the dependence of the cross 
section on the top mass. A fit to the distribution shows that Aa/a ~ hAmt/mt. As a result, a 5% mea- 
surement of the total cross section is equivalent to a 1% determination of mt (approximately 2 GeV). 
As will be shown later on, 2 GeV is a rather safe estimate of the expected experimental accuracy in the 
determination of mt (1 GeV being the optimistic ultimate limit). It follows that 5% should be a minimal 
goal in the overall precision for the measurement of a(tt). The scale uncertainty of the theoretical pre- 
dictions is shown in the lower inset of Fig. ||. The dashed lines refer to the NLO scale dependence, which 
is of the order of ±12%. The dotted lines refer to the inclusion of the NLL corrections, according to the 
results of [47]. The solid lines include the resummation of NLL effects, but assume a different structure 
of yet higher order (NNLL) corrections, relative to those contained in the reference NLL results (this is 
indicated by the value of the A parameter equal to 2, see [47] for the details). The scale uncertainty, after 
inclusion of NLL corrections, is significantly reduced. In the most conservative case of A = 2, we have a 
±6% variation. A detailed breakdown of the NLO 0{a 2 s + a|) and higher-order C(a| 4 ) contributions, 
as a function of the scale and of the value of the parameter A, is given in Table |3[ A recent study [59] 
of resummation effects on the total cross section for photo- and hadro-production of quarkonium states 
indicates that allowing /j, R ^ \i F increases the scale dependence of the NLL resummed cross-sections 
to almost match the scale dependence of the NLO results [|60t]. Preliminary results of this study also 
suggest a similar increase of scale dependence in the case of it production, if \i R and \i F are varied inde- 
pendently. This dependence can however be reduced by replacing \i R with fj, F as the argument of a s in 
the sub-leading coefficients of the resummed exponent [|6l|]. 

The PDF dependence is shown on the right hand side of Fig. |3|, and given in detail for m t = 
175 GeV in Table (|. The current uncertainty is at the level of ±10%. Notice that the largest deviations 
from the default set occur for sets using different input values of a s (Mz)- The difference between the 
reference sets of the two groups (MRST and CTEQ5M) is at the level of 3%. It is interesting to explore 
potential correlations between the PDF dependence of top production, and the PDF dependence of other 
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Table 3: Resummation contributions to the total tt cross-sections (m t = 175 GeV) in pb. PDF set MRST. 







NLL resummed, A=2 


NLL resummed, A=0 


Mil = fJ-F 


NLO 
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NLO+NLL 




NLO+NLL 


mt/2 
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-12 


878 


m t 


796 


29 


825 


63 


859 


2m t 


705 


77 


782 


148 


853 



Table 4: Total tt cross-sections (m t = 175 GeV) in pb. NLO+NLL (A = 0). 



PDF 
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877 
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876 
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MRST a s t 
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942 


934 


CTEQ5M 


904 


886 


881 


CTEQ5HJ 


905 


886 


881 



processes induced by initial states with similar parton composition and range in x. One such example is 
given by inclusive jet production. Fig. |] shows the initial-state fraction of inclusive jet final states (with 
\rj\ < 2.5) as a function of the jet-E^ threshold. For values of E T ~ 200 GeV, 90% of the jets come 
from processes with at least one gluon in the initial state. This fraction is similar to that present in tt 
production, where 90% of the rate is due to gg collisions. On the right side of Fig. |] we show the double 
ratios: 



[a(ti)/a(jet,E T >E™*)] pDF 



(10) 



[a(tt)/a(jtt,E T >Ef™)] MmjT 

As the plot shows, there is a strong correlation between the PDF dependences of the two processes. The 
correlation is maximal for EJp m ~ 200 GeV, as expected, since for this value the flavour composition of 
the initial states and the range of partonic momentum fractions probed in the two production processes are 
similar. In the range 180 < Ef in < 260 GeV the PDF dependence of the ratio a{tt) /a (jet, E T > Ef in ) 
is reduced to a level of ±1%, even for those sets for which the absolute top cross-section varies by ±10%. 
The jet cross-sections were calculated [52] using a scale fi jet = E T = fi J et . If we vary the scales for tt 
and jet production in a correlated way (i.e. selecting ^ et / 1 [i 3 ^ 1 = pL tl / ' Hq), no significant scale depen- 
dence is observed. There is however no a-priori guarantee that the scales should be correlated. Unless 
this correlation can be proved to exist, use of the inclusive-jet cross section to normalise the tt cross sec- 
tion will therefore leave a residual systematic uncertainty which is no smaller than the scale dependence 
of the jet cross section. We do not expect this to become any smaller than the PDF dependence in the 
near future. 

Combining in quadrature the scale and PDF dependence of the total tt cross section, we are left 
with an overall 12% theoretical systematic uncertainty, corresponding to a 4 GeV uncertainty on the 
determination of the top mass from the total cross section. 



3.3 Kinematical properties of tt production 

We start from the most inclusive quantity, the top p T spectrum. The NLO predictions are shown in Fig. |5[ 
Here we also explore the dependence on scale variations and on the choice of PDF. The uncertainties are 
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Fig. 4: Left: initial state composition in inclusive jet events, as a function of the jet E T (\t]\ < 2.5). Right: PDF dependence of 
the top-to-jet cross-section ratio, as a function of the minimum jet E T . 




±15% and ±10%, respectively. The reconstruction of top quarks and their momenta, as well as the 
determination of the reconstruction efficiencies and of the possible biases induced by the experimental 
selection cuts, depend on the detailed structure of the final state. It is important to verify that inclusive 
distributions as predicted by the most accurate NLO calculations are faithfully reproduced by the shower 
Monte Carlo calculations, used for all experimental studies. This is done in Fig. ||, where the NLO 
calculation is compared to the result of the HERWIG Monte Carlo, after a proper rescaling by an overall 
constant K-factor. The bin-by-bin agreement between the two calculations is at the level of 10%, which 
should be adequate for a determination of acceptances and efficiencies at the percent level. 

Similar results are obtained for the invariant mass distribution of top quark pairs, shown in the 
plot of Fig. ||. The scale and PDF dependence of the NLO calculation are similar to those found for the 
inclusive px spectrum, and are not reported in the figure. 

Contrary to the case of inclusive p T and M tf spectra, other kinematical distributions show large 
differences when comparing NLO and Monte Carlo results [42]. This is the case of distributions which 
are trivial at LO, and which are sensitive to Sudakov-like effects, such as the azimuthal correlations or the 
spectrum of the ti pair transverse momentum p". These two distributions are shown in the two plots of 
Fig. 0. Notice that the scale uncertainty at NLO is larger for these distributions than for previous inclusive 
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Fig. 6: Comparison of NLO and (rescaled) HERWIG spectra. Left: inclusive top pr spectrum. Right: inclusive M t i spectrum. 
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Fig. 7: Left: azimuthal correlation between the t and t. Right: integrated transverse momentum spectrum of the top quark pair. 
Continuous lines correspond to the parton-level NLO calculation, for different scale choices; the plots correspond to the result 
of the HERWIG Monte Carlo. 



quantities. These kinematical quantities are in fact trivial at 0(ot 2 s ) (proportional to S-f unctions), and 
their evaluation at 0(a s ) is therefore not a true NLO prediction. The regions p"— >0 and Aifi—nr are 
sensitive to multiple soft-gluon emission, and the differences between the NLO calculation (which only 
accounts for the emission of one gluon) and the Monte Carlo prediction (which includes the multi-gluon 
emission) is large. The region p" 3> m t is vice-versa sensitive to the emission of individual hard gluons, 
a process which is more accurately accounted for by the full 0{o? s ) matrix elements included in the 
NLO calculation than by the Monte Carlo approach. Notice that the average value of p" is quite large, 
above 50 GeV. This is reasonable, as it is of the order of a s times the average value of the hardness of 
the process ((M ti ) ~ 540 GeV). It is found that this large transverse momentum is compensated by the 
emission of a jet recoiling against the top pair, with a smaller fraction of events where the comes 
from emission of hard gluons from the final state top quarks. The large-p" discrepancy observed in 
Fig. [7] should be eliminated once the matrix element corrections to top production will be incorporated 
in HERWIG, along the lines of the work done for Drell-Yan production in [63]. 

Emission of extra jets is also expected from the evolution of the decay products of the top quarks 
(6's, as well as the jets from the hadronic W decays). Gluon radiation off the decay products is included 
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AR n iets 
Fig. 8: Left: Distributions of the minimum invariant opening angle AR among the three hardest jets, with (HERWIG 6.1, 
solid line) and without (6.0, dotted) matrix element corrections in the t—>bW decay. Right: jet multiplicity distributions. 



in the shower Monte Carlo calculations. In the case of the latest version of HERWIG (v6. 1) [51], the 



emission of the hardest gluon from the b quarks is evaluated using the exact matrix elements [|54]]. This 
improvement, in addition to a few bug fixes, resolve the discrepancies uncovered in [|^] between an 
exact parton level calculation and previous versions of HERWIG The matrix-element corrections do not 
alter significantly most of the inclusive jet observables. As examples, we show in Fig. |8] the AR and the 
jet multiplicity distributions for events where both W's decay leptonically. More details can be found 

1. As can be seen, 



in [|64J. Jets are defined using the krp algorithm Q65Q, with radius parameter R 
the impact of the exact matrix element corrections is limited, mostly because the extra-jet emission is 
dominated by initial-state radiation. 

The impact on quantities which more directly affect the determination of the top mass remains 
to be fully evaluated. Given the large rate of high-E^ jet emissions, their proper description will be a 
fundamental ingredient in the accurate reconstruction of the top quarks from the final state jets, and in the 
determination of the top quark mass. A complete analysis will only be possible once the matrix element 
corrections to the tt production will be incorporated in the Monte Carlos. Work in this direction is in 
progress (G. Corcella and M.H. Seymour). 



3.4 Non-QCD radiative corrections to tt production 

The production and decay of top quarks at hadron colliders is a promising environment for the detection 
and study of loop induced SUSY effects: at the parton level there is a large center of mass energy s 
available and owing to its large mass, the top quark strongly couples to the (virtual) Higgs bosons, a 
coupling which is additionally enhanced in SUSY models. Moreover, it might turn out that SUSY loop 
effects in connection with top and Higgs boson interactions less rapidly decouple than the ones to gauge 
boson observables. 



To fully explore the potential of precision top physics at the LHC and at the Tevatron [66] to detect, 
discriminate and constrain new physics, the theoretical predictions for top quark observables need to be 
calculated beyond leading order (LO) in perturbation theory. Here we will concentrate on the effects of 
non-QCD radiative corrections to the production processes gg— >ti and qq—>ti, including supersymmetric 
corrections. When searching for quantum signatures of new physics also the SM loop effects have to be 
under control. The present SM prediction for tt observables includes the QCD corrections as discussed 
above and the EW one-loop contributions to the QCD tt production processes [67, 68, |69l]. The latter 



modify the gtt(qqg) vertex by the virtual presence of the EW gauge bosons and the SM Higgs boson. 
At the parton level, the EW radiative corrections can enhance the LO cross sections by up to « 30% 
close to the threshold ~ 2m t when the SM Higgs boson is light and reduce the LO cross sections 
with increasing s by up to the same order of magnitude. After convoluting with the parton distribution 
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functions (PDF's), however, they only reduce the LO production cross section a(pp^tiX) at the LHC 
by a few percent [67]: up to 2.5(1.8)% for the following cuts on the transverse momentum pj* and the 
pseudo rapidity rj of the top quark: px > 100(20) GeV and \rj\ < 2.5. 

So far, the studies of loop induced effects of BSM physics in ti production at hadron colliders 
include the following calculations: 

The 0(a) corrections within a general two Higgs doublet model (G2HDM) (=SM with two Higgs 
doublets but without imposing SUSY constraints) to qq-^ti [J7(| [Tl]] and gg—>tt [|7~ij]. In addition to the 
contribution of the W and Z, the gti(qqg) vertex is modified by the virtual presence of five physical 
Higgs bosons which appear in any G2HDM after spontaneous symmetry breaking: H°, h°, A , H^. 
Thus, the G2HDM predictions for ti observables depend on their masses and on two mixing angles, j3 
and q. The G2HDM radiative corrections are especially large for light Higgs bosons and for very small 
(< 1) and very large values of tan/3 due to the enhanced Yukawa-like couplings of the top quark to 
the (virtual) Higgs bosons. Moreover, there is a possible source for large corrections due to a threshold 
effect in the renormalised top quark self-energy, i.e. when m t ~ M H ± + m&. In [[7l]] the s-channel 
Higgs exchange diagrams in the gluon fusion subprocess, gg—>h°,H°—>ti, had been included. For 
this workshop we also considered the gg—>A°^>ti contribution [72]. A study of the s-channel Higgs 



exchange diagrams alone, can be found in [ }73| ] (H°) and [74, 75] (H° and A ). They are of particular 
interest, since they can cause a peak-dip structure in the invariant tt mass distribution for heavy Higgs 
bosons, M H o > 2m t , when interfered with the LO QCD tt production processes. 

The SUSY EW 0(a) corrections within the MSSM to qq^ti [[H], |77|, [7§] and gg^ti [fTI 



79]. In [ [Til ] also the squark loop contribution to the gg^h , H° production process in the s channel 
Higgs exchange diagrams has been taken into account. The SUSY EW corrections comprise the con- 
tributions of the supersymmetric Higgs sector, and the genuine SUSY contributions due to the virtual 
presence of two charginos x 1 * 1 , four neutralinos x°, two top squarks tz,R and two bottom squarks bi,R- 
The MSSM input parameters can be fixed in such a way that the tt observables including MSSM loop 
corrections depend on a relatively small set of parameters [|7T|]: tan /3, M^o , , , <3? 4 ~, fi, Mi, where 
LR mixing is considered only in the top squark sector, parametrized by the mixing angle <E>£. mr and 
m l L = m l R denote the mass of the lighter top squark and the bottom squark, respectively. The effects 
of the supersymmetric Higgs sector tend to be less pronounced than the ones of the G2HDM: since su- 
persymmetry tightly correlates the parameters of the Higgs potential, the freedom to choose that set of 
parameters which yield the maximum effect is rather limited. On the other hand, they can be enhanced 
by the genuine SUSY contribution depending on the choice of the MSSM input parameters. The SUSY 
EW corrections can become large close to the threshold for the top quark decay t—>i + x • They are 
enhanced for very small (< 1) and very large values of tan (3 and when there exists a light top squark 
{m- ti m 100 GeV). 



The SUSY QCD 0(a s ) corrections to qq^tt [{7J| g gjj, g §3p and gg^tt @]. So far, 
there are only results available separately for the qq—>ti (Tevatron) and the gg^tt (LHC) production 
processes. The combination of both is work in progress and will be presented in [|5|]. The SUSY QCD 
contribution describes the modification of the gti(qqg) vertex and the gluon vacuum polarisation due to 
the virtual presence of gluinos and squarks. Thus, additionally to the dependence on squarks masses (and 
on mixing angles if LR mixing is considered) the SUSY QCD corrections introduce a sensitivity of tt 
observables on the gluino mass nig. As expected, the effects are the largest the lighter the gluino and/or 
the squarks. Again, there are possible enhancements due to threshold effects, for instance close to the 

anomalous threshold m? = m? + mf . 

L y ti 

The ti observables under investigation so far comprise the total ti production cross section a t f, 
the invariant ti mass distribution da/dM t i and parity violating asymmetries Alr in the production 
of left and right handed top quark pairs. At present, the numerical discussion is concentrated on the 
impact of BSM quantum effects on ti observables in p { p ^tiX. A parton level Monte Carlo program 
for pp ' ^ti^W + W~bb^(fif-)(fjfj)bb is presently under construction [72]. This will allow a more 
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Table 5: The relative corrections to pp^ttX at the LHC when only including SUSY QCD one-loop corrections [ p4| (with 
p T >20 GeV, \rj\ < 2.5) or only the EW one-loop corrections within the G2HDM and the MSSM §^ (p T > 100 GeV). For 
comparison the SM prediction is also listed. 





SM (M H = 100 GeV) 


G2HDM 


SUSY EW 


SUSY QCD 


\ a NLO _ a LOu a Lp 

\ u tt u tt \l u tt 


2.5% 


<4% 


< 10% 


< 4% 



realistic study of the sensitivity of a variety of kinematical distributions to SUSY quantum signatures in 
the tt production processes, for instance by taking into account detector effects. 

In the following we give an overview of the present status of BSM quantum effects in tt observ- 
ables at the LHC: 

a t i : In Table || we provide the relative corrections for a t i at the LHC for different BSM physics 
scenarios. They reflect the typical maximum size of the radiative corrections within the models under 
consideration. As already mentioned there are possible enhancements due to threshold effects, which 
can yield much larger relative corrections. However, they only arise for very specific choices of the 
MSSM input parameters. The SUSY EW one-loop corrections always reduce the LO production cross 
sections and range from SM values, to up to « —5% for heavy squarks and up to « —20% close to 
mt = mf + m^o. The SUSY QCD one-loop corrections, however, can either reduce or enhance a t j. 
The relative corrections are negative for small m g and increase with decreasing gluino and/or squark 
masses. They change sign when approaching the threshold for real sparticle production and reach a 
maximum at m g 200 GeV of about +2% [84]. Again, very large corrections arise in the vicinity of a 
threshold for real sparticle production, m t = m g + . The SUSY EW and QCD one-loop corrections, 
so far, have only been combined for the qq—>ti production process and numerical results are provided 
for the Tevatron pp collider in [78, 83]. To summarise, apart from exceptional regions in the MSSM 
parameter space, it will be difficult to detect SUSY through loop contributions to the tt production 
rate. If light sparticles exist, they are most likely directly observed first. Then, the comparison of the 
precisely measured top production rate with the MSSM predictions will test the consistency of the model 
under consideration at quantum level and might yield additional information on the parameter space, for 
instance constraints on tan and <I> t ~. 

dcr / dM t t : More promising are the distributions of kinematic variables. Here we will concen- 
trate on the impact of SUSY quantum signatures on the invariant tt mass distribution. Results for the 
effects of EW one-loop corrections within the G2HDM and the MSSM on da/dM t j at the LHC are 
provided in [f7l[]. So far, the impact of the SUSY QCD one-loop contribution on da/M t i has only 
been discussed for the Tevatron pp collider [|8T|], where it turned out that they can significantly change 
the normalisation and distort the shape of da / dM t f. As already mentioned, there is the possibility for 
an interesting peak-dip structure due to a heavy neutral Higgs resonance in gg—>ti within two Higgs 
doublet models. The potential of the LHC for the observation of such resonances has been studied 
in [74, 86]. In Section 3^ the results of an ATLAS analysis of the observability of the H/A^tt chan- 
nel for different luminosities are presented. In Fig. |9] we show preliminary results for the invariant tt 
mass distribution to pp^tt^W + W~bb— >(v e e + )(du)bb at the LHC when including MSSM EW one- 
loop corrections [72]. When M A o > 2m t the gg^>H° , A°—>tt contributions can cause an excess of ti 
events at M t i slightly below M A o , when the Higgs bosons are not too heavy, and a dip in the distribution 
slightly above M t i = M A o. For the choice of MSSM parameters used in Fig. ^ the peak vanishes for 
M A o > 400 GeV and only a deficiency of events survives which decreases rapidly for increasing M A o . 
These effects can be enhanced when the SUSY QCD contributions are taken into account. 

-4lr ; Parity violating asymmetries in the distribution of left and right-handed top quark pairs 
at hadron colliders directly probe the parity non-conserving parts of the non-QCD one-loop corrections 
to the ti production processes within the model under consideration and have been studied at the Teva- 
tron [0 1|, I?), ||, HI Hi and at the LHC Q. In Fig. [jj] we show the left-right asymmetries Alr in 
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Fig. 9: The LO and NLO invariant mass distributions da/dMu (left) and the relative corrections (right) to the reaction 
pp^ti^{v e e + ){du)bb at the LHC with pr(e,q) > 20 GeV, E T {v) > 20 GeV and |f?(e,g)| < 3.2 for different values 
of M A o (tan/3 = 1.6, $ t - = tt/8, m h = 160 GeV, m- bL = 500 GeV, (j, = 120 GeV and M 2 = 3|/i|.). For comparison 
the relative correction when only taking into account the EW SM one-loop contribution is also shown. The CTEQ3LO set of 
PDF's is used and m t = 174 GeV. 
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Fig. 10: The left-right asymmetry Alr in the invariant ti mass distribution to pp—>ttX at the LHC with pr > 100 GeV 

for different values of $j and and for two extreme choices of tan/?: tan/3 = 0.7 (left) and tan/3 = 50 (right) (with 

m h = 90 GeV, M H ± = 110 GeV, fj, = 120 GeV and M 2 = 3|/x|). 

the invariant mass distribution of (longitudinally) polarised top quark pairs in pp—>tL,RtL,RX, induced 
by SM and MSSM EW one-loop corrections []88|]. The parity violating asymmetry within the MSSM 
results from the interplay of the supersymmetric Higgs sector (M H ±) and the genuine SUSY contribu- 
tions (x^iX )- The contribution from the charged Higgs boson can either be enhanced or diminished 
depending on the values of m~ t and Within the G2HDM the loop-induced asymmetries are most 
pronounced for a light charged Higgs boson and very small and very large values of tan j3. At the LHC, 
the G2HDM and MSSM EW one-loop corrections induce asymmetries in the total production rate of left 
and right-handed top quark pairs of up to about 2.5% and 3.2%, respectively, and thus can be consider- 
ably larger than the SM expectation (SM: 1.2%). When the squarks are non-degenerate in mass also the 
SUSY QCD one-loop corrections induce parity violating asymmetries in strong ti production. So far, 
there exist only studies for the Tevatron [|2l 181], |3f|. 



3.5 Measurement of tt production properties 

According to the SM, the top quark decays almost exclusively via t — > Wb. The final state topology 
of tt events then depends on the decay modes of the W bosons. In approximately 65.5% of tt events, 
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both W bosons decay hadronically via W —> jj, or at least one W decays via W — > tv. These events 
are difficult to extract cleanly above the large QCD multi-jet background, and are for the most part 
not considered further. Instead, the analyses presented here concentrate on leptonic tt events, where at 
least one of the W bosons decays via W — > iv (£ = e,p). The lepton plus large E™ iaa , due to the 
escaping neutrino(s), provide a large suppression against multi-jet backgrounds. The leptonic events, 
which account for approximately 34.5% of all tt events, can be subdivided into a "single lepton plus 
jets" sample and a "di-lepton" sample, depending on whether one or both W bosons decay leptonically. 
As discussed below, the selection cuts and background issues are quite different for the various final state 
topologies. 

An important experimental tool for selecting clean top quark samples is the ability to identify b- 
jets. Techniques for 6-tagging, using secondary vertices, semi-leptonic 6-decays, and other characteristics 
of 6-jets, have been extensively studied. Both ATLAS and CMS expect to achieve, for a 6-tagging 
efficiency of 60%, a rejection of at least 100 against prompt jets (i.e. jets containing no long-lived 
particles) at low luminosity. At high luminosity, a rejection factor of around 100 can be obtained with a 
somewhat reduced 6-tagging efficiency of typically 50%. 

All the results presented in this section are obtained using for the signal the PYTHIA Monte Carlo 
program. Most background processes have also been generated with PYTHIA, with the exception of Wbb, 
which has been produced using the HERWIG implementation [89] of the exact massive matrix-element 
calculation. 



3.51 Single lepton plus jets sample 

The single lepton plus jets topology, tt — ► WWbb —> (£u)(jj)bb arises in 2 x 2/9 x 6/9 as 29.6% of 
all tt events. One expects, therefore, production of almost 2.5 million single lepton plus jet events for 
an integrated luminosity of 10 fb _1 , corresponding to one year of LHC running at 10 33 cm~ 2 s . The 
presence of a high px isolated lepton provides an efficient trigger. The lepton and the high value of E™ isa 
give a large suppression of backgrounds from QCD multi-jets and bb production. 

For the single lepton plus jets sample, it is possible to fully reconstruct the final state. The four- 
momentum of the missing neutrino can be reconstructed by setting M v = 0, assigning E T = E™ iaa , and 
calculating p u z , with a quadratic ambiguity, by applying the constraint that M = M\y- 

An analysis by ATLAS [3C] examined a typical set of selection cuts. First, the presence of an 
isolated electron or muon with pr > 20 GeV and \r)\ < 2.5 was required, along with a value of E™ ' aa > 
20 GeV. At least four jets with px > 20 GeV were required, where the jets were reconstructed using a 
fixed cone algorithm with cone size of AR = 0.7. After cuts, the major sources of backgrounds were 
VF+jet production with W — ► Iv decay, and Z+jet events with Z — ► £ + l~ . Potential backgrounds 
from WW, WZ, and ZZ gauge boson pair production have also been considered, but are reduced to a 
negligible level after cuts. 

A clean sample of it events was obtained using ^-tagging. Requiring that at least one of the jets be 
tagged as a 6-jet yielded a selection efficiency (not counting branching ratios) of 33.3%. For an integrated 
luminosity of 10 fb -1 , this would correspond to a signal of 820,000 tt events. The total background, 
dominated by VK+jet production, leads to a signal-to-background ratio (S/B) of 18.6. Tighter cuts can be 
used to select a particularly clean sample. Examples of this will be given in Section |] 



3.52 Di-lepton sample 

Di-lepton events, where each W decays leptonically, provide a particularly clean sample of tt events, 
although the product of branching ratios is small, 2/9 x 2/9 rj 4.9%. With this branching ratio, one 
expects the production of over 400,000 di-lepton events for an integrated luminosity of 10 fb _1 . 

The presence of two high px isolated leptons allows these events to be triggered efficiently. Back- 
grounds arise from Drell-Yan processes associated with jets, Z — > r + r~ associated with jets, VFl^+jets, 
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and bb production. Typical selection criteria [ ft (4 190| ] require two opposite-sign leptons within \rj\ < 2.5, 
with pt > 35 and 25 GeV respectively, and with E™ iss > 40 GeV. For the case of like-flavour leptons 

(e + e~ and 

/i+/0> an additional cut \M U - M z \ > 10 GeV was made on the di-lepton mass to remove 
Z candidates. Requiring, in addition, at least two jets with pr > 25 GeV produced a signal of 80,000 
events for 10 fb _1 , with S/B around 10. Introducing the requirement that at least one jet be tagged as a 
6-jet reduced the signal to about 58,000 events while improving the purity to S/B k, 50. 

3.53 Multi-jet sample 

The largest sample of tt events consists of the topology tt — ► WWbb — > (jj)(jj)bb. The product of 
branching ratios of 6/9x6/9 « 44.4% implies production of 3.7 million multi-jet events for an integrated 
luminosity of 10 fb^ 1 . However, these events suffer from a very large background from QCD multi-jet 
events. In addition, the all-jet final state poses difficulties for triggering. For example, the trigger menus 



examined so far by ATLAS Q30| ] consider multi-jet trigger thresholds only up to four jets, for which a jet 
Et threshold of 55 GeV is applied at low luminosity. Further study is required to determine appropriate 
thresholds for a six-jet topology. 

At the Fermilab Tevatron Collider, both the CDF and DO collaborations have shown that it is pos- 
sible to isolate a tt signal in this channel. The CDF collaboration has obtained a signal significance over 
background of better than three standard deviations |Q] by applying simple selection cuts and relying on 
the high 6-tagging efficiency (~ 46%). To compensate for the less efficient ^-tagging, the DO collab- 



oration has developed a more sophisticated event selection technique [jlCy. Ten kinematic variables to 
separate signal and background were used in a neural network, and the output was combined in a second 
network together with three additional variables designed to best characterise the tt events. 

ATLAS has made a very preliminary investigation [3(J 91 ] of a simple selection and reconstruction 
algorithm for attempting to extract the multi-jet tt signal from the background. Events were selected by 
requiring six or more jets with pr > 15 GeV, and with at least two of them tagged as 6-jets. Jets were 
required to satisfy \rj\ < 3 (\rj\ < 2.5 for &-jet candidates). In addition, the scalar sum of the transverse 
momenta of the jets was required to be greater than 200 GeV. The tt signal efficiency for these cuts was 
19.3%, while only 0.29% of the QCD multi-jet events survived. With this selection, and assuming a 
QCD multi-jet cross-section of 1.4 x 10~ 3 mb for p? (hard process) > 100 GeV, one obtains a signal-to- 
background ratio S/B 1/57. 

Reconstruction of the tt final state proceeded by first selecting di-jet pairs, from among those jets 
not tagged as 6-jets, to form W — ► jj candidates. A Xyv was calculated from the deviations of the two 
Mjj values from the known value of Mvy- The combination which minimised the value of Xw was 
selected, and events with x%/ > 3.5 were rejected. For accepted events, the two W candidates were then 
combined with fo-tagged jets to form top and anti-top quark candidates, and a x\ wa s calculated as the 
deviation from the condition that the top and anti-top masses are equal. Again, the combination with the 
lowest Xt wa s selected, and events with Xt > 7 were rejected. After this reconstruction procedure and 
cuts, the value of S/B improved to 1/8 within the mass window 130-200 GeV. Increasing the p? threshold 
for jets led to some further improvement; for example, requiring pj, > 25 GeV yielded S/B = 1/6. 

The isolation of a top signal can be further improved in a number of ways, such as using a multi- 
variate discriminant based on kinematic variables like aplanarity, sphericity or Ai?(jet-jet), or restricting 
the analysis to a sample of high p l ^ v events. These techniques are undergoing further investigation, but 
it will be very difficult to reliably extract the signal from the background in this channel. In particular, 
the multi-jet rates and topologies suffer from very large uncertainties. 

3.54 Measurement of the tt invariant mass spectrum 

As discussed previously, properties of events provide important probes of both SM and BSM physics. 
For example, a heavy resonance decaying to tt might enhance the cross-section, and might produce a 
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peak in the M# invariant mass spectrum. Deviations from the SM top quark branching ratios, due for 
example to a large rate of t — ► H + b, could lead to an apparent deficit in the ti cross-section measured 
with the assumption that BR(t — > Wb) 1. 

Due to the very large samples of top quarks which will be produced at the LHC, measurements of 
the total cross-section a(tt) will be limited by the uncertainty of the integrated luminosity determination, 
which is currently estimated to be 5%-10%. The cross-section relative to some other hard process, such 
as Z production, should be measured more precisely. 



Concerning differential cross-sections, particular attention has thus far been paid by ATLAS [ pOy 
to measurement of the il% invariant mass spectrum. A number of theoretical models predict the existence 
of heavy resonances which decay to ti. An example within the SM is the Higgs boson, which will decay 
to ti provided the decay is kinematically allowed. However, the strong coupling of the SM Higgs boson 
to the W and Z implies that the branching ratio to ti is never very large. For example, for Mh = 500 
GeV, the SM Higgs natural width would be 63 GeV, and BR(iT -> ti) « 17%. The resulting value of 
ctxBR for H — > ti in the SM is not sufficiently large to see a Higgs peak above the large background 
from continuum ft production. In the case of MSSM, however, if Mjj,a > 2mt, then BR(H/A — > ti) 
« 100% for tan j3 « 1. For the case of scalar or pseudo-scalar Higgs resonances, it has been pointed 
out J73| , f74| ] that interference can occur between the amplitude for the production of the resonance via 
gg — > H/A — > ti and the usual gluon fusion process gg — ► ti. The interference effects become stronger 
as the Higgs' mass and width increase, severely complicating attempts to extract a resonance signal. 

The possible existence of heavy resonances decaying to ti arises in technicolor models [^] as well 
as other models of strong EW symmetry breaking ||93|]. Recent variants of technicolor theories, such as 
Topcolor J94]], posit new interactions which are specifically associated with the top quark, and could give 
rise to heavy particles decaying to ti. Since ti production at the LHC is dominated by gg fusion, colour 
octet resonances ("colourons") could also be produced [|95|]. 



Because of the large variety of models and their parameters, ATLAS performed a study [ |30j , |96J 
of the sensitivity to a "generic" narrow resonance decaying to ti. Events of the single lepton plus jets 
topology ti — ► WWbb — > (£u)(jj)bb were selected by requiring E™ tss > 20 GeV, and the presence 
of an isolated electron or muon with px > 20 GeV and \r)\ < 2.5. In addition, it was required that 
there were between four and ten jets, each with pr > 20 GeV and \rj\ < 3.2. At least one of the jets 
was required to be tagged as a 6-jet. After these cuts, the background to the ti resonance search was 
dominated by continuum ti production. 

The momentum of the neutrino was reconstructed, as described previously, by setting M v = 0, 
assigning E T = E™' ss , and calculating p v z (with a quadratic ambiguity) by applying the constraint that 
Mi v = My/- The hadronic W — > jj decay was reconstructed by selecting pairs of jets from among 
those not tagged as 6-jets. In cases where there were at least two 6-tagged jets, candidates for t — ► Wb 
were formed by combining the W candidates with each fc-jet. In events with only a single 6-tagged jet, 
this was assigned as one of the 6-quarks and each of the still unassigned jets was then considered as a 
candidate for the other 5-quark. 

Among the many different possible jet-parton assignments, the combination was chosen that min- 
imised the following x 2 : 

X 2 = {M jjb - m t ) 2 /a 2 (M jjb ) + {M M - m t f j a 2 {M lvh ) + (M# - M w f /^{Mjj) 

Events were rejected if either or Mjjb disagreed with the known value of mt by more than 30 GeV. 

For events passing the reconstruction procedure, the measured energies were rescaled, according 
to their resolution, to give the correct values of M\y and mt for the appropriate combinations. This 
procedure improved the resolution of the mass reconstruction of the ti pair to a{M t i) /M t i 6.6%. As 



an example, Fig. 1 1 shows the reconstructed M t j distribution for a narrow resonance of mass 1600 GeV. 



The width of the Gaussian core is well described by the resolution function described above. The size 
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Fig. 1 1 : Measured ti invariant mass distribution for reconstruction of a narrow resonance of mass 1 600 GeV decaying to ti . 
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Fig. 12: Value of a x BR required for a 5a discovery potential for a narrow resonance decaying to ti, as a function of M t t, and 
for an integrated luminosity of either 30 or 300 fb -1 . 



of the tails, which are dominated by incorrect jet-parton assignments, is such that approximately 65% of 
the events are contained within ±2<t of the peak. 

The reconstruction efficiency, not including branching ratios, for ti — > WWbb — ► (£u)(jj)bb was 
about 20% for a resonance of mass 400 GeV, decreasing gradually to about 15% for M t j = 2 TeV. 

For a narrow resonance X decaying to ti, Fig. 12 shows the required crx BR(X — > ti) for dis- 
covery of the resonance. The criterion used to define the discovery potential was observation within a 
±2c7 mass window of a signal above the ti continuum background, where the required signal must have 
a statistical significance of at least 5a and must contain at least ten events. Results are shown versus Mx 
for integrated luminosities of 30 fb -1 and 300 fb -1 . For example, with 30 fb -1 , a 500 GeV resonance 
could be discovered provided its crx BR is at least 2560 fb. This value decreases to 830 fb for Mx = 
1 TeV, and to 160 fb for Mx = 2 TeV. The corresponding values for an integrated luminosity of 300 fb -1 
are 835 fb, 265 fb, and 50 fb for resonances masses Mx = 500 GeV, 1 TeV, and 2 TeV, respectively. 

Once predictions from models exist for the mass, natural width, and crx BR for a specific reso- 
nance, the results in Fig. 12 can be used to determine the sensitivity and discovery potential for those 
models. As discussed above, for the case of scalar or pseudo-scalar Higgs resonances, extra care must be 
taken due to possible interference effects. While such effects are small for the case of a narrow resonance, 
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Fig. 1 3 : For various integrated luminosities, 5a discovery contours in the MSSM (M A , tan ff) plane for the channel H, A — » ti. 



they can be significant once the finite widths of heavy resonances are taken into account. For example, 
ATLAS has performed an analysis [|(| |7|] of the decays H/A —* tt in MSSM with tan/3 = 1.5 and 
Mh,a > 2m$- Assuming the ti continuum background is well known, a combined H + A signal would 
be visible for Higgs masses in the range of about 370 - 450 GeV. However, the interference effects pro- 
duce an effective suppression of the combined H + A production rates of about 30% for Mh,a = 370 
GeV, increasing to 70% for masses of 450 GeV, essentially eliminating the possibility to extract a sig- 
nal for higher Higgs masses, and thereby severely limiting the MSSM parameter space for which this 
channel has discovery potential (see Fig. 13). 



4. TOP QUARK MASS[] 

As discussed in Section |2/J one of the main motivations for top physics at the LHC is an accurate 
measurement of the top mass. Currently the best Tevatron single-experiment results on m t are obtained 
with the lepton plus jets final states. These yield: m t = 175.9 ± 4.8 (stat.) ± 5.3 (syst.) (CDF) |@] 
and 173.3 ± 5.6 (stat.) ± 5.5 (syst.) (D0) [Qj. The systematic errors in both measurements are largely 
dominated by the uncertainty on the jet energy scale which amounts to 4.4 GeV and 4 GeV for CDF 
and D0, respectively. On the other hand, the systematic errors in the di-lepton channels are somewhat 
less, but the statistical errors are significantly larger, by a factor of > 2, as compared to the lepton 
plus jets final states. Future runs of the Tevatron with an about 20-fold increase in statistics promise a 
measurement of the top mass with an accuracy of up to ~ 3 GeV [ 98 ] ; in the lepton plus jets channel 
the error is dominated by the systematics while in the di-lepton channels the limiting factor is still the 
statistics. 

Several studies of the accuracy which can be expected with the LHC experiments have been per- 
formed in the past ||99|]. It is interesting to see whether one can use the large statistics available after a 
few years of high-luminosity running to push the precision further. In particular, it is interesting to study 
the ultimate accuracy achievable at a hadronic collider, and the factors that limit this accuracy. 

In the following subsections, we begin with general remarks on the top quark mass and a very 
brief review of the present status of the theoretical understanding of top quark mass measurement in the 
threshold scan at a future e + e~ collider. We then present the results of a recent studies of top mass 
reconstruction at the LHC. The techniques used include the study of the lepton plus jets final states 
(inclusive, as well as limited to high-p T top quarks), di-lepton final states (using the di-leptons from the 
leptonic decay of both W's, as well as samples where the isolated W lepton is paired with a non-isolated 
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lepton from the decay of the companion b hadron). A very promising analysis using the J ftp from the 
b hadron decay paired with the lepton from the leptonic decay of the W is discussed at the end. The 
conclusions of these studies indicate that an accuracy of 2 GeV should be achievable with the statistics 
available after only 1 year of running at low luminosity. An accuracy of 1 GeV accuracy could be 
achieved after the high luminosity phase. 



4.1 General remarks and the top mass measurement in e + e annihilation 

Although one speaks of "the" top quark mass, one should keep in mind that the concept of quark mass 
is convention-dependent. The top quark pole mass definition is often implicit, but in a confining theory 
it can be useful to choose another convention. This is true even for top quarks when one discusses mass 
measurements with an accuracy of order of or below the strong interaction scale. Since different mass 
conventions can differ by 10 GeV (see Section 2. 1 ), the question arises which mass is actually determined 
to an accuracy of 1-2 GeV by a particular measurement. 

The simple answer is that a particular measurement determines those mass parameters accurately 
in terms of which uncalculated higher order corrections to the matrix elements of the process are small. 
This in turn may depend on the accuracy one aims at and the order to which the process has already 
been calculated. To clarify these statements we briefly discuss the top quark mass measurement at a high 
energy e + e~ collider. 

"The" top quark mass can be measured in e + e~ collisions by recontructing top quark decay prod- 
ucts in much the same way as at the LHC. In addition, there exists the unique possibility of determining 
the mass in pair production near threshold. This is considered to be the most accurate method [ 10C ] and 
it appears that an uncertainty of 5m t ~ 0.15 GeV can be achieved for the top quark M S mass with the 
presently available theoretical input [101]. This is a factor two improvement compared to the accuracy 
that could be achieved with the same theoretical input if the cross section were parametrised in terms 
of the top quark pole mass. The fundamental reason for this difference is the fact that the concept of 
a quark pole mass is intriniscally ambiguous by an amount of order Aqcd [ 102 ] and this conclusion 
remains valid even if the quark decays on average before hadronisation [103]. In the context of pertur- 
bation theory this ambiguity translates into sizeable higher order corrections to the matrix elements of a 
given process renormalized in the pole mass scheme. This makes it preferable to choose another mass 
convention if large corrections disappear in this way as is the case for the total cross section in e + e~ an- 
nihilation, because the total cross section is less affected by non-perturbative effects than the pole mass 
itself. Note, however, that despite this preference the position of the threshold is closer to twice the pole 
mass than twice the MS mass, hence a leading order calculation determines the pole mass more naturally. 
It is possible to introduce intermediate mass renormalizations that are better defined than the pole mass 
and yet adequate to physical processes in which top quarks are close to mass shell [101, 104]. The con- 
clusion that the top quark pole mass is disfavoured is based on the existence of such mass redefinitions 
and the existence of accurate theoretical calculations. 

The situation with mass determinations at the LHC appears much more complicated, since the 
mass reconstruction is to a large extent an experimental procedure based on leading order theoretical 
calculations, which are not sensitive to mass renormalization at all. Furthermore the concept of invari- 
ant mass of a top quark decay system is prone to "large" non-perturbative corrections of relative order 
AqcdA^, because the loss or gain of a soft particle changes the invariant mass squared by an amount 
of order mt Aqcd- The parametric magnitude of non-perturbative corrections is of the same order of 
magnitude as for the top quark pole mass itself and cannot be decreased by choosing another mass renor- 
malization prescription. For this reason, top mass measurements based on reconstructing m t from the 
invariant mass of the decay products of a single top quark should be considered as measurements of the 
top quark pole mass. From the remarks above it follows that there is a limitation of principle on the ac- 
curacy of such measurements. However, under LHC conditions the experimental systematic uncertainty 
discussed later in this section is the limiting factor in practice. A potential exception is the measure- 
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Table 6: Efficiencies (in percent) for the inclusive ti single lepton plus jets signal and for background processes, as a function 
of the selection cuts applied. No branching ratios are included in the numbers. The last column gives the equivalent number of 
events for an integrated luminosity of 10 fb _1 , and the signal-to-background ratio. 





p l T > 20 GeV 


as before, 


as before, 


events, 


Process 


E™' ss > 20 GeV 


plus N jet > 4 


plus N b _ jet > 2 


per 10 fb _1 


ti signal 


64.7 


21.2 


5.0 


126000 


W + jets 


47.9 


0.1 


0.002 


1658 


Z + jets 


15.0 


0.05 


0.002 


232 


WW 


53.6 


0.5 


0.006 


10 


wz 


53.8 


0.5 


0.02 


8 


zz 


2.8 


0.04 


0.008 


14 


Total background 








1922 


S/B 


65 



ment of nit in the decay mode £J/ipX discussed at the end of this section, since the systematic error 
is estimated to be below 1 GeV and since the systematic error is to a large extent theoretical. It would 
be interesting to investigate non-perturbative power corrections and principle obstructions to an accurate 
mass measurement for this process. This analysis has however not yet been carried out in any detail, 
comparable to the threshold scan in e + e~ annihilation. 



4.2 m t in the lepton plus jets channel. Inclusive sample 

The inclusive lepton plus jets channel provides a large and clean sample of top quarks for mass recon- 
struction. Considering only electrons and muons, the branching ratio of this channel is 29.6%. Therefore, 
one can expect more than 2 millions events for one year of running at low luminosity. ATLAS performed 
an analysis in that channel using events generated using PYTHIA [52] and the ATLAS detector fast simu- 
lation package ATLFAST [105]. The top mass is determined using the hadronic part of the decay, as the 
invariant mass of the three jets coming from the same top: m t = rrijjb. The leptonic top decay is used to 
tag the event with the presence of a high p T lepton and large E™ lss . For the background processes, the 
HERWIG [Ml, 89] generator was used for the background process Wbb. 

The following background processes have been considered: bb, W + jets with W — > £v, Z + jets 
with Z -> £ + £~, WW with one W — ► tv and the other W -> qq, WZ with W -> £v and Z -> qq, 
ZZ with one Z — ► £ + £~ and Z — ► qq, and Wbb with W — ► £v. Events are selected by requiring an 
isolated lepton with p T > 20 GeV and |ry|< 2.5, E™ iss > 20 GeV, and four jets with p T > 40 GeV and 
\w\< 2.5, of which two of them were required to be tagged as 6-jets. Jets were reconstructed using a 
fixed cone algorithm with AR= 0.4. Although at production level the signal over background is very 
unfavourable, after the selection cuts and for an integrated luminosity of 10 fb _1 , 126000 signal events 
and 1922 background events were kept, yielding a value of S/B = 65 (see Table |J). 

The reconstruction of the decay W — > jj is first performed. The invariant mass nijj of all the 
combinations of jets (with p T > 40 GeV and \rj\< 2.5) that were not tagged as 6-jets is computed and 
the jet pair with an invariant mass closest to mw is selected as the W candidate. Fig. 14 represents the 
invariant mass distribution of the selected jet pairs. The reconstructed W mass is consistent with the 
generated value, the mass resolution being 7.8 GeV. Within a window of ±20 GeV around the W mass, 
the purity (P) and the overall efficiency (E) of the W reconstruction are respectively P=67% and E=1.7%. 
Additional pair association criteria, such as requiring the leading jet to be part of the combination, did 
not improve significantly the purity and have not been considered further in the analysis. W candidates, 
retained if | mjj — Mw | < 20 GeV, have then to be associated with one 6-tagged jet to reconstruct the 
decay t — ► Wb. To reconstruct the right combination, some association criteria have been tried, such as 
choosing the 6-jet furthest from the isolated lepton, the 6-jet closest to the reconstructed W, and choosing 
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Fig. 14: Left: invariant mass distribution of the selected jj pairs. Right: invariant mass distribution of the selected jjb com- 
bination. Both distributions are normalised to an integrated luminosity of lOfb -1 . The shaded area shows the combinatorial 
background. 
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Fig. 15: Left: dependence of the reconstructed top mass on the generated value. Right: dependence of the reconstructed top 
mass on the transverse momentum(p T ) of the reconstructed top. 



the jjb combination having the highest p T for the reconstructed top. These various methods gave similar 



results. Fig. |14| presents the invariant mass distribution of the reconstructed top when the jjb combi- 
nation having the highest p T has been used as association criteria. No Myy constraint is applied for the 
light quark jets. For an integrated luminosity of 10fb -1 , the total number of reconstructed top is 32000 
events, of which 30000 are within a window of ±35 GeV around the generated top mass mt= 175 GeV 
The total number of combinatorial events is 34000, of which 14000 are within the mass window. The 
number of background events coming from other processes is negligible. The rrijjb distribution fitted 
by a Gaussian plus a third order polynomial yields a top mass consistent with the generated value of 
175 GeV and a top mass resolution of 11.9 GeV The resulting statistical uncertainty for an integrated 
luminosity of lOftr 1 is 5m t = 0.070 GeV 

The dependence of the top reconstruction algorithm on the top mass has been checked using 
several samples of tt events generated with different values of m ( ranging from 160 to 190 GeV The 



results, shown in Fig. [15J, demonstrate a linear dependence of the reconstructed top mass on the generated 
value: the data points are fitted to a linear function with x 2 /ndf = 6.7/8. The stability of the mass value 
as a function of the transverse momentum of the reconstructed top (p T (top)) was also checked. As shown 
in Fig. |l5j, no significant p T (top) dependence is observed: the data points are fitted to a constant with 
X 2 /ndf = 6.25/5. For more details of this analysis, see [|106p. 



25 



> 
a 

3t 



a 

> 



300 



200 



LOO 



i 1 

J 


a = 11.4 GeV 




L 


iii Ma ■ i 





> 

^ 200 



a 

> 



150 — 



100 



50 



100 



200 



300 



m jjb (GeV) 



400 




m jjb (GeV) 



400 



Fig. 16: Invariant jjb mass distributions. Left: from fast simulation. Right: from full simulation. 

The results presented above, obtained with a fast simulation package, have been cross-checked 
with 30000 events passed through the ATLAS GEANT-based full simulation package [107]. In full 
simulation, in order to save computing time, events have been generated under restrictive conditions at 
the generator level. The comparison is done by using the same generated events which have been passed 
through both the fast and full simulation packages. The results, in terms of purity, efficiency and mass 
resolutions show a reasonable agreement between fast and full simulation. In addition, as it is shown in 
Fig. 16, the shape and amount of the combinatorial background for the rrijjb distributions are in good 



agreement between the two types of simulations. 

It has to be noted that for this analysis as well as for the other top mass reconstruction studies 
performed within ATLAS, the jets were calibrated using the ratio p T (parton)/p T (jet) obtained from Monte 
Carlo samples of di-jet events or H — > bb with mu = 100 GeV. In that aspect this calibration does not 



include all possible detector effects and corrections. More details can be found in Chapter 20 of Q30| ] and 
in Appendix 



4.3 m t in the lepton plus jets channel. High px sample 

An interesting possibility at the LHC, thanks to the large ti production rate, is the use of special sub- 
samples, such as events where the top and anti-top quarks have high p T . In this case, they are produced 
back-to-back in the lab-frame, and the daughters from the two top decays will appear in distinct "hemi- 
spheres" of the detector. This topology would greatly reduce the combinatorial background as well as the 
backgrounds from other processes. Furthermore, the higher average energy of the jets to be reconstructed 
should reduce the sensitivity to systematic effects due to the jet energy calibration and to effects of gluon 
radiation. However, in this case a competing effect appeals which can limit the resulting precision: as the 
top p T increases, the jet overlapping probability increases as well, which again affects the jet calibration. 
ATLAS performed a preliminary study of this possibility using two different reconstruction methods: 

• in the first one an analysis similar to the inclusive case is done, with m t being reconstructed from 
the three jets in the one hemisphere (mt=mjjb); 

• in the second one, m t is reconstructed summing up the energies in the calorimeter towers in a large 
cone around the top direction. 

In the following paragraphs, highlights of these analyses are discussed. 



4.31 Jet Analysis 



High p T tt events were generated using PYTHIA 5.7 Q52[ ] with a p T cut on the hard scattering process 
above 200 GeV. The expected cross-section in this case is about 120 pb, or about 14.5% of the total 
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Fig. 17: Left: invariant mass distribution of the selected di-jet combinations for the high p T (top) sample. Right: invariant 
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tt production cross-section. The selection cuts required the presence of an isolated lepton with p T > 
30 GeV and \rj\ < 2.5, and E™ iss > 30 GeV. The total transverse energy of the event was required to be 
greater than 450 GeV. Jets were reconstructed using a cone algorithm with radius Ai?=0.4. The plane 
perpendicular to the direction of the isolated lepton was used to divide the detector into two hemispheres. 
Considering only jets with p T > 40 GeV and \t]\< 2.5, the cuts required one 6-tagged jet in the same 
hemisphere as the lepton, and three jets, one of which was 6-tagged, in the opposite hemisphere. Di-jet 
candidates for the W — ► jj decay were selected among the non-6-tagged jets in the hemisphere opposite 
to the lepton. The resultant mjj invariant mass distribution is shown in Fig. 17 (left). Fitting the six bins 



around the peak of the mass distribution with a Gaussian, yielded a W mass consistent with the generated 
value, and a rrijj resolution of 7 GeV, in good agreement with that obtained for the inclusive sample. 
Di-jets with 40 GeV<mjj<120 GeV were then combined with the 6-tagged jet from the hemisphere 
opposite to the lepton to form t — > jjb candidates. Finally, the high p T (top) requirement was imposed by 
requiring p T (jjb)> 250 GeV. With these cuts, the overall signal efficiency was 1.7%, and the background 
from sources other than tt was reduced to a negligible level. The invariant mass distribution of the 



accepted jjb combinations is shown in Fig. |17| (right). Fitting the six bins around the peak of the mass 
distribution with a Gaussian, yielded a top mass consistent with the generated value of 175 GeV, and a 
mjjb mass resolution of 11.8 GeV. For an integrated luminosity of 10 fb _1 , a sample of 6300 events 
would be collected in ATLAS, leading to a statistical error of 5m 4 (stat.) = ±0.25 GeV, which remains 
well below the systematic uncertainty. As in the case of the inclusive sample, no strong p T dependence 
was observed and the reconstructed mass depends linearly on the Monte Carlo input value. 



4.32 Using a large calorimeter cluster 

For sufficiently high p T (top) values, the jets from the top decay are close to each other with a large 
possibility of overlap. In such a case it might be possible to reconstruct the top mass by collecting all the 
energy deposited in the calorimeter in a large cone around the top quark direction. Such a technique has 
the potential to reduce the systematic errors, since it is less sensitive to the calibration of jets and to the 
intrinsic complexities of effects due to leakage outside the smaller cones, energy sharing between jets, 
etc. Some results from a preliminary investigation of the potential of this technique are discussed here. 
More details of the analysis can be found in [30, 108 ]. 

Similar event selection criteria as in the previous case were used: an isolated lepton with p T > 
20 GeV and \rj\ < 2.5, E™ lsa > 20 GeV, one 6-tagged jet (with A.R=0.4 and p T > 20 GeV) in the lepton 
hemisphere, and at least 3 jets in the hemisphere opposite to the lepton (Ai?=0.2, p T > 20 GeV) with 
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one of them 6-tagged. For the accepted events, the two highest p T non-6-tagged jets were combined 
with the highest p T 6-jet candidate in the hemisphere opposite to the lepton to form candidates for the 
jjb hadronic top decay. The selected jjb combination was required to have p T > 150 GeV and \rj\< 2.5. 
With these selection criteria, about 13000 events would be expected in the mass window from 145 to 
200 GeV, with a purity of 90%, for an integrated luminosity of 10fb _1 . The reconstructed invariant mass 



of the jjb combination is shown in Fig. 18 (left). The direction of the top quark was then determined 



from the jet momenta. Figure [jj (right) shows the distance AR in (77, <fi) space between the reconstructed 
and the true top direction at the parton level, demonstrating good agreement. 

A large cone of radius AR was then drawn around the top quark direction, and the top mass was 
determined by adding the energies of all calorimeter "towers" within the cone. A calorimeter tower 
has a size of Srj x 5(j>= 0.1x0.1, combining the information of both the EM and hadronic calorimeters. 
The invariant mass spectrum is shown in Fig. 19 (left) for a cone size AR= 1.3, and exhibits a clean 



peak at the top quark mass. The fitted value of the reconstructed top mass is shown in Fig. |19| (right), 
where it displays a strong dependence on the cone size. If initial (ISR) and final (FSR) state radiation in 
PYTHIA are turned off, the fitted mass remains constant (to within 2%), independently of cone size. 

The large dependence of the reconstructed top mass on the cone size can be attributed to the 
underlying event (UE) contribution. A method was developed to evaluate and subtract the underlying 
event contribution using the calorimeter towers not associated with the products of the top quark decay. 
The UE contribution was calculated as the average E T deposited per calorimeter tower, averaged over 
those towers which were far away from the reconstructed jets of the event. As expected, the average 
E T per calorimeter tower increases as more activity is added, especially in the case of ISR. However, 
only a rather small dependence is observed on the radius AR used to isolate the towers associated with 
the hard scattering process. The resulting value of the reconstructed mass (jn cone ), with and without 



UE subtraction, is also shown in Fig. 19 (right) as a function of the cone radius. As can be seen, after 
the UE subtraction, the reconstructed top mass is independent of the cone size used. As a cross-check, 
the mean E T per cell subtracted was varied by ±10% and the top mass recalculated in each case. As 
shown superimposed on Fig. |^ (right), these "miscalibrations" lead to a re-emergence of a dependence 
of rat on the cone size. While the prescription for the UE subtraction does lead to a top mass which is 
independent of the cone size, it should be noted that the reconstructed mass is about 15 GeV (or 8.6%) 
below the nominal value, mt= 175 GeV, implying that a rather large correction is needed. 
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Fig. 19: Left: reconstructed t jjb mass spectrum obtained using cells in a single cone of size AR= 1.3, normalised to an 
integrated luminosity of 10 fb -1 . Right: the fitted top mass using cells in a single cone, before and after the underlying event 
(UE) subtraction and as a function of the cone size. 



To investigate if this correction can be extracted from the data without relying on Monte Carlo 
simulations, the same procedure was applied to a sample of W+ jet events generated with a range of 
p T comparable to that of the top sample. The W was forced to decay hadronically into jets. The UE 
contribution was estimated with the same algorithm as described above. The results agreed within 1% 
with the values determined for the high p T (top) sample. As in the case of the top events, the reconstructed 
W mass after UE subtraction is independent of the cone size. The average value of rrijj after the UE 
subtraction is about 8.5 GeV (or 10.6%) below the nominal value of mw- The fractional error on rrijj, 
as measured with the P^+jet sample, was used as a correction factor to m cone in the high p T (top) sample. 
For a cone of radius AR= 1.3, the top mass after UE subtraction increases from 159.9 GeV to 176.0 GeV 
after rescaling. The rescaled values of m cone are about 1% higher than the generated top mass. This 
over-correction of m t using the value of mw measured with the same method, is mainly due to ISR 
contributions. If ISR is switched off, the rescaling procedure works to better than 1%. 

4.4 Systematic uncertainties on the measurement of m t in the single lepton plus jets channel 

For the analyses presented above within ATLAS, a number of sources of systematic error have been 
studied using samples of events generated with PYTHIA and simulated mainly with the fast detector 
package ATLFAST, but also using a relatively large number of fully simulated events in order to cross- 
check some of the results. The results of these studies are summarised in Fig. ^ and discussed below. 

Jet energy scale: The measurement of mt via reconstruction of t — > jjb relies on a precise knowl- 
edge of the energy calibration for both light quark jets and 6-jets. The jet energy scale depends on a 
variety of detector and physics effects, including non-linearities in the calorimeter response, energy lost 
outside the jet cone (due, for example, to energy swept away by the magnetic field or to gluon radiation 
at large angles with respect to the original parton), energy losses due to detector effects (cracks, leak- 
age, etc.), and "noise" due to the underlying event. Preliminary studies done in ATLAS indicate that 
a jet energy scale calibration at the level of 1% for both light quark and 6-jets would be feasible at the 
LHC (see discussion in the Appendix |A|). In the case of the mt reconstructed from the invariant mass 
of the three jets (mjjb) the 6-jet energy scale enters directly in the measurement and therefore it must be 
calibrated from other sources, while the energy of the two light quark jets can be calibrated event-by- 
event using the W mass constraint. This would work quite well at least for the inclusive sample, where 
the jets are well separated. In the high p T case, energy sharing algorithms and corrections for the two 
jets are needed, and therefore in order to be conservative we assume in the following that no such an 
event-by-event correction can be made. To estimate the effect of an absolute jet energy scale uncertainty, 
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Fig. 20: Summary of systematic errors in the mt measurement. Left: the observed mass shifts for different effects studied. The 
dashed lines indicates a 1% mass window. Right: the quoted error in the mt measurement. 

different "miscalibration" coefficients were applied to the measured jet energies. A linear dependence 
was observed. 

b-quark fragmentation: The fraction of the original 6-quark momentum which will appear as vis- 
ible energy in the reconstruction cone of the corresponding 6-jet depends on the fragmentation function 
of the 6-quark. This function is usually parametrised in PYTHIA in terms of one variable, e&, using the 
Peterson fragmentation function Q1Q9Q. To estimate the systematic error in m t , the "default" value for 



e b (=-0.006) was varied within its experimental uncertainty (0.0025) [ |1 lOj , |19| ] and the difference in the 
reconstructed mt was taken as the systematic error 5mt- 

Initial and final state radiation: The presence of ISR or FSR can impact the measurement of mt. 
To estimate the systematic error due to these, data samples were generated where ISR or FSR in the 
PYTHIA generator were switched off. In the case of FSR, a large mass shift was observed for a jet cone of 
Ai?=0.4. This is reduced as expected when a larger cone is used. Clearly this case is rather pessimistic 
since the knowledge in both ISR and FSR is typically at the level of 10%. Therefore as a conservative 
estimate of the resultant systematic errors in m t , 20% of the mass shifts were used. 

An alternative approach uses the measured jet multiplicity to search, event-by-event, for the pres- 
ence of hard gluon radiation. Following the convention for this approach adopted at the Tevatron [||, §], 
the mass shift would be defined not by comparing events with radiation switched on and events with 
radiation switched off, but by the difference, Am t , between the value of m t determined from events 
with exactly four jets and that determined from events with more than four jets. The systematic error 
due to effects of initial and final radiation would then be considered as 5mt= Am t /Vl2. Such a calcu- 
lation would yield systematic errors of approximately 0.4-1.1 GeV, smaller than the more conservative 
approach adopted here. 

Background: Uncertainties in the size and shape of the background, which is dominated by "wrong 
combinations" in ti events, can affect the top mass reconstruction. The resultant systematic uncertainty 
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on rrit was estimated by varying the assumptions on the background shape in the fitting procedure. Fits of 
the rrijjb distribution were performed assuming a Gaussian shape for the signal and either a polynomial 
or a threshold function for the background. Varying the background function resulted in a systematic 
error on nit of 0.2 GeV. The structure of the UE can affect the top mass reconstruction. However, as 
discussed above, it is possible to estimate and correct for this effect using data. Given the large statistics 
available at the LHC, it is assumed that the residual uncertainty from the underlying event will be small 
compared to the other errors (note that the UE denotes here a minimum bias event, since the impact of 
ISR has already been accounted for). 

For the particular case of the m t reconstructed using a large calorimeter cluster, similar procedures 
were adopted to estimate the the systematic errors. It is important to notice that, as expected, the use of a 
large cone substantially reduces the effects of FSR and fr-quark fragmentation, each of which gives rise 
to a systematic error of 0. 1 GeV. The uncertainty arising from ISR, which can affect the determination 
of the UE subtraction, is about 0. 1 GeV as well. However, the main uncertainty in this technique comes 
from the calibration procedure. The calibration with the W+ jet sample produces a value of m t which 
is about 1 % above the generated value. Furthermore, the W — ► jj events would suffer from background 
from QCD multi-jet events. On-going studies suggest that one could calibrate using W — ► jj decays 
from the high p T (top) events themselves, selecting those events in which the 6-tagged jet is far away 
from the other two jets of the W decay and then reconstructing the W — ► jj decay using a single cone 
of size AR= 0.8. Further study is required to reliably estimate the potential of this calibration procedure, 
and therefore a conservative systematic uncertainty of 1 % is assigned to it. 

4.5 m t in the di-lepton channel 

Di-Lepton events can provide a measurement of the top quark mass complementary to that obtained 
from the single lepton plus jets mode. The signature of a di-lepton event consists of two isolated high 
p T leptons, high E™ lss due to the neutrinos, and two jets from the 6-quarks. The measurement of m t using 
di-lepton events is not a direct measurement as in the previous case but it relies on the relation between 
the kinematical distributions of the top decay products and mt, and on how they can be reproduced 
by the Monte Carlo simulation. About 400000 di-lepton ti events are expected to be produced in a 
data sample corresponding to an integrated luminosity of 10 fb _1 . Backgrounds arise from Drell-Yan 
processes associated with jets, Z —> tt associated with jets, WW+ jets and bb production. 

Of the many possible kinematic variables which could be studied, ATLAS performed a preliminary 
study using: the mass m^ of the lepton+6-jet system, the energy of the two highest E T jets, and the 
mass mu of the di-lepton system formed with both leptons originating from the same top decay (i.e. 
t — > Ivb followed by b —> ivc). The event selection criteria required two opposite-sign leptons within 
\r]\ < 2.5, with p T > 35 and 25 GeV respectively, and with E™ iss > 40 GeV. Two jets with p T > 25 GeV 
were required in addition. After the selection cuts, 80000 signal events survived, with S/B around 10. 

4.51 Top mass measurement using mil, 

In this analysis, the value of m t was estimated using the expression: 

mt 2 = M^ + 2(m 2 eb }/[l-(cose a }} (11) 

Here, (m| fe ) is the squared mean invariant mass of the lepton and t»-jet from the same top decay. The 
mean value of (cos 6 a), the angle between the lepton and the 6-jet in the W rest frame, can be regarded 
as an input parameter to be taken from Monte Carlo. To obtain a very clean sample, the two highest 
p T jets were required to be tagged as 6-jets, leaving a total of about 15200 signal events per 10 fb -1 . 
One cannot determine, in general, which lepton should be paired with which fe-jet. The pairing which 
gave the smaller value of (m 2 eb ) was chosen, and checking the parton-level information showed that this 
criterion selected the correct pairing in 85% of the cases, for a generated top mass of 175 GeV. The mean 
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Fig. 21: Template distributions of the total energy of the two leading jets in ti events for top quark masses of 165 and 175 GeV. 
The two distributions are normalised to the same area. 

value (m 2 eb ) was measured for samples generated with different input top masses m, and then mt was 
calculated from the expression above. For an integrated luminosity of 10 fb -1 , the expected statistical 
uncertainty on m t using this method is ±0.9 GeV. Major sources of systematics include uncertainty on 
the 6-quark fragmentation function, which produces a systematic error on mt of 0.7 GeV if defined as 
described in Section fO] . Systematic errors due to the effects of FSR and ISR together are about 1 GeV, 
while those due to varying the jet energy scale by 1% are 0.6 GeV. Further studies are required to estimate 
the uncertainties due to the reliance upon the Monte Carlo modelling of the ti kinematics. 



4.52 Top mass measurement using the energy of the two leading jets 

Increased sensitivity could be obtained with a technique which utilises not only the mean, but also the 
shape of the kinematic distribution. As an example, a study has been made of the sensitivity to mt ob- 
tained by comparing to "template" distributions the energy of the two highest-i? T jets. The template 
distributions were made by generating PYTHIA samples of ti events with different values of m t in the 
range 160-190 GeV, in steps of 5 GeV. Figure 21 shows, as an example, the templates obtained for 
mt = 165 GeV and 175 GeV. For each possible top mass value m, a x 2 ( m ) was obtained by compar- 
ing the kinematical distribution of the simulated data with the templates of mass m. The best value for 
the mass was the value which, for the "data" set, generated with m t = 175 GeV, gave the minimum y 2 - 
For an integrated luminosity of 10 fb^ 1 , the expected statistical sensitivity on m t corresponds to about 
±0.4 GeV. Varying the calorimeter jet energy scale by 1% produced a systematic error on m t of 1.5 GeV. 
Other sources of systematic error result from the dependence of the method on the Monte Carlo mod- 
elling of the ti kinematics, and require further study. As an example, changing the choice of the structure 
functions used in the Monte Carlo simulation (for example, from CTEQ2L to CTEQ2M or EHQL1) led 
to differences in the top mass of ±0.7 GeV. 



4.53 Top mass measurement using mu in tri-lepton events 

The invariant mass distribution of the two leptons from the same top quark decay {i.e. t — ► Ivb followed 
by b — ► ivc) is quite sensitive to m t . It has been shown that the mass distribution of lepton pairs from the 
same top quark decay is much less sensitive to the top quark transverse momentum distribution than that 
of lepton pairs from different top quarks Signal events are expected to contain two leptons from 
the decay of the W bosons produced directly in the top and anti-top quark decays, and one lepton from 
the 6-quark decay. In addition to the cuts described above, one non-isolated muon with p T > 15 GeV 
was required. For an integrated luminosity of 10 fb _1 , the expected signal would be about 7250 events, 
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Fig. 22: Schematics of the top decay to leptonic final states with J/xj). 

yielding a statistical uncertainty on the measurement of rat of approximately ±1 GeV. This technique 
is insensitive to the jet energy scale. The dominant uncertainties arise from effects of ISR and FSR and 
from the 6-quark fragmentation, which sum up to about 1.5 GeV. 



4.6 nit from t^l + J/V> + X decays 

An interesting proposal [ 111 ] by CMS, explored in detail during the workshop [112], is to take advantage 
of the large top production rates and exploit the correlation between the top mass and the invariant mass 
distribution of the system composed of a J/ip (from the decay of a b hadron) and of the lepton (£ = e,fi) 
from the associated W decay (see Fig. f22|). 

The advantage of using a J/ip compared to the other studies involving leptons as presented above 
is twofold: first, the large mass of the J/ip induces a stronger correlation with the top mass (as will be 
shown later). Second, the identification of the J/ip provides a much cleaner signal. In order to uniquely 
determine the top decay topology one can tag the charge of the b decaying to J/ip by requiring the 
other fe-jet to contain a muon as well. The overall branching ratio is 5.3 x 1CT 5 , taking into account the 
charge conjugate reaction and W — > ev decays. In spite of this strong suppression, we stress that these 
final states are experimentally very clean and can be exploited even at the highest LHC luminosities. 
Furthermore, one can also explore other ways to associate the J/ip with the corresponding isolated 
lepton - for example by measuring the jet charge of identified 6's. One should say that all these methods 
of top mass determination essentially rely on the Monte-Carlo description of its production and decay. 
Nonetheless the model, to a large extent, can be verified and tuned to the data. 



4.61 Analysis 

In the following we assume a tt production cross-section of 800 pb for mt = 175 GeV. Events are sim- 
ulated with the PYTHIA5.7 [52] or HERWIG 5.9 [|T]] event generators. Particle momenta are smeared 
according to parameterisations obtained from detailed simulation of the CMS detector performance. 
Four-lepton events are selected by requiring an isolated lepton with p T >15 GeV and |ry|<2.4, and three 
non-isolated, centrally produced muons of p T >4 GeV and |ry|<2.4, with the invariant mass of the two 
of them being consistent with the J / ip mass. These cuts significantly reduce the external (non-tt) back- 
ground, mainly Wbb production^ which can be further reduced by employing, in addition, two central 
jets from another W. The resulting kinematical acceptance of the selection criteria is 30%; this rather 
small value is largely due to soft muons from J/ip and b. In one year high luminosity running of LHC, 
corresponding to an integrated luminosity of 100 fb^ 1 , and assuming trigger plus reconstruction effi- 
ciency of 0.8, we expect about 10 5 x 800 x 5.3 • 10~ 5 x 0.3 x 0.8 = 1000 events. 

An example of the I J/ip mass distribution with the expected background is shown in Fig. ^3|. The 
background is internal (from the tt production) and is due to the wrong assignment of the J/ip to the 
corresponding isolated lepton. These tagging muons of wrong sign are predominantly originating from 

4 PYTHIA results indicate that with the above cuts this source of the background can be kept at a per cent level. 
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B°/B° oscillations, b— >c— >fj, transitions, W(— >c, r)— >/i decays, 7r/K decays in flight and amount to 
~ 30% of the signal combinations. The shape of the signal £J/ip events (those with the correct sign 
of the tagging muon) is consistent with a Gaussian distribution over the entire mass interval up to its 
kinematical limit of ~ 175 GeV. The background shape is approximated by a cubic polynomial. The 
parameters of this polynomial are determined with "data" made of the wrong combinations of £J/ip with 
an admixture of signal. In such a way the shape of the background is determined more precisely and in 
situ. Thus, when the signal distribution is fitted, only the background normalisation factor is left as a free 



parameter along with the three parameters of a Gaussian. The result of the fit is shown in Fig. 23. We 
point out that this procedure allows to absorb also the remaining external background (if any) into the 
background fit function. 
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Fig. 24: Correlation between M max and the top quark mass 
in isolated lepton plus J/ip (solid line) and isolated lepton 
plus /i-in-jet (dashed line) final states. 



As a measure of the top quark mass we use the mean value (position of the maximum of the 
distribution) of the Gaussian, MJfj^. In four years running at LHC with high luminosity the typical 
errors on this variable, including the uncertainty on the background, are about 0.5 GeV. It is composed 
of < 0.5 GeV statistical error and < 0.15 GeV systematics contribution due to the uncertainty on the 
measurement of the background shapej^] 

The measurement of the MYf/l can then be related to the generated top quark mass. An example 



of the correlation between the M max and m t is shown in Fig. [24J along with the parameters of a linear 
fit. For comparison, we also show the corresponding dependence in a more traditional isolated lepton 
plus /u-in-jet channel. Not surprisingly, the stronger correlation, and thus a better sensitivity to the top 
mass, is expected in the I J '/ip final states as compared to the isolated lepton plus /i-in-jet channel. This 
is because, in the former case, we pickup a heavy object (the J/ip) which carries a larger fraction of the 
6-jet momentum. The measurement error, statistical and systematic, scales as the inverse slope 

value of the fit, which is a factor of 2 in our case. Hence the statistical error on the top mass in this 
particular example is ~ 1 GeV. 

It is appropriate to comment on the ways to obtain a larger event sample. Encouraging results 
have been obtained in [113] to reconstruct the b— >J/ip— >e + e~ decays for low luminosity runs. The 
extension of these studies for a high luminosity environment is very desirable. Another possibility would 
be to relax the kinematical requirements. The choice of pt cut on soft muons is not dictated by the 



s The statistical power of the sample can be further improved by exploiting full spectrum, rather than its Gaussian part. 
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background considerations but by the trigger rates, and is set here to 4 GeV rather arbitrarily. For 
example, the di-muon trigger with ry-dependent thresholds which is available in CMS for low luminosity 



runs [114] allows to significantly increase the kinematical acceptance, practically to the limit determined 
by muon penetration up to the muon chambers. Therefore, the assessment of the trigger rates at high 
luminosity with lower px thresholds and in multi-lepton events clearly deserves a dedicated study. 

An even larger event sample can be obtained in three lepton final states, using instead the jet-charge 
technique to determine the tt decay topology instead of the tagging muon. The jet charge is defined as 
a pr-weighted charge of particles collected in a cone around the J/ip direction. Obviously, this kind 
of analysis requires detailed simulations with full pattern recognition which are under way. However, 
particle level simulations performed with PYTHIA and with realistic assumptions on track reconstruc- 
tion efficiency give event samples comparable to the muon-tag performance, with about 10 times less 
integrated luminosity. In any case, through the LHC lifetime, one can collect enough events so that the 
overall top mass measurement accuracy would not be hampered by the lack of statistics; it would rather 
be limited by the systematic uncertainties which are tightly linked with the Monte-Carlo tools in use, as 
will be argued in the following section. 

4.62 Systematics 

An essential aspect of the current analysis is to understand limitations which would arise from the Monte- 
Carlo description of the top production and decay. It is important to realize that the observable used in 
this study enjoys two properties: it is Lorentz invariant an it does not depend on the detailed structure of 
the jets, but only on the momentum spectrum of the 6-hadron and of the J/ip from its decay. 

As a result, were it not for distortions of the £J/ip mass distribution induced by acceptance effects 
and by the presence of an underlying background, the measurement would be entirely insensitive with 
respect to changes in the top production dynamics, and in the structure of the underlying event. As a 
result, typical systematics such as those induced by higher-order corrections to the production process, 
or by the ISR and by the structure of the minimum bias event, are strongly reduced relative to other 
measurements of m t . This expectation will be shown to be true in the following of this section. 

The main limitations to an accurate extraction of the top mass using this technique are expected 
to come from: i) the knowledge of the fragmentation function of the b hadrons contained in the 6-jet 
and, ii) the size of the non-perturbative corrections to the relation between the top quark mass and the 
iJ/ip mass distribution. The J/ip spectrum in the decay of the 6-hadrons will be measured with high 
accuracy in the next generation of S-factory experiments. It should be pointed out, however, that the 
composition of 6-hadrons measured at the T(45) and in the top decays will not be the same. In this 
second case, one expects a non-negligible contribution from baryons and from B s states. The size of the 
relevant corrections to the inclusive J/ip spectrum in top decays is not known, and, although expected to 
be small, it needs to be studied. Additional effects, such as QED corrections to the W leptonic decay, W 
polarisation and spin correlation effects can all be controlled and included in the theoretical simulations. 



The rest of this section presents the results of a detailed study [112] of the systematics, mostly 
based on PYTHIA. 

Detector resolution: Here we have considered only Gaussian smearing of particle momenta and the effect 
on the MJf 1 ^ measurement uncertainty is negligible. A possible nonlinearity of the detector response 
can be well controlled with the huge sample of J/ip, T and Z leptonic decays that will be available. 

Background: The uncertainty would be mainly due to an inaccurate measurement of the background 
shape and the systematics contribution of < 0.15 GeV quoted in previous section would scale down with 
increasing statistics. For example, already with ~ 10 4 events the induced uncertainty is < 0.1 GeV. 

PDF: Depending on the relative fraction of gluon/quarks versus x in various PDF's the top produc- 
tion kinematics might be different. No straightforward procedure is available for the moment to evalu- 
ate uncertainties due to a particular choice of PDF. We compared results obtained with the default set 
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CTEQ2L [115] and a more recent CTEQ4L [116] parameterisations of PDF's. The observed change in the 
Mflffij, value is well within 0.1 GeV. 

Top p T spectrum: As shown in Section p.3[ one does not expect significant uncertainties in the prediction 
of the top p T spectrum. However, to see an effect we have artificially altered the top p T spectrum by 
applying a cut at the generator level. We found that even requiring all top quarks to have p T > 100 GeV 
gives rise to only a la change (±0.7 GeV) in the fitted value of mt. 

Initial state radiation: The MffiK value is unchanged even switching off completely the ISR. 

Top and W widths: Kinematical cuts that are usually applied affect the observed Breit-Wigner shape 
(tails) of decaying particles. Conversely, poor knowledge of the widths may alter the generated IJ/ip 
mass spectrum depending on the cuts. In our case, only a small change in the MJf 1 ^ value is seen 
relative to the zero-width approximation. 

W polarisation: A significant shift is found for the isotropic decays of W when compared to the SM 
expectation of its ~ 70% longitudinal polarisation. In future runs of the Tevatron the W polarisation 
will be measured with a ~ 2% accuracy [J)8|], and at the LHC this would be further improved, so that it 
should not introduce additional uncertainties in simulations. 

ti spin correlations: A "cross-talk" between t and t decay products is possible due to experimental cuts. 
To examine this effect in detail the 2^6 matrix elements have been implemented in PYTHIA preserving 
the spin correlations [117]. No sizeable difference in the MJ 1 ^ value is seen compared to the default 
2^2 matrix elements. 

QED bremsstrahlung: Only a small effect is observed when it is switched off. Furthermore, QED radia- 
tion is well understood and can be properly simulated. 

Final State Radiation: A large shift of ~ 7 GeV is observed when the FSR is switched off. This is due 
to the absence of evolution for the b quark, whose fragmentation function will be unphysically hard. 
To evaluate the uncertainty we varied the partem virtuality scale m m i n , the invariant mass cut-off below 
which the showering is terminated. A ±50% variation of it around the default (tuned to data) value of 
1 GeV induces an uncertainty of *o" 15 GeV. 

b fragmentation, except FSR: As a default, in PYTHIA we have used the Peterson form for the 6-quark 
fragmentation function with Eb = 0.005. Variation of this value by ±10% [118] leads to an uncertainty 
of +025 GeV. (The ±10% uncertainty on £j is inferred from LEP/SLD precision of ~ 1% on the average 
scaled energy of i?-hadrons.) It should be pointed out that recent accurate measurements of the 6-quark 
fragmentation function [ |119| ] are not well fitted by the Peterson form. 

The last two items of this list deserve some additional comments. While the separation between 
the FSR and the non-perturbative fragmentation phases seems unnecessary, and liable to lead to an over- 
estimate of the uncertainty, it is important to remark that our knowledge of the non-perturbative hadroni- 
sation comes entirely from the production of 6-hadrons in Z° decays at LEP and SLC. It is important to 
ensure that the accuracy of both perturbative and non-perturbative effects is known, since the perturbative 
evolution of b quarks from Z° and top decays are not the same owing to the different scales involved. An 
agreement between data and Monte Carlo calculations for the 6-hadron fragmentation function at the Z° 
does not guarantee a correct estimate of the 6-hadron fragmentation function in top decays. 

To be specific, we shall consider here the effects induced by the higher-order matrix element 
corrections to the radiative top decays t^bWg [j54]]. These effects cannot be simulated by a change 
in the virtuality scale m m i n as explored above in the study based on PYTHIA, as they have a different 
physical origin. The extended phase-space available for gluon emission after inclusion of the matrix- 
element corrections leads to a softening of the 6-quark, and, as a result, of the £J/ip spectrum. For 
simplicity, we study here the invariant mass of the system B£. The resulting invariant mass distributions, 
for mt= 175 GeV, with (HERWIG 6.1) and without (HERWIG 6.0) matrix element corrections are shown 
in Fig. The averages of the two distributions, as a function of the top mass, are given on the right 
of the figure, and the difference of the averages are given in Table ^. Given the slopes of the correlation 
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Fig. 25: Left: invariant mass of the B-lepton system for mt = 175 GeV, according to HERWIG 6.0 (dotted) and 6.1 (with 
matrix element corrections, solid). Right: linear fits to the average invariant mass {tubi) as a function of mt. 

Table 7: Negative shift in the average invariant mass (mse) after inclusion of matrix element corrections for the top decay in 
HERWIG. Left: average over all values of msi- Right: average over the sample with mse > 50 GeV. 



m t 


<mj$) - (m B f) (all m Bi ) 


( m Be) ~ ( m jgi) i m Bi > 50 GeV) 


171 GeV 


(0.891 ± 0.038) GeV 


(0.479 ± 0.036) GeV 


173 GeV 


(0.844 ± 0.038) GeV 


(0.479 ± 0.034) GeV 


175 GeV 


(0.843 ± 0.039) GeV 


(0.510 ±0.035) GeV 


177 GeV 


(0.855 ± 0.039) GeV 


(0.466 ± 0.035) GeV 


179 GeV 


(0.792 ± 0.040) GeV 


(0.427 ± 0.036) GeV 



between (niBe) and mt, we see that the corrections due to inclusion of the exact matrix elements are 
between 1 GeV (for rnsi > 50 GeV) and 1.5 GeV (for the full sample). 

More details of the analysis will be found in [54]. It is also found there that the dependence of 
{mBi) on the hadronic center of mass energy, or on the partonic initial state producing the ti pair, is no 
larger than 100 MeV We take this as an indication that the effects of non-factorisable non-perturbative 
corrections (such as those induced by the neutralisation of the colour of the top quark decay products) 
are much smaller than the 1 GeV accuracy goal on the mass. 

A summary of these studies is given in Fig. |2^. One sees an impressive stability of the re- 
sults for reasonable choices of parameters. The expected systematic error in the M^Jj^ determination 

is ^ -o '1 GeV which translates into a systematic error on the top mass of 5m t < tol GeV 

In addition to the above studies, we also compared directly the results of HERWIG (v5.9) and 
PYTHIA. With HERWIG we have tried various tunings from LEP experiments as well as its default 
settings [51]. They all yield comparable results to each other and to PYTHIA results, and are within 
< 0.5 GeV This corresponds to a systematic uncertainty 5m t < 1 GeV 



4.7 Conclusions for the top mass measurement at the LHC 

The very large samples of top quark events which will be accumulated at the LHC lead to a precision 
measurement of the top quark mass. Different statistically independent channels have been investigated 
and from the studies so far a precision of better than 2 GeV in each case can be obtained. In particular 
for the lepton plus jets channel where the m t is measured directly reconstructing the invariant mass of 
the mjjb candidates, such a precision can be achieved within a year or running at low luminosity. For 
the channels involving two or more leptons, data from several years have to be combined to limit the 
statistical error in the measurement beyond the expected systematic errors. 
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Fig. 26: Observed m t shifts for the various systematic effects studied for the £ + J/ip channel. 



With the statistical error not being a problem, the emphasis of the work was devoted to estimate 
the systematic error involved in each method. For each sample, the contributing systematic errors are 
different, a fact which will allow important cross-checks to be made. The results indicate that a total 
error below 2 GeV should be feasible. In the case of the lepton plus jet channel the major contribution 
to the uncertainty is identified in the jet energy scale (in particular for the 6-jets) and in the knowledge 
of FSR. When a special sub-sample of high p T top events is used and the m t is reconstructed using a 
large calorimeter cluster the FSR sensitivity is reduced, but further work is required to validate it. For 
the channels using two or more leptons for the top decay, the major contribution in the systematic error 
comes from the Monte Carlo and from how well the kinematic observable used for the mass measurement 
is related to the mass of the top quark. 

In U/t/j final states the top mass can be determined with a systematic uncertainty of < 1 GeV. 
These final states are experimentally very clean and can be exploited even at highest LHC luminosities. 
The precision would be limited by the theoretical uncertainties which is basically reduced to the one 
associated with the t—>B meson transition. This method of top mass determination looks very promising, 
and a final definition of its ultimate reach will rely on a better understanding of theoretical issues, and on 
the possibility to minimise the model dependence using the LHC data themselves. 



5. SINGLE TOP PRODUCTION]] 

At the LHC, top quarks are mostly produced in pairs, via the strong process gg—>ti (and, to a lesser 
extent, qq—>ti). However, there are a significant number of top quarks that are produced singly, via the 
weak interaction. There are three separate single-top quark production processes of interest at the LHC, 
which may be characterised by the virtuality of the W boson (of four-momentum q) in the process: 

• t-channel: The dominant process involves a space-like W boson (q 2 < 0), as shown in Fig. |27|(a) 
[120]. The virtual W boson strikes a b quark in the proton sea, promoting it to a top quark. This 



6 Section coordinators: S. Willenbrock, D. O'Neil (ATLAS), J. Womersley (CMS). 
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b^ 

(a) (b) (c) 

Fig. 27: Feynman diagrams for single-top quark production in hadron collisions: (a) t-channel process; (b) s-channel process; 
(c) associated production (only one of the two diagrams for this process is shown). 

Table 8: Total cross sections (pb) for single-top quark production and top quark pair production at the LHC, for m t =175 



±2 GeV. The NLO t-channel cross section is from [125]. The NLO s-channel cross section is from [ 12C ]. The cross section 
for the Wt process is from [ 124 ]; it is leading order, with a subset of the NLO corrections included. The uncertainties are due 
to variation of the factorisation and renormalisation scales; uncertainty in the parton distribution functions; and uncertainty in 
the top quark mass (2 GeV). 



process: 


t-channel s-channel Wt tt 




245 ±27 10.2 ±0.7 51 ± 9 ~ 800 



process is also referred to as VF-gluon fusion, because the b quark ultimately arises from a gluon 
splitting to bb. 

s-channel: If one rotates the t-channel diagram such that the virtual W boson becomes time-like, 
as shown in Fig. |7|(b), one has another process that produces a single top quark [ |121 , 122 ]. The 
virtuality of the W boson is q 2 > (m t + mi,) 2 . 

Associated production: A single top quark may also be produced via the weak interaction in 



association with a real W boson (q 2 = Myy), as shown in Fig. 27(c) [123, 124]. One of the initial 



partons is a b quark in the proton sea, as in the t-channel process. 

The total cross sections for these three single-top quark production processes are listed in Table |8L 
along with the cross section for the strong production of top quark pairs. The t-channel process has 
the largest cross section; it is nearly one third as large as the cross section for top quark pairs. The s- 
channel process has the smallest cross section, more than an order of magnitude less than the t-channel 
process. The Wt process has a cross section intermediate between these two. We will argue that all three 
processes are observable at the LHC. The t-channel and s-channel processes will first be observed at the 



Fermilab Tevatron Q127D ; the Wt process will first be seen at the LHC. 

There are several reasons for studying the production of single top quarks at the LHC: 

• The cross sections for single-top quark processes are proportional to \ V t iy\ 2 . These processes pro- 
vide the only known way to directly measure Vtb at hadron colliders. 

• Single-top quark events are backgrounds to other signals. For example, single-top quark events 



are backgrounds to some signals for the Higgs boson [128]. 



Single top quarks are produced with nearly 100% polarisation, due to the weak interaction [ |123 



12S , |130| , |1 3 lp - This polarisation serves as a test of the V — A structure of the top quark charged- 
current weak interaction. 

• New physics may be discernible in single-top quark events. New physics can influence single-top 
quark production by inducing non-SM weak interactions [ 129 , 132| , 133, 134| , |135[ ], via loop effects 
J136| , |137| , |138| , 139, 140], or by providing new sources of single-top quark events [133, 137[ , 141 , 

In the next three subsections we separately consider the three single-top quark production pro- 
cesses. The subsection after these discusses the polarisation of single top quarks. In the concluding 
section, we discuss the accuracy with which V t b can be measured in single-top quark events at the LHC. 
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5.1 t-channel single-top production 

5.11 Theory 

The largest source of single top quarks at the LHC is via the t-channel process, shown in Fig. |27|(a) 
[ [i20| , |23| , p5j |l29i (T43J, [144[ p5| |. A space-like (g 2 < 0) W boson strikes a 6 quark in the proton sea, 
promoting it to a top quark. As shown in Table H the cross section for this process is about one third that 
of the strong production of top quark pairs. Thus there will be an enormous number of single top quarks 
produced via the t-channel process at the LHC. 

It is perhaps surprising that the cross section for the weak production of a single top quark, of 
order ctyy, is comparable to that of the strong production of top quark pairs, of order a 2 . There are 
several enhancements to the t-channel production of a single top quark that are responsible for this: 

• The differential cross section for the t-channel process is proportional to da/dq 2 ~ l/(q 2 —M^ v ) 2 , 
due to the VF-boson propagator. The total cross section is therefore dominated by the region 
\q 2 \ < Myy, and is proportional to 1/My^. In contrast, the total cross section for the strong 
production of top quark pairs is proportional to 1/s, where s > Am 2 is the parton center-of-mass 
energy. 

• Since only a single top quark is produced, the typical value of the parton momentum fraction x is 
half that of top quark pair production. Since parton distribution functions scale roughly like l/x 
at small values of x, and there are two parton distribution functions, this leads to an enhancement 
factor of roughly four. 

The fact that the total cross section is dominated by the region | q 2 \ < also has the implication that 
the final-state light quark tends to be emitted at small angles, i.e., high rapidities. This characteristic 
feature of the signal proves to be useful when isolating it from backgrounds. 

The b distribution function in the proton sea arises from the splitting of virtual gluons into nearly- 
collinear bb pairs. Thus it is implicit that there is a b in the final state, which accompanies the top quark 
and the light quark. The final-state b tends to reside at small px, so it is usually unobservable. 

The total cross section for the t-channel production of single top quarks has been calculated at 
NLO [ 125| , 143]; the result is given in Table H. A subset of the NLO corrections is shown in Fig. 28 |a). 
This correction arises from an initial gluon which splits into a bb pair. If the bb pair is nearly collinear, 
then this process contributes to the generation of the b distribution function, which is already present at 
leading order; hence, one does not include this kinematic region as a contribution to the NLO correction. 
This is indicated schematically in Fig. |28|(b). Only the contribution where the bb pair is non-collinear is a 
proper NLO correction to the total cross section.^] The other corrections to this process, due to final-state 
and virtual gluons, as well as corrections associated with the light quark, are also included in the cross 
section given in Table |g. 

The central value for the cross section is obtained by setting the factorisation scale^] of the b distri- 
bution function equal to /x 2 = — q 2 + m 2 . The uncertainty in the NLO cross section due to the variation 
of the factorisation scale between one half and twice its central value is 4%. Due to the similarity with 
deep-inelastic scattering, the factorisation scale of the light quark is /i 2 = —q 2 , and is not varied [ 125]. 

Since the b tends to reside at low pr, the dominant final state is Wbj, where the Wb are the decay 
products of the top quark, and the jet is at high rapidity. However, the b is at pr > 20 GeV in roughly 40% 
of the events, in which case the final state is Wbbj. From a theoretical perspective, the optimal strategy 
is to isolate both final states and thereby measure the total cross section, which has an uncertainty of 
only 4% from varying the factorisation scale, as mentioned above. However, the Wbbj final state has a 
large background from tt, and it has not yet been established by ATLAS or CMS that this signal can be 
isolated, although the analysis of [145] gives cause for optimism. Thus we focus on the Wbj final state, 

7 The formalism for s epar ating the neaiiy-collinear and non-collinear regions, and for generating the b distribution function, 
was developed in Refs. [146, 58j |, 

8 The factorisation and renormalisation scales are set equal. 
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Fig. 28: (a) Initial-gluon correction to single-top quark production via the f-channel process (the diagram with the W and 
gluon lines crossed is not shown); (b) the kinematic region in which the gluon splits to a nearly-collinear bb pair (the double 
line through the b propagator indicates that it is nearly on shell) is subtracted from the correction, as it is already included at 
leading order. 

demanding that the b have px < PTcut- For pxcut = 20 GeV,f] the cross section for this semi-inclusive 



process is 164 pb, with an uncertainty of 10% from varying the factorisation scale [144], about twice 
the uncertainty of the total cross section. Work is in progress to calculate the differential cross section 
da/dprpi at NLO with the goal of reducing this uncertainty [ ]147| ]. It would also be desirable to calculate 
the total cross section at next-to-next-to-leading order (NNLO). 

Additional theoretical uncertainties stem from the top quark mass and the partem distribution func- 
tions. An uncertainty in the top quark mass of 2 GeV yields an uncertainty of only 2% in the cross section, 
which is negligible. This is due to the fact that the cross section scales like 1/M^ rather than 1/s. The 
uncertainty in the cross section due to the partem distribution functions is estimated in [148] to be 10%. 
That analysis suggests that the uncertainty can be reduced below this value. Combining all uncertainties 
in quadrature, we conclude that the total theoretical uncertainty is presently 15% in the Wbj cross sec- 
tion (11% in the total cross section). The discussion above suggests that this can be significantly reduced 
with further effort. 



5.12 Phenomenology 

Studies of the i-channel process have been carried out by both ATLAS and CMS. We will first describe 
the CMS study, and then that of ATLAS. 

In order to reject the large it background in this channel, it is necessary to impose a cut on jet 
multiplicity. Accurate modelling of jet response and resolution is therefore desirable, and so CMS [ 149 ] 
used a full GEANT calorimeter simulation of the detector. The GEANT simulation also allows a more 
realistic modelling of the missing-p^ response of the detector, which is important in understanding the 
mass resolution which can be obtained on the reconstructed t quark. The detailed calorimeter simulation 
was combined with a parameterised 6-tagging efficiency. 

Signal events were generated using PYTHIA 5.72 [p2|], with m t = 175 GeV and the CTEQ2L 
parton distribution functions. Events were preselected at the generator level to have one and only one 
charged lepton (with px > 25 GeV and \t]\ < 2.5) and one or two jets (generator-level jets were found 
using the LUCELL clustering algorithm, which is part of PYTHIA). Generated events were then passed 
through the parameterised 6-tagging and the GEANT detector simulation. The CMS 6-tagging perfor- 
mance is taken from a study which used a detailed detector simulation combined with existing CDF data 
on impact-parameter resolutions. The tagging efficiency for pt > 50 GeV is typically 50% for 6-jets, 
10% for c-jets, and 1-2% for light quarks and gluons. These efficiencies fall quite rapidly for lower 
transverse momenta, and it was assumed no tagging could be performed for px < 20 GeV or \rj\ > 2.4. 
The generated luminosity corresponded to about 100 pb -1 - only 30 hours of running at 10 33 cm~ 2 s _1 . 

The ti and WZ backgrounds were also generated using PYTHIA 5.72. The same pre-selections 
were applied at the generator level. The W+ jets backgrounds were generated using the VECBOS 

'The CMS analysis presented below uses prcut = 20 GeV; the ATLAS analysis uses prcut = 15 GeV. 
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Top moss (GeV) 

Fig. 29: Reconstructed top mass for signal plus backgrounds (open histogram) and backgrounds only (shaded). The back- 
grounds considered are ii,W + 2 jets and W + 3 jets. The vertical scale is events per 6 GeV mass bin per pb~ 1 of luminosity. 



generator [150], combined with HERWIG 5.6 [51] to fragment the outgoing partons.[^] W + 2 jets and 
W + 3 jets processes were generated separately. Again, events were preselected to have a charged lepton 
with px > 25 GeV and \r)\ < 2.5, and to have a (parton-level) pp > 15 GeV for the final-state jets. 

Events were then selected which passed the following requirements: 

• One and only one isolated lepton (£ = e or fj) with px > 20 GeV and \rj\ < 2.5. This allows the 
events to pass a reasonable lepton trigger. 

• Missing pr > 20 GeV, and transverse mass (of the lepton and missing pr) 50 < mx < 100 GeV. 
These two requirements select W^iv candidates. 

• Exactly two jets withpr > 20 GeV and \r\\ < 4. Requiring at least two jets reduces the W+ jets 
background, while requiring no more than two jets rejects the tt background which naively would 
produce four jets in the final state. 

• One jet withp T > 20 GeV and \rj\ < 2.5, the other jet with p T > 50 GeV and 2.5 < \rj\ < 4.0. 
The requirement that the second jet be at forward rapidities tends to select the desired t-channel 
process. 

• Leading jet pp < 100 GeV. This helps to reduce the tt background. 

• Exactly one 6-tagged jet (given the 6-tagging acceptance, this is always the central jet). This 
requirement again reduces tt, and of course rejects W+ jets processes with light-quark or gluon 
jets. 

• Invariant mass of the two jets in the 80 — 100 GeV range. This rejects WZ events with Z— >bb. 

The single-top signal is then searched for in the invariant mass of the W and the fe-tagged jet (which 
should peak at the top quark mass). The mass was reconstructed assuming the solution for the W 
kinematics which yields the lower \p v z \. (It is possible to use other choices, for example the solution 
which gives the Wb mass closest to m t . This would result in an apparently better top mass resolution 
but would also severely bias the background shape; the statistical significance of the signal would not be 
improved.) 



Figure |29| shows the reconstructed mass distribution for signal and background combined. The 
signal is apparent as an excess over the background (the shaded histogram) around 160 GeV. (Since jet 



10 The version of VECBOS used here, and its interface to HERWIG, were developed for use in CDF [151], and were adapted 
for CMS by R. Vidal. 
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energy scale corrections have not been applied to the simulated events, the top mass reconstructs to less 
than its true value.) The signal-to-background ratio in a window of 160 ± 20 GeV is 3.5 with a clear peak 
visible in the Wb invariant-mass distribution. The number of signal events is 66 in 100 pb _1 , giving a 
signal efficiency of 1.2% (after the W^iv branching ratio). We then find that 10 fb _1 would yield 6600 
signal events (S) and 1900 background (B), sufficient for a statistical accuracy on the number of signal 
events of y/S + B/S = 1.4%. 

The largest background comes from Wcj with the charm jet mistagged as a fe-jet. It would be 
worthwhile to develop a 6-tagging algorithm having greater rejection against such mistags, even at the 
cost of some signal efficiency. The Wbb background was found to be a small contribution to the W + 2 
jets background at the parton level for the selection cuts employed here, and was therefore not explicitly 
included in the analysis. 

The use of the forward jet tag substantially improves the signal-to-background ratio, and allows 
a clear reconstructed top-mass peak to be seen. However, it does not significantly improve V 'S + B / ' S 
[144]. One could therefore imagine omitting the forward jet requirement if the systematic uncertainty 
could thereby be reduced. 



Compared with earlier studies (for example [144]), this analysis uses more realistic jet and mis- 
sing-pr resolutions, and includes initial- and final-state gluon radiation. As a result, the top-mass res- 
olution is worsened; but the resolution found here compares well with the result of a full simulation of 
single-top production in CDF. 

A study of the cross-section measurement for the t-channel process was also carried out by ATLAS 
[152]. Signal events were generated using the ONETOP parton-level Monte Carlo [153] with fragmen- 
tation, radiation, and underlying event simulated by PYTHIA 5.72. Backgrounds containing top quarks 
(ti and other single-top production) were also generated using ONETOP, while W+ jets and Wbb back- 
grounds were generated by HERWIG 5.6.[]] These events were processed by the ATLAS parameterised 
detector simulation assuming a 60% 6-tagging efficiency for 6-jets, 10% for c-jets, and 1% for light 
quarks and gluons. The events were then analysed with a view towards separating t-channel single top 
from background and measuring its cross section. 

Event selection criteria were divided into two types: pre-selection and selection cuts. The pre- 
selection criteria were as follows: 

• at least one isolated lepton with pt > 20 GeV; 

• at least one 6-tagged jet with pr > 50 GeV; 

• at least one other jet with pr > 30 GeV. 
These were followed by the selection cuts: 

• two and only two jets in the event (a jet has pt > 15 GeV); 

• one jet is a central 6-tagged jet; 

• the other jet is a forward (\rj\ > 2.5) untagged jet with pr > 50 GeV. 

The application of these cuts, and also the requirement of a reconstructed top mass between 150 
and 200 GeV, yields the number of events shown in Table ^. The final signal efficiency is 3% and the 
signal-to-background ratio is 2.4. This implies a statistical precision on the cross-section measurement of 
VS 1 + B/S = 0.9% with 10 fb _1 of data. Introducing other event selection variables (see [fo[ IBij [l55| ]) 
it is possible to improve the signal-to-background ratio to nearly 5, but this does not improve the cross- 
section measurement due to the small remaining signal efficiency. 

Both the CMS and ATLAS studies indicate that it will be possible to observe t-channel single-top 
production with a good signal-to-background ratio and a statistical uncertainty in the cross section of 
less than 2% with 10 fb -1 . Thus the uncertainty in the extracted value of V t b will almost certainly be 
dominated by systematic uncertainties, as discussed in the conclusions. 

"The Wbb background was generated using the matrix element from m% interfaced to HERWIG 5.6. 
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Table 9: Cumulative effect of cuts on f-channel signal and backgrounds. The first four rows of this table refer to cumulative 
efficiencies of various cuts. The last two rows refer to the number of events for 10 fb _1 . Only events in which W — > ev or fxv 
are considered in this table. Uncertainties quoted in this table are due entirely to Monte Carlo statistics. 



cut 


f-channel 

eff(%) 


tt 
eff(%) 


Wbb 
eff(%) 


IV + jets 
eff(%) 


pre-selection 


18.5 


44.4 


2.53 


0.66 


njets=2 


12.1 


0.996 


1.55 


0.291 


fwd jet 

M > 2.5 
p T >50 GeV 


4.15 


0.035 


0.064 


0.043 


M M 
150-200 GeV 


3.00 


0.017 


0.023 


0.016 


events/10 fbr 1 
(before cuts) 


5.43 x 10 5 


2.40 x 10 6 


6.67 x 10 5 


4.00 x 10 7 


events/10 fbr 1 
(after cuts) 


16515±49 


455±74 


155±17 


6339±265 



5.2 s-channel single-top production 

5.21 Theory 



The s-channel production of single top quarks is shown in Fig. |27|(b) Jl21| , |122[ , |123| , |126|, |144| , |145J. The 
cross section is much less than that of the i-channel process because it scales like 1/s rather than 
However, the s-channel process has the advantage that the quark and antiquark distribution functions are 
relatively well known, so the uncertainty from the partem distribution functions is small. Furthermore, 
the parton luminosity can be constrained by measuring the Drell-Yan process qq^W*^£i>, which has 
the identical initial state | ff22| , p36|].PI 

The total cross section for the s-channel process has been calculated at NLO [126]; the result is 
given in Table || The factorisation and renormalisation scales are set equal to fi 2 = q 2 ; varying each, 
independently, between one-half and twice its central value yields uncertainties in the cross section of 2% 
from each source. The uncertainty in the cross section from the parton distribution functions is estimated 
to be 4%. The largest single source of uncertainty is the top quark mass; an uncertainty of 2 GeV yields 
an uncertainty in the cross section of 5%. The relatively large sensitivity of the cross section to the top 
quark mass is a manifestation of the 1/s scaling. Combining all theoretical uncertainties in quadrature 
yields a total uncertainty in the cross section of 7%. This is much less than the present theoretical 
uncertainty in the i-channel cross section. 

The Yukawa correction to this process, of order cn^mf /M^, is less than one percent [ |1 26| ] . 
However, this correction could be significant in a two-Higgs-doublet model for low values of tan (5, 
in which the Yukawa coupling is enhanced [|138p. 



5.22 Phenomenology 

In order to evaluate the potential to separate the s-channel signal from its backgrounds, Monte Carlo 
events have been processed by a fast (parameterised) simulation of an LHC detector. At parton level the 
signal and the ti background were generated by the ONETOP Monte Carlo [153]. Radiation, showering, 
and the underlying event were added by PYTHIA 5.72 [52]. The W+ jets and Wbb backgrounds were 
generated using HERWIG 5.6 Table || presents the cross sections assumed for the processes 

12 The parton luminosity can only be constrained, not directly measured, with this process. Since the neutrino longitudinal 
momentum is unknown, the q of the virtual W cannot be reconstructed. 

B The Wbb background was generated using the matrix element from [p9|] interfaced to HERWIG 5.6. 



44 



Table 10: Cumulative effect of cuts on s-channel signal and backgrounds. The first five rows of this table refer to cumulative 
efficiencies of various cuts. The last two rows refer to the number of events for 30 fb~ . Only events in which eu or \iv 
are considered in this table. Uncertainties quoted in this table are due entirely to Monte Carlo statistics. 



cut 


s-channel 

eff(%) 


i-channel 

eff(%) 


Wt 
eff(%) 


tt 
eff(%) 


Wbb 
eff(%) 


W+jets 
eff (%) 


pre-selection 


27.0 


18.5 


25.5 


44.4 


2.53 


0.667 


njets=2 


18.4 


12.1 


4.03 


0.996 


1.55 


0.291 


nbjet=2 

p T > 75 GeV 


2.10 


0.035 


0.018 


0.023 


0.034 


0.0005 


E ]ets PT 

>175 GeV 


1.92 


0.031 


0.016 


0.021 


0.028 


0.0005 


150-200 GeV 


1.36 


0.023 


0.006 


0.012 


0.0097 


0.00014 


events/30 fbr 1 
(before cuts) 


6.66 x 10 4 


1.63 x 10 6 


4.5 x 10 5 


6.9 x 10 6 


2.0 x 10 6 


1.2 x 10 8 


events/30 fbr 1 
(after cuts) 


908 ± 35 


375 ± 13 


27 ± 15 


853 ± 175 


194 ± 34 


169 ± 76 



containing top quarks. The cross section for the W + jets background is normalised to that predicted 
by the VECBOS Monte Carlo [|l5Q|] and is taken to be 18000 pb.0 The Wbb cross section is taken 
from Jig to be 300 pb. 

From a phenomenological standpoint the most important distinction between the s-channel and 
t-channel sources of single top is the presence of a second high-py b-jet in the s-channel process. As 
mentioned previously, in i-channel events the second 6-jet tends to be at low pr and is often not seen. 
Therefore, requiring two 6-jets above 75 GeV pt will eliminate most of the t-channel background. Re- 
quiring two high-pr b-jets m the event also suppresses the W+ jets background relative to the signal. 

In addition to suppressing the t-channel background it is also necessary, as in other single-top 
signals, to design cuts to reduce the W+ jets and tt backgrounds. In order to reduce contamination by 
W+ jets events, the reconstructed top mass in each event must fall within a window about the known top 
mass (150-200 GeV), and the events must have a total transverse jet momentumP] above 175 GeV. Only 
events containing exactly two jets (both tagged as b's) are kept in order to reduce the it background. 



Table |10J presents the cumulative effect of all cuts on the s-channel signal and on the backgrounds. 
Events from t-channel single-top production are included in this table as a background to the s-channel 
process. From this table the predicted signal-to-background ratio for the s-channel signal is calculated 
to be 0.56. The results also imply a signal statistical significance (S/y/B) of 23 with an integrated 
luminosity of 30 fb _1 . The statistical precision on the cross section, calculated from y/S + B/S, is 
5.5% with 30 ft.- 1 . 

This study indicates that, despite the large anticipated background rate, it should be possible to 
perform a good statistical measurement of the s-channel single-top cross section. The accuracy with 
which Vtb can be measured is discussed in the conclusions. 



5.3 Associated production 

5.31 Theory 



Single top quarks may also be produced in association with a W boson, as shown in Fig. 27 [c) [123, 124 



145]. Like the t-channel process, one of the initial partons is a b quark. However, unlike the t-channel 



14 This cross section is defined for events containing at least two jets, each withpr > 15 GeV and |r;| < 5. 
15 Scalar sum of the transverse momentum of all jets in the event. 
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process, this process scales like 1/s. This, combined with the higher values of x needed to produce both 
a top quark and a W boson, leads to a cross section for associated production which is about a factor of 
five less than that of the the t-channel process, despite the fact that it is of order a s aw rather than a^y. 

The total cross section for associated production has been calculated at leading order, with a subset 



of the NLO corrections included J124] , |145| ] ; the result is given in Table |[ This subset is analogous to the 
initial-gluon correction to the t-channel process, discussed previously. The other corrections have not 
yet been evaluated.^ The initial-gluon correction contains an interesting feature which has no analogue 
in the t-channel process. One of the contributing diagrams to the initial-gluon correction (gg^Wtb) 
corresponds to gg—>tt, followed by i— >Wb. This should not be considered as a correction to associated 
production, but rather as a background (it is in fact the dominant background, as discussed below). Thus, 
when evaluating the initial-gluon correction, it is necessary to subtract the contribution in which the t is 



on shell. This is done properly in Q124Q . 

The cross section is evaluated with the common factorisation and renormalisation scales set equal 
to /i 2 = s. The uncertainty in the cross section due to varying these scales between one half and twice 
their central value is 15%. This uncertainty would presumably be reduced with a full NLO calculation. 
The uncertainty in the cross section from the parton distribution functions is estimated to be 10% | ]148| ],[ rT 1 
although this could be improved with further study. The uncertainty in the cross section due to an 
uncertainty in the top quark mass of 2 GeV is 4%, relatively large due to the 1/s scaling of the cross 
section. Combining all theoretical uncertainties in quadrature yields a total uncertainty at present of 18%, 
the largest of the three single-top processes. 

5.32 Phenomenology 

The strategy for measuring the cross section for associated production (Wt mode) is similar to that for 
the t-channel process, as they share the same backgrounds. However, the nature of associated production 
makes it relatively easy to separate from W+ jets and difficult to separate from tt events. This difficulty 
in removing the tt background does not preclude obtaining a precise cross-section measurement in this 
channel, assuming the rate for tt can be well measured at the LHC. 

Two studies designed to separate signal from background have been performed using two different 
final states. The first is a study by ATLAS [00] which attempts to isolate Wt signal events in which one 



W decays to jets and the other decays to leptons. The second study, which is presented in [ |124| ], attempts 
to isolate signal events in which both W's decay leptonically. 

The first study presented here was done by ATLAS using the same event sample described in 



Section p.l[ Since the presence of a single isolated high-p^ lepton is one of the preconditions of this 
study, the second W must decay to two jets to be accepted by the event pre-selection. Therefore requiring 
a two-jet invariant mass within a window around the W mass will serve to eliminate most events that 



do not contain a second W. The two-jet invariant-mass distribution is shown in Fig. 30 and clearly 
demonstrates the presence of a sharp peak in the associated-production signal and the tt background. 
This effectively leaves tt as the only background to Wt events. 

In addition to these special distinguishing features of the Wt signal, there are several simple 
kinematic requirements which can be employed to reduce the tt background. By choosing events with 
exactly three jets and with exactly one of them tagged as a 6-jet, some rejection of the tt background is 
possible. Some further rejection is obtained by limiting the selection to events with invariant mass less 
than 300 GeV, where the invariant mass of an event is defined as the invariant mass obtained by adding 
the four-vectors of all reconstructed jets and charged leptons (e and /x). However, even with these cuts 
the tt background is significantly larger than the Wt signal. 



Table 1 1 presents the cumulative effect of all cuts on the Wt signal and on the tt and W+ jets 



The analogous calculation for Wc production has been performed in [157] 



7 This is the uncertainty in the gluon-gluon luminosity at yfr = (m t + Mw)/\/~S » 0.02, where \f§ = 14 TeV. 
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Fig. 30: The normalised two-jet invariant-mass distribution. For each event the two-jet combination with mass closest to the 
W mass is plotted. This clearly shows a peak in the distribution for Wt and ti which is not present for the other backgrounds. 

backgrounds. The Wbb and f-channel single-top backgrounds are virtually eliminated by the cuts and 
so are not included in the table. From this table the predicted signal-to-background ratio for the Wt 
signal is calculated to be 0.24. After three years of running at low luminosity (30 fb -1 ), this implies a 
signal statistical significance (S/ ^/B) of 25 and a statistical error on the Wt cross section (y/S + B / S) 
of 4.4%. 



The second study J124] ] was done at parton level and involved the separation of signal from back- 
ground in the mode in which both W's decay to leptons. This signal contains two high px leptons and 
only one jet (the 6-jet produced from the top decay). In this decay channel it was found that, after apply- 
ing detector acceptance cuts, requiring precisely one 6-tagged jet with px > 15 GeV is enough to yield 
a signal-to-background ratio of nearly unity. Also, the signal efficiency is significantly higher than in the 
ATLAS analysis, allowing more total signal events to pass the cuts despite the lower branching ratio for 
this decay mode. The statistical precision on the cross section measured in this analysis is 1.3% with 
an integrated luminosity of 30 fb _1 . The accuracy with which Vtb can be extracted is discussed in the 
conclusions. 



5.4 Polarisation in single-top production 

5.41 Theory 



Because single top quarks are produced through the weak interaction, they are highly polarised [123 



129|, N_30J, 131, 144]. In the ultra-relativistic limit, the top quarks are produced in helicity eigenstates 



with helicity —1/2 (the top antiquarks have helicity +1/2), because the V — A structure of the weak 
interaction selects quarks of a definite chirality. However, if the top quarks are not ultra-relativistic, 
chirality is not the same as helicity. Nevertheless, it was shown in [ |130| ] that there is a basis in which 
the top quark is 100% polarised, regardless of its energy. The top quark spin points along the direction 
of the d-type or d-type quark in the event, in the top quark rest frame (the t spin points opposite this 
direction). In i-channel production, this is the direction of the final-state light quark (ub—>dt) or the 
beam direction (db—>ut). In s-channel production, this is the beam direction (ud—>tb). In associated 
production (gb— >Wt), this is the direction of the d quark (or charged lepton) from the W decay. 

We focus our attention on the i-channel single-top process for the remainder of this section. The 
top quark polarisation in the t-channel process has been calculated at NLO [131 1; the results below are 
taken from this study. In the case of t production, 80% of the events have the d-type quark in the final 
state. This suggests using the direction of the light-quark jet, as observed in the top quark rest frame, to 



47 



Table 11: Cumulative effect of cuts on Wt signal and backgrounds. Pre-selection cuts are defined in the same way as for the 
ATLAS t-channel analysis described earlier in this report. The first five rows of this table refer to cumulative efficiencies of 
various cuts. The last two rows refer to the number of events for 30 fb _1 . Only events in which e,u or fiv are considered 
in this table. Uncertainties quoted in this table are due entirely to Monte Carlo statistics. 



cut 


Wt 
eff(%) 


tt 
eff(%) 


W+ jets 
eff(%) 


pre-selection 


25.5 


44.4 


0.66 


njets=3 
p T > 50 GeV 


3.41 


4.4 


0.030 


nbjet=l 

p T > 50 GeV 


3.32 


3.24 


0.028 


Invariant Mass 
< 300 GeV 


0.55 


0.36 


0.00051 


65 < M j:j < 95 


0.49 


0.14 


0.000085 




events/30 fb _1 
(before cuts) 


5.3 xlO 5 


7.2xl0 6 


1.2xl0 8 


events/30 ftr 1 
(after cuts) 


2608 ± 166 


10616 ± 625 


102 ± 59 



measure the spin. This has been dubbed the "spectator basis" [|130p . The polarisation of the top quark 
in this basis (defined as P = (JVj — N^)/(N^ + N\)) is 0.89. However, the polarisation is increased to 
nearly 100% when the cuts used in the i-channel analysis are imposed. This is because the polarisation is 
diluted by events in which the b is produced at high pt', but such events are eliminated by the requirement 
of only two jets. 

In the case of t production, 69% of the events have the <i-type quark in the initial state. This 
suggests using the beam direction to measure the t spin. However, it turns out that the spectator basis 
again yields the largest polarisation, P = —0.87. This polarisation is increased to P = —0.96 when cuts 
are applied.^ 

Since the top quark decays via the weak interaction, its spin is analysed by the angular distribution 
of its decay products. The most sensitive spin analyser in top decay is the charged lepton, which has a 
(leading order) angular distribution with respect to the top quark spin of 

i t^r i 

p-j ^- = -(l + cos^) (12) 

1 a cos Of 2 



in the top quark rest frame [158]. Hence the charged lepton tends to point along the direction of the 
spectator jet. 



5.42 Phenomenology 

The goal of this analysis is to estimate the sensitivity of ATLAS and CMS to the measurement of the 
polarisation of the top quarks produced by the f-channel single-top process. The f-channel process was 
chosen due to the large statistics available in this channel and the relative ease with which it is separated 
from its backgrounds. The t-channel events produced by the ONETOP generator and passed through 
PYTHIA and a parameterised detector simulation are analysed to attempt to recover the predicted SM 
top polarisation in the presence of background and detector effects. Details of the study are presented 
in 

18 With cuts applied, the polarisation in the so-called "77-beamline basis" is slightly higher, P = —0.97. 
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The experimental measurement of the polarisation of the top quark is essentially a measurement 
of the angular distribution of its decay products in the top quark rest frame. As explained above, the 
most sensitive angle is between the charged lepton from top decay and the direction of the spectator jet, 
in the top quark rest frame. In the absence of background or detector effects the angular distribution of 
the charged lepton is given by 

/(cos0*) = ~(l + i>cos0<) (13) 

where P is the polarisation of the sample and can range from —1 to 1. 

Experimentally, in order to measure the angular distribution of the charged lepton in the top quark 
rest frame, it is necessary to first reconstruct the four-momentum of the top quark. However, the recon- 
struction of the top four-momentum suffers from an ambiguity due to the unknown longitudinal momen- 
tum of the neutrino produced in the top decay. Using the W and top masses as constraints,^] one can 
reconstruct the top four-momentum, but the quality of the reconstruction is degraded by this ambiguity. 
Once the top four-momentum has been reconstructed, one can determine the direction of the spectator 
jet and the charged lepton in the top quark rest frame. The angle between these two directions is 6g. 

In order to extract the value of the top polarisation from the angular distribution, reference event 
samples were created with 100% alignment with the polarisation axis (spin up, P = +1) and with 
100% anti-alignment with the polarisation axis (spin down, P = —1). These reference distributions 
were compared to a statistically-independent data set with the predicted SM top quark polarisation. This 
comparison was done by minimising 



(cos 8), °th i ^°d i 



where the subscript d represents quantities calculated for the data distribution and the subscript th refers 
to the generated reference distribution. The theoretical value /th(cos Qi) is calculated via 

/thCcos e ) = i((l - P)/d(cos 6 t ) + (1 + P)/c/(cos 6 e )) (15) 

where fjj and fjj refer to the value of the generated theoretical distribution for the 100% spin-down and 
the 100% spin-up tops, respectively, and P is the polarisation of the top sample. The procedure returns 
an estimate of the top polarisation and an error on that estimate. In this way the sensitivity to changes in 
top polarisation can be quantified. 

Moving from the parton-level simulation to a simulation which includes both hadronisation and 
detector effects is certain to complicate the measurement of the polarisation of the top quark. In ad- 
dition, the signal could be biased by an event selection designed to eliminate background and will be 
contaminated by residual background events. 



The first histogram in Fig. [31] shows the angular distribution for signal only, at parton-level. The 
second histogram in Fig. |T] shows the angular distribution of the charged lepton after detector effects 
have been simulated. In addition to effects associated with detector energy smearing, jet and cluster 
definitions, etc. , this distribution includes the effects of ambiguities in reconstructing the top quark due to 
the absence of information about the neutrino longitudinal momentum. It does not, however, contain the 
effects of any event selection in order to separate signal from background. This histogram demonstrates 
that the effect of hadronisation and detector resolution changes the shape of the angular distribution but 
still produces a highly asymmetric distribution. 

In addition to the effects introduced by the detector resolution, the effect of applying the event- 
selection criteria can be evaluated by applying them one at a time and observing the change in shape of 
this distribution. For the purposes of the p olarisation analysis the event-selection criteria are: 

9 The W mass can be used to calculate the neutrino longitudinal momentum to within a two-fold ambiguity. Of these two 
solutions the one which produces the best top mass is chosen. 
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Fig. 31: Angular distribution of charged lepton in top rest frame for various data samples. The histograms progress from left- 
to-right, top-to-bottom. The first histogram shows the parton-level distribution. The second histogram is after the simulation of 
detector and reconstruction effects. The final four histograms illustrate the influence of event selection criteria on the angular 
distribution. The effects of the cuts are cumulative and are the result of adding pre-selection cuts, a jet-multiplicity requirement, 
a forward jet tag, and a top mass window, respectively. 



• Pre-selection (trigger) cuts as in ATLAS i-channel analysis described previously; 

• number of jets = 2; 

• forward jet (\r]\ > 2.5) withp r > 50 GeV; 

• reconstructed top mass in the range 150-200 GeV. 

This set of criteria leads to a signal efficiency of 3.0%, corresponding to more than 16000 events in 
10 fb _1 of integrated luminosity. Fig. 31 demonstrates the effect of applying these cuts in a cumulative 



manner. Again the asymmetry of the f-channel angular distribution is preserved, though more degrada- 
tion is clearly evident, in particular near cos 0% = 1. The degradation is worse at these values of cos 0£ 
because the leptons from these events are emitted in the direction opposite to the top boost. This reduces 
the momentum of the leptons causing more of them to fail p^-based selection criteria. 

Since W+ jet events dominate the background remaining after cuts, they are taken as the only 



background in this analysis. Fig. [32| shows the cumulative effect of cuts on the angular distribution of 
the charged lepton from W+ jets events. A peculiar feature of these events is evident in all of these 
distributions. This is the tendency for events to be grouped near cos Qi = 1. The events which populate 
this region tend to be the highest pt events. This shows that even basic jet and isolated-lepton definitions 
and pre-selection cuts bias the angular distribution of W+ jets events. 

When the event-selection criteria described in the previous sections are applied, the signal-to- 
background ratio (treating W+ jets as the only background) is found to be 2.6. Using the methods 
described earlier it is possible to estimate the polarisation of a mixed sample of t-channel signal and 
W+ jets background. The reference distributions for 100% spin-down and 100% spin-up top quarks 



mixed with background in a ratio of 2.6 are shown in Fig. [33J. Also shown is the angular distribution 
corresponding to a statistically-independent data sample with SM polarisation mixed with background 



50 




I 1£ 



14 
12 
1 

OjB 
0£ 
OA 
02 






4 




fwd jet 




+ 




+- 




++++ + 


: , , , , 1 , , , 


, i , , , , i , , , , 



-05 



05 1 

C03(Q ) 




Fig. 32: The effect of event selection cuts on the angular distribution of the charged lepton in Wjj events. The effects of the 
cuts are cumulative. The first distribution is the result of applying the pre-selection (trigger) cuts only. Further cuts are applied 
cumulatively from left-to-right, top-to-bottom. 

in the ratio 2.6. The \ 2 function presented in ( fH] ) is minimised to obtain an estimate of the polarisation 
of the top. To estimate the precision for one year of data-taking, the fit was done with 3456 signal events 
and 1345 background events, corresponding to 2 fb _1 of integrated luminosity (~ 1/5 of a year). For 
this integrated luminosity the error on the polarisation measurement is 4.0%. Then, assuming the statis- 
tics on the reference distributions, /£>(cos#£) and fjj{cos6p), will lead to a negligible source of error, 
this precision improves to 3.5%. Projecting these results to one year of data-taking at low luminosity 
(10 fb _1 ), assuming that the errors scale as the square root of the number of events, yields a predicted 
statistical precision of 1.6% on the measurement of the top polarisation. 



5.5 Conclusions on single top production 

As mentioned in the introduction, single-top quark production is the only known way to directly measure 
Vtb at a hadron collider. In this section we estimate the accuracy with which Vtb can be extracted at the 
LHC, and discuss what will be required to achieve that accuracy. 

There are four sources of uncertainty in the extraction of \ V t b\ 2 from the single-top cross section: 
theoretical, experimental, statistical, and machine luminosity. As we have seen, the statistical uncertainty 
with 30 fb _1 of integrated luminosity is less than 2% for both the f -channel process^] and associated 
production, and is 5.5% for the s-channel process (3% with 100 fb" 1 ). It will be a challenge to reduce 
the other sources of uncertainty to 5%, so we regard the statistical accuracy as being sufficient in all three 
processes. 

The traditional uncertainty in the machine luminosity is about 5% [159]. It may be possible to 
reduce the uncertainty below this value using Drell-Yan data, but this relies on accurate knowledge of 
the quark distribution functions. However, the process qq-^W* — >£p involves the identical combination 

20 Only 10 fb _1 are required to achieve this accuracy in the t-channel process. 
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Fig. 33: The first histogram shows the reference distribution for 100% spin-up top quarks after detector effects and event- 
selection criteria have been applied and the appropriate level of background has been mixed in. The second histogram shows 
the reference distribution for 100% spin-down top quarks. The third histogram represents the expected SM distribution for a 
statistically-independent sample of signal and background. 



of partem distribution functions as the s-channel process, so it can be used to almost directly measure the 
relevant parton luminosity, thereby avoiding the need to measure the machine luminosity [ |156[ ] . 

The theoretical uncertainty is under the best control in the s-channel process. The theoretical 
uncertainty is dominated by the uncertainty in the top quark mass; an uncertainty of 2 GeV yields an 
uncertainty of 5%. This is cut in half if the uncertainty in the top mass is reduced to 1 GeV. The small 
uncertainty due to variation of the factorisation and renormalisation scales can be reduced to a negligible 
amount by calculating the cross section at NNLO order, which should be possible in the near future. 
The small uncertainty from the parton distribution functions can be further reduced as described in the 
previous paragraph; this also obviates the need for a measurement of the machine luminosity. 

The theoretical uncertainty in the t-channel process is presently dominated by the factorisation- 
scale dependence and the parton luminosity. Although the scale dependence of the total cross section 
is small (4%), the uncertainty in the semi-inclusive cross section (a(p T i) < 20 GeV) is about 10%. 
This can be reduced by calculating the px spectrum of the b at NLO. It may also prove possible to 
measure the total cross section, although this has yet to be demonstrated. It is therefore plausible that 
the factorisation-scale dependence will be about 5% once the LHC is operating. It is also likely that the 
uncertainty from the parton distribution functions will be reduced below its present value of 10%. The 
parton luminosity could be directly measured using Wj production, which is dominated by gq-^Wq, 
and therefore involves the identical combination of parton distribution functions as the t-channel process. 
Again, this has the desirable feature of eliminating the need to measure the machine luminosity. 

The theoretical uncertainty in the associated-production cross section can be reduced far below 
its present value of 18%. A full NLO calculation should reduce the factorisation-scale dependence to 
roughly 5%. It is likely that the uncertainty from the parton distribution functions will also be reduced. 
Unless it is possible to measure the gg luminosity directly, the uncertainty from the parton distribution 
functions will be augmented by the uncertainty in the machine luminosity. 

As far as experimental systematic uncertainties are concerned, the extraction of a signal cross 
section requires knowledge of the backgrounds and of the efficiency and acceptance for the signal. These 
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analyses require hard cuts on both signal and background, and so the processes need to be modelled and 
understood very well. 

For all of these processes, the major backgrounds are tt and W+ jets. The largest background 
for the s-channel process (where a double 6-tag is employed) and associated production is tt. The tt 
process can be isolated in other decay modes and in principle well measured. In the t-channel process 
the biggest background comes from Wcj with the charm jet mistagged as a 6-jet. Obviously it would be 
worthwhile to develop a 6-tagging algorithm having greater rejection against such mistags, even at the 
cost of some signal efficiency, given that the signal rate is large. It may be possible to understand the 
W+jets backgrounds by comparing with a sample of Z+ jets events after applying similar selections to 
those used to select the single-top sample in W+ jets. The Z+ charm rate will be suppressed compared 
to the W+ charm rate since the latter is mostly produced from the strange sea, which is bigger than the 
charm sea; nonetheless, the cross section, kinematics, jet multiplicities and so on can all be compared to 
our simulations using the Z+ jets sample. 

The forward jet tag is very effective in enhancing the signal-to-background ratio in the t-channel 
process. This means that jets need to be found with good efficiency up to large rapidities, at least \rj\ ~ 4 
in the calorimeter. Unfortunately these observations also imply that the background estimate is very 
sensitive to the Monte Carlo predicting the correct mix of jet flavours and jet rapidities in the W+ jets 
events. (We note that VECBOS generates very few jets in the tagging region, and so far there is no 
collider data on forward jets in vector-boson events which could verify whether this is correct.) Of 
course, effort applied to understanding W+ heavy-flavour jets backgrounds will pay off in many other 
searches besides this one, and will be a very worthwhile investment. We also look forward to the results 
of ongoing efforts to improve the Monte Carlo simulation of vector-boson plus jet production [ ]160| ]. 
Requiring exactly two jets (as was done here to reject the tt background) also means that we will be very 
sensitive to our knowledge of jet efficiencies, QCD radiation, etc. The cross-section measurement also 
requires knowledge of the 6-tagging efficiency. This should be measurable at the few-percent level using 
control samples of tt events selected with kinematic cuts alone. 

As mentioned above, the purely statistical uncertainty in the cross-section measurement will be 
less than 5%, as will most of the theoretical uncertainties. It will be a considerable challenge to reduce 
the experimental systematic uncertainty to this level. At the present time, the experimental systematic 
uncertainty in the tt cross section at the Tevatron (which is a similar challenge in many respects, in- 
volving jets, 6-tagging, and background subtraction) is about 19% [[To[]. This total is made up of many 
components which are each at the 5% level, so while it will be a lot of work to reduce them, there is no 
obvious "brick wall" that would prevent this. 

Many of these systematic issues can also be addressed by comparing the t-channel and s-channel 
single-top processes. It will be a powerful tool to be able to measure V t b in two channels which have 
different dominant backgrounds, different selection cuts, and a different balance between theoretical and 
experimental systematic uncertainties. 

We are only just now entering the era of precision top physics with Run II at the Tevatron. Single- 
top production has not yet even been observed. We will learn a great deal over the next few years about 
how to model top events and their backgrounds, and how to understand the systematic uncertainties. The 
LHC will undoubtedly benefit from all this experience. 

If all sources of uncertainty are kept to the 5% level or less, it should be possible to measure 
| Vtb | 2 to 10% or less. We therefore regard the measurement of Vtb with an accuracy of 5% or less as an 
ambitious but attainable goal at the LHC. We have also seen that a measurement of the polarisation of 
single top quarks produced via the i-channel process will be possible with a statistical accuracy of 1.6% 
with 10 fb _1 . We have not attempted to estimate the systematic uncertainty in this measurement. 
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6. tt SPIN CORRELATIONS AND CP VIOLATION^] 

For ti production at the LHC quantities associated with the spins of the top and antitop quark will 
be "good" observables as well. The reason for this is well known. Because of its extremely short 



lifetime r t (see Section 2.1) the top quark decays before it can form hadronic bound states. Thus the 
information on the spin of the top quark does not get diluted. As the spin-flip time is much larger 
than Tt it is, moreover, very unlikely that the top quark changes its spin-state by emitting gluon(s) via a 
chromomagnetic dipole transition before it decays. In any case this amplitude is calculable with QCD 
perturbation theory. Hence by measuring the angular distributions and correlations of the decay products 
of t and i the spin-polarisations and the spin-spin correlations that were imprinted upon the ti sample by 
the production mechanism can be determined and compared with predictions made within the SM or its 
extensions. Therefore these spin phenomena are an additional important means to study the fundamental 
interactions involving the top quark. 

In this section we are concerned with the production and decay of top-antitop pairs. At the LHC 
the main tt production process is gluon-gluon fusion, qq annihilation being sub-dominant. As the main 
SM decay mode is t^W + b we shall consider here the parton reactions 

gg,qq->it + X->M + 4:f + X, (16) 

where / denotes either a quark, a charged lepton or a neutrino. If the final state in (|l^) contains two, 
one, or no high pt charged lepton(s) then we call these reactions, as usual, the di-lepton, single lepton, 
and non-leptonic tt decay channels, respectively. To lowest order QCD the matrix elements for (|l6|), 
including the complete ti spin correlations and the effects of the finite top and W widths, were given in 



[161, 162]. Spin correlation effects in tt production in hadron collisions were studied within the SM in 



[ |T62| , [T63i |16§ |T65|, [166|, |I67| |l68|l. 



In order to discuss the top spin-polarisation and correlation phenomena that are to be expected 
at the LHC it is useful to employ the narrow-width approximation for the t and i quarks. Because 
Tt/m t <C 1 one can write, to good approximation, the squares of the exact Born matrix elements M.^ x \ 
A = gg, qq, in the form 

| M (X) \ 2 oc Tr [pR^p] = Pa'aR ( al',^Pp'P- (17) 

The complete spin information is contained in the (unnormalised) spin density matrices for the 
production of on-shell ti pairs and in the density matrices p, p for the decay of polarised t and t quarks 
into the above final states. The trace in ( |l7| ) is to be taken in the t and i spin spaces. The decay density 
matrices will be discussed below. The matrix structure of R^ is 

Ralw = ^ (A) <w<w + + #f } <w(^w + 4 A Vw<yw > (is) 

where a 1 are the Pauli matrices. Using rotational invariance the "structure functions" , b\ ^ and 
can be further decomposed. A general discussion of the symmetry properties of these functions is 
given in [ ]169| ]. The function A^ x \ which determines the cross section, is known in QCD at NLO [38]. 
Because of parity (P) invariance the vectors B( A ) , B( A ) can have, within QCD, only a component normal 
to the scattering plane. This component, which amounts to a normal polarisation of the t quark, V\ , is 
induced by the absorptive part of the respective scattering amplitude and was computed for the above 



LHC processes to order o? s Ql70| ]. (V] = V\ if CP invariance holds.) The size of the normal polarisation 



depends on the top quark scattering angle and on the cm. energy. In the gluon-gluon fusion process V\_ 
reaches peak values of about 1.5%. In production at the LHC the polarisation of the top quark within 
the partonic scattering plane, which is P-violating, is small as well within the SM. Therefore the t and t 
polarisations in the scattering plane are good observables to search for P-violating non-SM interactions 



in the reactions (16) - see Section 3.4. 



21 Section coordinators: W. Bernreuther, A. Brandenburg, V. Simak (ATLAS), L. Sonnenschein (CMS). 
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The tt production by the strong interactions leads, on the other hand, to a significant correlation 
between the t and i spins. This correlation is encoded in the functions . Using P- and charge- 
conjugation (C) invariance they have, in the case of a tt final state, the structure 



C 



(A) _ JA) 



1 



Sij + c 



(19) 



where p and k t are the directions of flight of the initial quark or gluon and of the t quark, respectively, in 



the partem cm. frame. So far the functions dfi are known to lowest-order QCD only (see, e.g., [1164J]). 



For a ttX final state a decomposition similar to ( |19| ) can be made. 

From ( |l"9| ) one may read off the following set of spin-correlation observables [164] 



(k t ■ Si)(k r s(), 

(P • St)(p • Sf), 
St • Sf, 

(p • St) (Iff • s f ) + (p • Sf)(k t • s t ) 



(20) 

(21) 
(22) 
(23) 



where s t , Sf are the t and i spin operators, respectively. The observables (^0|), and (23) determine 
the correlations of different t, i spin projections. Eq. ( |20| ) corresponds to a correlation of the t and i 
spins in the helicity basis, whi le (|2~ij ) correlates the spins projected along the beam line. We note that the 
"beam-line basis" defined in [ |166| ] refers to spin axes being parallel to the left- and right-moving beams 
in the t and t rest frames, respectively. The tt spin correlation in this basis is a linear combination of 
©, ©, and ©. 

A natural question is: what is - assuming only SM interactions - the best spin basis or, equiva- 
lently, the best observable for investigating the ti spin correlations? For quark-antiquark annihilation, 
which is the dominant production process at the Tevatron, it turns out that the spin correlation (pTb 



[ |64j , [168J] and the correlation in the beam-line basis [|166j is stronger than the correlation in the helicity 



basis. In fact, for qq annihilation a spin-quantisation axis was constructed in [ |167[ ] with respect to which 
the t and i spins are 100% correlated. At the LHC the situation is different. For gg—>ti at threshold 
conservation of total angular momentum dictates that the tt is in a state. Choosing spin axes parallel 
to the right- and left-moving beams this means that we have t^t^ and t^t^ states at threshold. On the 
other hand at very high energies helicity conservation leads to the dominant production of unlike helicity 
pairs tjitL and t^tR. One can show that no spin quantisation axis exists for gg^ti with respect to which 
the t and i spins are 100% correlated. The helicity basis is a good choice, but one can do better. This 
is reflected in the above observables. Computing their expectation values and statistical fluctuations one 



finds [[164J that ( |22| ) has a higher statistical significance than the helicity correlation (gOJ) which in turn is 



more sensitive than (21 ) or the correlation in the beam-line basis. 

The spins of the t and t quarks are to be inferred from their P-violating weak decays, i.e., from 
t-^bW + -^>bl + or bqq' and likewise for i if only SM interactions are relevant. As already mentioned 
and used in previous sections, in this case the charged lepton from W decay is the best analyser of the top 
spin. This is seen by considering the decay distribution of an ensemble of polarised t quarks decaying 
into a particle / (plus anything) with respect to the angle between the polarisation vector £ t of the top 
quark and the direction of flight q/ of the particle / in the t rest frame. This distribution has the generic 
form 

' ^ + (24) 



T d cos I 



7 



where the magnitude of the coefficient Kf signifies the spin-analyser quality of /. The SM values for 
some /, collected from [171, 172, 173, 174 ], are given in Table 12. The corresponding t decay density 
matrix in the t rest frame is read off from (24) to be p a ' a = (11 + Kf a ■ q/) Q ' Q . The distributions for 
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Table 12: Correlation coefficient Kf for V — A charged current. In the last column l.e.j. stands for least energetic jet in the t 
rest frame. 





e + ,d,s 


V£,U,C 


b 


W + 


l.e.j. from qq' 




1 


-0.31 


-0.41 


0.41 


0.51 



the decay of polarised antitop quarks are obtained by replacing Kf — > — Kf in (24). The order a s QCD 
corrections to the decays t—>b£v and t—>Wb of polarised t quarks were computed in [171] and [175], 
respectively. For t, i polarisation observables these corrections are small. 

From the above table it is clear that the best way to analyse the tt spin correlations is through 
angular correlations among the two charged leptons £ + £'~ in the di-lepton final state. Using the produc- 
tion and decay density matrices in (17), neglecting the 1-loop induced QCD normal polarisation, and 
integrating over the azimuthal angles of the charged leptons one obtains the following normalised double 
distribution, e.g. in the helicity basis 

1 d 2 a 1 + Ck£+k#- cos 6 + cos 9- 



a dcos 9+dcos 



(25) 



where k#+K£- = —1 and is the angle between the tit) direction in the tt cm. frame and the 

direction of flight in the t(t) rest frame. The coefficient C, which is the degree of the spin 
correlation in the helicity basis, results from the q in (|9]) and it is related to [|165|] : 



c = N{t L i L + t R t R ) - N(t L t R + t R t L ) 
N{t L i L + t R t R ) + N(t L t R + t R t L ) ' 

For partonic final states and to lowest order in a s one gets C = 0.332 for the LHC. (The number 
depends somewhat on the partem distributions used. Here and below the set CTEQ4L [116] was used.) 
The optimum would be to find a spin axis with respect to which |C| = 1. But, as stated above, this is not 
possible for gg fusion. In addition to (25), analogous correlations among £ + from t and jets from i decay 
(and vice- versa) in the single lepton channels, and jet-jet correlations in the non-leptonic decay channels 
should, of course, also be studied. While the spin-analysing power is lower in these cases, one gains in 
statistics. 

From the above example is quite obvious that, for a given tt decay channel, the tt spin correlation 
will be most visible when the angular correlations among the t and t decay products are exhibited in 
terms of variables defined in the t and t rest frames. An important question is therefore how well the 
4-momenta of the t and i quarks can be reconstructed experimentally? We briefly discuss the results of 
a simulation of the single lepton and di-lepton channels [176] which includes hadronisation and detector 



effects using PYTHIA [ ]52[ ] and the ATLFAST Q105| ] software packages. The transverse momentum of 
every reconstructed object like a jet, a charged lepton, or the missing transverse energy of an event has to 
exceed a certain minimum value p™ n . The detector acceptances impose further restrictions on the phase 
space of the detected objects in pseudo-rapidity. 

In the case of the single lepton decay channels one isolated lepton (e or fj,^) is required. From 
the missing transverse energy of the event and the W mass constraint the longitudinal momentum p z 
of the neutrino can be determined up to a twofold ambiguity. It turns out that in most cases the lower 
solution of p z is the correct one. To complete the event topology, four jets are demanded. Two of them 
have to be identified as 6-jets coming from top decay. 

The two non-tagged jets are often misidentified due to additional activity in the detector from 
initial and final state radiation. To suppress the QCD background the invariant mass of the two jets has to 
lie in a narrow mass window around the known mass of the W boson. After this cut the two-jet system 
is rescaled to the W mass. Finally there is a twofold ambiguity when the 6-jets are combined with the 
reconstructed W bosons. The combination which yields the lower reconstructed top mass turns out to be 
the correct one most of the time. 
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Fig. 34: Joint distribution d 2 N/dcos8 + dcos9- generated Fig. 35: Same distribution as in the figure to the left, but 
with default PYTHIA. The detector response was simulated including the SM ti spin correlations. The detector response 
with CMSJET. was simulated with CMSJET. 

In the case of the di-lepton decay channels two isolated oppositely charged leptons are requested. 
Moreover two jets have to be detected and tagged as 6-jets. With the known top and W masses and with 
the missing transverse energy of the event the unknown 3-momenta of the neutrino and anti-neutrino can 
be computed using the kinematic constraints of the event. These result in a system of two linear and 
four quadratic equations. The equations can be solved numerically and usually several solutions arise. 
Since the experimentally determined momenta do not coincide with the corresponding variables at the 
parton level the kinematic constraints have to be relaxed somewhat in order to improve the reconstruction 



efficiency. The algorithm set up in [ |176[ ] was used to solve these equations. The best solution can be 



obtained by computing weights from known distributions. Following [176] the highest efficiency was 
obtained using the weight given by the product of the energy distributions of and V£ and the cos Q\ 
distribution in the ti cm. frame. 

For the LHC running at low luminosity (C = 10 33 cm _2 s _1 ), about 4 x 10 5 ti events per year 
are expected in the di-lepton decay channels (£ = e,(i). A further simulation of these channels was 

was used for the event generation, 
for the reconstruction of the i, i 



performed in order to study the joint distribution (|25j). PYTHIA 5.7 [| 



CMSJET [ ]177| ] for the detector response and the algorithm of [ |17f. 
momenta. The transverse momenta of the two isolated, oppositely charged leptons and of the two jets 
were required to exceed 20 GeV. The minimal missing transverse energy of the event was chosen to be 
40 GeV. A further selection criterion was that each jet provides at least two tracks with a significance 
of the transverse impact parameter above 3.0 to be tagged as 6-jet. The processes were simulated in 
two different ways. First the SM matrix elements of [75] for the reactions (|T6|), which contain the ti 
spin correlations, were implemented into PYTHIA. For comparison these channels were also simulated 
with the PYTHIA default matrix elements for gg, qq^-ti which do not contain spin correlations. In both 
simulations initial and final state radiation, multiple interactions, and the detector response was included. 
In Figs. 34, 35 we have plotted the resulting double distributions d 2 N/d cos 0+dcos They have been 



corrected for the distortions of the phase space due to the cuts. A fit to the distribution Fig. |35| according 
to (25) yields the correlation coefficient C = 0.331 ± 0.023, in agreement with the value C = 0.332 
obtained at the parton level without cuts. A fit to Fig. 34 leads to C = —0.021 ± 0.022 consistent with 
C = 0. Systematic errors, for instance due to background processes, e.g., Z*—>£ + l~ accompanied by 
two 6-jets, remain to be investigated. 
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From these double distributions one may form one- or zero-dimensional projections, for instance 



asymmetries as considered in [ ]16q , |165| , |168[ ]. Another approach is to study distributions and expectation 
values of angular correlation observables which would be zero in the absence of the ti spin correlations. 
A suitable set of observables is obtained by transcribing, for instance, the spin observables given above 
into correlations involving the directions of flight of those final state particles that are used to analyse 
the t and i spins. As an example we discuss the case of the single lepton channels t—>bqq', i^b£~vg. 
One may choose to analyse the t spin by the direction of flight of the 6-jet in the rest frame of the t 
quark and the t spin by the momentum direction q of the t~ in the laboratory frame. The latter is rather 



conservative in that no reconstruction of the t momentum is necessary. Then (|20|)-([22|) are translated into 
the observables 

Ox = (q^-p P )(q--p P ), (27) 

2 = (q£-k t )(q_-k t ), (28) 

03 = Q6-Q-, (29) 

where p p refers to the beam direction. The pattern of statistical sensitivities of the spin observables 
(PQ)-(p2|) stated above is present also in these angular correlations. Computing the expectation values 
(Oi) and the statistical fluctuations AO, and those of the observables for the corresponding charge 
conjugated channels, one gets for the statistical significances of these observables at the parton level 



f6l: Si « 0.007 y/N M -, S 2 ~ 0.025^7^, and S 3 as 0.055y/N M -, where N u - is the number of 



reconstructed events in the specific single lepton channel. The linear combination 

04 = 03 " Ox (30) 

has a still higher sensitivity than 0\, namely £3 0.073i/7V^. Even with 10 4 reconstructed b£~ and 
b£ + events each one would get a 7.3a spin-correlation signal with this observable. The significance of 
these observables after the inclusion of hadronisation and detector effects remains to be studied. 

The results of the above simulations are very encouraging for the prospect of t, i spin physics. On 
the theoretical side the NLO QCD corrections to the helicity amplitudes, and to the spin density matrices 
should be computed in order to improve the precision of the predictions and simulation tools. 

If ti production and/or decay is affected by non-SM interactions then the correlations above will 
be changed. One interesting possibility would be the existence of a heavy spin-zero resonance Xq (for 
instance a heavy (pseudo)scalar Higgs boson as predicted, e.g., by SUSY models or some composite 
object) that couples strongly to top quarks. For a certain range of masses and couplings to ti such an 



object would be visible in the tt invariant mass spectrum [74, 75]. Suppose one will be fortunate and 
discover such a resonance at the LHC. Then the parity of this state may be inferred from an investigation 
of tt spin correlations. This is illustrated by the following example. As already mentioned above, close 
to threshold gluon-gluon fusion produces a ti pair in a 1 Sq state. On the other hand if the pair is produced 
by the Xq resonance, gg-+Xo—>-ti, then for a scalar (pseudo-scalar) Xq the ti pair is in a s Pq (^o) 
state and has therefore characteristic spin correlations. Let us evaluate, for instance, the observable 
(E2). Its expectation value at threshold is (SfSj) = 1/4 (—3/4) if ti is produced by a (pseudo)scalar 
spin-zero boson, ignoring the gg^ti background. An analysis which includes the interference with the 
QCD ti amplitude shows characteristic differences also away from threshold. By investigating several 
correlation observables (i.e., employing different spin bases) one can pin down the scalar/pseudo-scalar 



nature of such a resonance for a range of Xq masses and couplings to top quarks [|75|]. 

Another effect of new physics might be the generation of an anomalously large chromomagnetic 
form factor k (see Section 7.1) in the it production amplitude which would change the spin correlations 



with respect to the SM predictions JT7^,[T79|] (see also [fiM pi]]). For the LHC with 100 ftr 1 integrated 



luminosity one obtains from a study of asymmetries (that were also used in | |179[ ]) at the parton level a 
statistical sensitivity of 5k ~ 0.02. 
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The top quark decay modes t^b£ + ue,bqq' might also be affected by non-SM interactions, for 
instance by right-handed currents or by charged Higgs-boson exchange, and this would alter the angular 
correlations discussed above as well. A Michel-parameter type analysis of the sensitivity to such effects 
at the LHC remains to be done. 

The large ti samples to be collected at the LHC offer, in particular, an excellent opportunity to 
search for CP-violating interactions beyond the SM in high energy reactions. (The Kobayashi-Maskawa 
phase induces only tiny effects in ti production and decay.) We mention in passing that such interactions 
are of great interest for attempts to understand the baryon asymmetry of the universe. Many proposals 
and phenomenological studies of CP symmetry tests in ti production and decay at hadron colliders have 



been made. The following general statements apply []169j] : A P- and CP- violating interaction affecting 
ti production induces additional terms in the production density matrices which generate two types 
of CP-odd spin-momentum correlations, namely 

k*-(st-st), (31) 

and 

kt • fa x sf) , (32) 

and two analogous correlations where k t is replaced by p. The longitudinal polarisation asymmetry ( pi] ) 
requires a non-zero CP- violating absorptive part in the respective scattering amplitude. In analogy to the 
SM spin correlations above, ( |3l| ) and (32) can also be transcribed into angular correlations among the t 



and t decay products, which may serve as basic CP observables (see below). 

As to the modelling of non-SM CP violation two different approaches have been pursued. One is to 
parameterise the unknown dynamics with form factors or, neglecting possible dependences on kinematic 
variables, with couplings representing the strength of effective interactions [ 180| , 182 , 173 , 183 , 178 , 179 



184 [185 1, and compute the effects on suitable observables. This yields estimates of the sensitivities to 
the respective couplings. For instance if ti production is affected by a new CP- violating interaction with 
a characteristic energy scale Ac p > V§ then this interaction may effectively generate a chromoelectric 
dipole moment (CEDM) dt of the top quark (see Section 7.1). Assuming 10 7 non-leptonic, 6x 10 6 single 
lepton, and 10 6 ti di-lepton events, the analysis of [185], using optimal CP observables, comes to the 



conclusion that a la sensitivity of 5(Re d t ) ~ 5x 10 g s cm may be reached at the LHC. A detector 



level study of CP violation in tt decays with di-lepton final states was performed in | |186[ ] . 

Alternatively one may consider specific extensions of the SM where new CP-violating interac- 
tions involving the top quark appear and compute the induced effects in ti production and decay, in 
particular for the reactions (|I6|). We mention two examples. In supersymmetric extensions of the SM, 
in particular in the minimal one (MSSM), the fermion-sfermion-neutralino interactions contain in gen- 
eral CP-violating phases which originate from SUSY-breaking terms. These phases are unrelated to the 
Kobayashi-Maskawa phase. The interaction Lagrangian for the top quark coupling to a scalar top t\^ 
and a gluino G reads in the mass basis 

C 6ti = iV2g s (e- i H L T l G a T a t l + e**i R T , l &T a t l ) + h.c, (33) 

1=1,2 

where g s is the QCD coupling. A priori the phase (fi t is unrelated to the analogous phases in the light 
quark sector which are constrained by the experimental upper bound on the electric dipole moment of 



the neutron. The CP-violating one-loop contributions of ( |33| ) to gg, qq—>tt were computed in [ |187| , |185| ]. 
A non-zero CP effect requires, apart from a non-zero phase 4> t , also non-degeneracy of the masses of 
t\ 2- For fixed phase and t\ — t% mass difference the effect decreases with increasing gluino and scalar 



top masses. Assuming the same data samples as in the CEDM analysis above, [185] concludes from 
a computation of optimal CP observables that a sensitivity \cp t \ > 0.1 can be reached at the LHC if the 
gluino and squark masses do not exceed 400 GeV. 
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Fig. 36: Right: differential expectation value of Qi as a function of the tt invariant mass at \[a = 14 TeV for reduced Yukawa 
couplings at = 1, at = — 1, and a Higgs boson mass m v — 400 GeV. The dashed line represents the resonant and the solid 
line the sum of the resonant and non-resonant ip contributions. Left: same as figure to the left, but for the observable Q2 [^ojl. 

Another striking possibility would be CP violation by an extended scalar sector manifesting itself 
through the existence of non-degenerate neutral Higgs bosons with undefined CP parity. Higgs sector CP 



violation can occur already in extensions of the SM by an extra Higgs doublet (see, for instance [188]). It 



may also be sizable in the MSSM within a certain parameter range [ ]189j ]. The coupling of such a neutral 
Higgs boson ip with undefined CP parity to top quarks reads 



Cy = -{\/2GF) l l 2 m t (a t tt + a t ti^t) ip , 



(34) 



where a t and a t denote the reduced scalar and pseudo-scalar Yukawa couplings, respectively (in the 
SM a t = 1 and a t = 0). The CP-violating effects of ([34]) on gg, qq—>ti were investigated for light 
ip in [190] and for <p bosons of arbitrary mass in [191, |169|] (see also [185, 86]). The exchange of 92 



bosons induces, at the level of the tt states, both types of correlations (32), (31) (the CP asymmetry 
ANlr = [iV(tifc) — N (t Rt r)]/ (all ti) considered in Jl90| ] corresponds to the longitudinal polarisation 
asymmetry (kf • (sf — s t ))). If the mass of ip lies in the vicinity or above 2m t the s-channel ^-exchange 
diagram gg^ip^ti becomes resonant and is by far the most important ip contribution. 

Simple and highly sensitive observables and asymmetries were investigated for the different tt 
decay channels in [|86|]. For the di-lepton channels the following transcriptions of (31) and (32) may be 
used: 

Qi = k t • q+ - k t - • q_ , (35) 



Q 2 = (k t - k f ) • (q_ x q+)/2 , 



(36) 



where k t , kf are here the t, t momentum directions in the tt cm. frame and q+,q_ are the £ + , 
momentum directions in the t and t quark rest frames, respectively. Note that Qi = cos 9 + — cos I 



where 6± are defined after (25). When taking expectation values of these observables the channels 
£ + , with £,£' = e,fi are summed over. The sensitivity to the CP- violating product of couplings 
7cp = —atat of heavy Higgs bosons is significantly increased when expectation values of (|^), ( f36| ) are 
taken with respect to bins of the tt invariant mass M 4 f. Two examples of these "differential expectation 



values" are shown in Fig. 36. In order to estimate the measurement errors we have used a sample of di- 
lepton events, obtained from a simulation at the detector level using the same selection criteria as in the 
simulation described above, and determined the resulting error on the expectation value of Q\, choosing 
Mff bins with a width of 10 GeV. With 2 x 10 5 reconstructed di-lepton events in the whole M t f range 
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we find that the error on (Qx) M . is slightly below 1% for a bin at, say, M t i = 400 GeV. In addition one 
may employ the following asymmetries which are experimentally more robust than (Qi): 

Am = N u ( Qi >0)-N u ( Qi <0) t (3?) 

where i = 1, 2 and Nu is the number of di-lepton events. From an analysis of these observables and 
asymmetries and analogous ones for the single lepton channels at the level of partonic final states the 



conclusion can be drawn [ |86| ] that one will be sensitive to |7cp| £ 0.1 at the LHC. This will constitute 
rather unique CP tests. 

7. TOP QUARK ANOMALOUS INTERACTIONS 22 

In the SM the gauge couplings of the top quark are uniquely fixed by the gauge principle, the structure 
of generations and the requirement of a lowest dimension interaction Lagrangian. Due to the large top 
mass, top quark physics looks simple in this renormalisable and unitary quantum field theory. Indeed, 

• the top quark production cross section is known with a rather good accuracy (~ (10 — 15) %), 

• there are no top hadrons (mesons or baryons), 

• the top quark decay is described by pure (V — A) weak interactions, 

• only one significant decay channel is present: t—>bW + (other decay channels are very suppressed 
by small mixing angles). 

This simplicity makes the top quark a unique place to search for new physics beyond the SM. If anoma- 
lous top quark couplings exist, they will affect top production and decay at high energies, as well as 
precisely measured quantities with virtual top quark contributions. 

We do not know which type of new physics will be responsible for a future deviation from the SM 
predictions. However, top quark couplings can be parametrized in a model independent way by an effec- 



tive Lagrangian. The top quark interactions of dimension 4 can be written (in standard notation [192]): 



V 1 q=d,s,b 



ry q=u,c,t 

plus the hermitian conjugate operators for the flavour changing terms. T a are the Gell-Mann matrices 
satisfying Tr (T a T b ) = 5 ab /2. Gauge invariance fixes the strong and electromagnetic interactions in (38) 



and hemiticity implies real diagonal couplings vf t ,af t , whereas the non-diagonal ones v tq , a t ' can 
be complex in general. Within the SM = a^ q = with Vtq the Cabbibo-Kobayashi-Maskawa 
(CKM) matrix elements, vf t = | — Isin 2 ^, af t = |, and the non-diagonal Z couplings are equal 
to zero. Typically modifications of the SM couplings can be traced back to dimension 6 operators in 



the effective Lagrangian description valid above the EW symmetry breaking scale [193, 194, 195] (see 



also [196, 132, 197]). Hence, they are in principle of the same order as the other dimension 5 and 



6 couplings below the EW scale. However, in specific models the new couplings in Eq. ( |38| ) can be 



large [ ]198| ]. Moreover, the present experimental limits are relatively weak and these couplings can show 
up in simple processes and can be measured with much better precision at the LHC. 

The dimension 5 couplings to one on-shell gauge boson, after gauge symmetry breaking, have the 



generic form: [199] : 



9 W 
A = -9s £ ^ta^T a (f? q + ihl ql5 )qG^-j= £ -f-ta^US + <^ W ^ 

q=u,c,t V q=d,s,b 



22 Section coordinators: F. del Aguila, S. Slabospitsky, M. Cobal (ATLAS), E. Boos (CMS). 
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" e E -fia^(f? g + ihl ql5 ) q A^ 



9 



q=u,c,t 



2 COS Ow 



E 

q=u,c,t 



J [IV 



(39) 



plus the hermitian conjugate operators for the flavour changing terms. is d^G® 



dyG ^ (see, how- 
ever, below) and similarly for the other gauge bosons. We normalise the couplings by taking A = 1 TeV. 
k is real and positive and /, h are complex numbers satisfying for each term |/| 2 + |/i| 2 = 1- As in 
the dimension 4 case these dimension 5 terms typically result from dimension 6 operators after the EW 
breaking. They could be large, although they are absent at tree level and receive small corrections in 
renormalizable theories. At any rate the LHC will improve appreciably their present limits. 

There are also dimension 5 terms with two gauge bosons. However, the only ones required by 
the unbroken gauge symmetry SU(3)c X U (1)q, and taken into account here, are the strong couplings 
with two gluons and the EW couplings with a photon and a W boson. They are obtained including 
also the bilinear term g s f abc G b ^Gl, with f abc the SU(3)c structure constants, in the field strength G^ u 
in (39) and the bilinear term —ie(A^W^~ — A^W^) in W+,, respectively. We do not consider any other 



dimension 5 term with two gauge bosons for their size is not constrained by SU(3)c x U(1)q and/or 
they only affect to top quark processes with more complicated final states than those discussed here. We 
will not elaborate on operators of dimension 6, although the first q 2 corrections to dimension 4 terms 



could be eventually observed at large hadron colliders [|134j]. In this section we are not concerned with 
the effective top couplings to Higgs bosons either. 

In what follows we study the LHC potential for measuring or putting bounds on the top quark 
anomalous interactions in (| 



f), ( |39| ) through production processes. Results from top quark decays are 
presented is Section [8] The ti couplings to gluons are considered first, since they are responsible for 
tt production. Secondly we discuss the top quark couplings tbW. In the SM this coupling is not only 
responsible for almost 100% of the top decays but it also leads to an EW single top production mode, 
as reviewed in Section EH. Finally we deal with the t flavour changing neutral currents (FCNC). The yti 
and Ztt vertices have not been considered here because e + e~ and colliders can give a cleaner 

environment for their study. 



With the exception of the summary Table 23, we will quote limits from the literature without 



attempting to compare them. In Table |13| we illustrate statistics frequently used and which we will refer 
to in the text when presenting the bounds. As can be observed, the number of signal events, and the limit 
estimates, vary appreciably with the choice of statistics. We do correct for the different normalizations 
of the couplings used in the literature. 



Table 13: Limits on the number of signal events S obtained with different statistics. B is the number of background events. In 
the other columns we gather S for (1): 99% CL (3 a) measurement, > 3; (3): 99% CL (3 a) limit, -jjL > 3; and (5): 

99% CL for the Feldman-Cousins (FC) statistics [Eoofa; and similarly for (2), (4), and (6), for the 95% CL (1.96 a), respectively. 



B 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 





9 


3.84 








4.74 


3.09 


5 


12.57 


6.71 


6.71 


4.38 


8.75 


6.26 


10 


15 


8.41 


9.49 


6.20 


10.83 


7.82 


15 


16.96 


9.75 


11.62 


7.59 


12.81 


9.31 



7.1 Probes of anomalous gtt couplings 



The combination ^K 9 tt ftt ( see ©) 

can be identified with the anomalous chromomagnetic dipole mo- 
ment of the top quark, which, as is the case of QED, receives one-loop contributions in QCD. Therefore, 
its natural size is of the order of a s /ir. As we observed above, when this coupling is non-zero a direct 
ggit four-point vertex is induced as a result of gauge invariance. 
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Table 14: Attainable la limits on Re(d t ) and Im(d t ) 
luminosity (10 fb _1 ) 



through T33, Ae and Q33 for one year of the LHC running at low 



[204] 



Observable 


Attainable la limits 


T33 
A E 
Q33 


\Re(d t )\ = 0.899 x Kr^&cm 
\Im(d t )\ = 0.858 x 10- ls g s cm 
\Im(dt)\ = 0.205 x 10- 17 5 s cm 



On the other hand the combination ^^Kf t hf t can be identified as the anomalous chromoelectric 



dipole moment of the top quark. Within the SM this can arise only beyond two loops [ [20 1| ]. On the 
other hand it can be much larger in many models of CP violation such as multi-Higgs-doublet models and 



SUSY fl202| ]. Therefore, such a non- vanishing coupling would be a strong indication of BSM physics. 
Considering the gluonic terms in (38), (|39|) for the process of light quark annihilation into it one 



obtains Q18l| , |203[ ] 

da, 



•qq 



dt 



~9W 



/3 2 (1 



8m t 
A 



(ftt + ftt) + 



32mf 



)VfJ 2 + ^(4) 2 /3 2 (i 



(40) 

s being the incoming parton total energy squared, z being the cosine of the scattering angle 8* in the cms 
of the incoming partons, and (5 - 



1 — Amf/s. 



The squared matrix element for gg annihilation is a more complicated expression; we refer to | J 1 8 1 



204] ] for exact formulas. If the (anomalous) couplings are assumed to be functions (form-factors) of q 
and then corrected by operators of dimension higher than 5, the gg annihilation amplitude would be eval- 
uated at different scales (for the t(u) and s channels), and an additional violation of the SU(3)c gauge 
invariance could be made apparent. For a detailed discussion of this problem see, for example, [ |181| ] and 
references therein. 

The effects associated with K^ t ff t were examined in [181, 205, [206| ]. As shown in [134] they will 
be easily distinguishable from the effects of q 2 corrections to the strong coupling due to operators of 



dimension 6, which are relatively straightforward to analyse [ ]195| ] in tt production since the effective 
coupling would be a simple rescaling of the strength of the ordinary QCD coupling by an additional (in- 



dependent amount. It was shown in [ 206 ] that the high-end tail of the top quark pt and M tf distributions 



are the observables most sensitive to non-zero values of n 9 tt ff t , with a reach for k = ^x^/vltftt as small 



as ~ 0.03. For these values of k, only a minor change in the total tt rate is expected (see Fig. |37| ). 
The effect of a non-zero tz g tt h 9 tt was analysed, in particular, in [204, 207, 172]. It was shown in [204] 
that information on nf t hf t could be obtained by studying the following correlation observables between 
l + £~ lepton pairs produced in in di-lepton decays: 



T 33 = 2(pi - p e ) 3 (Pi x p e ) 3 , 

A E = E- t -E e , Q33 = 2 (jpi + Pe ) s ( Pi - p e ) 3 



'M 



p\) 



Table 14 shows the 1 a sensitivities of these correlations to Re(dt) and Im(d t ) (where, dt 
Quantitatively, T33 and Q33 enable us to probe Re(dt) and Im(dt) of the order of 10~ 17 
tively, and Ae allows us to probe Im(dt) down to the order of 10~ 18 g s cm (see [204] for details). 



g s cm, respec- 



7.2 Search for anomalous Wtb couplings 

The Wtb vertex structure can be probed and measured using either top pair or single top production 
processes. The total ti rate depends very weakly on the Wtb vertex structure, as top quarks are dom- 
inantly produced on-shell [208]. However, more sensitive observables, like C and P asymmetries, top 
polarisation and spin correlations provide interesting information, as discussed in Section ^] The single 
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K 



Fig. 37: Cross section for tt production (solid) at the LHC as a function of k. The part of the cross section arising from the 



j(qq) annihilation is shown by the dash-dotted (dotted) curve (see [ 206 1 for details) 



top production rate is directly proportional to the square of the Wtb coupling, and therefore it is poten- 
tially very sensitive to the Wtb structure. In single top events the study of the top polarisation properties 



potentially provides a way to probe a Wtb coupling structure [J209J] . The potential to measure anomalous 
Wtb couplings at LHC via single top from the production rate and from kinematical distributions has 
been studied in several papers [195, 210, 135 , |3( 



In the model independent effective Lagrangian approach [|193|, 194, [1951] there are four indepen 



dent form factors describing the Wtb vertex (see Q195| ] for details) 
unitary gauge [211, 208, 135] is given in (38), (39). As already mentioned the (V 



The effective Lagrangian in the 
A) coupling in the 

SM carries the CKM matrix element V t b which is very close to unity. The value of a (V + A) coupling 
is already bounded by the CLEO b — > data [212, 213] at a level [195, 213] such that it will be out 
of reach even at the high energy colliders. Since we are looking for small deviations from the SM, 

and an analysis is presented only for the two 



K 



Hb 

(-/, 



in the following and will be set to v ^ 
'magnetic' anomalous couplings Fi2 = 2M ^ 

Natural values for the couplings \F L ^ 2 \ are in the region of 
unitarity violation bounds for \F 1 



l 

2 

W* 
tb 



F R2 



2M W ,„H/ 



tb 



yjm h m, t 



0.1 ni9q] and do not exceed the 



L(R)2\ 



0.6 [194]. 



Calculations of the complete set of diagrams for the two main processes pp — > bbW and pp 



bbW + jet have been performed Q135Q for the effective Lagrangian in (|38|), (J39|), using the package Com 



pHEP [214]. The calculation includes the single-top signal and the irreducible backgrounds. Appropriate 
observables and optimal cuts to enhance the single-top signal have been identified through an analysis of 
singularities of Feynman diagrams and explicit calculations. The known NLO corrections to the single 



top rate [126, 125] have been included, as well as a simple jet energy smearing. The upper part of Fig. 



presents the resulting 2 a exclusion contour for an integrated luminosity of 100 fb , assuming e, /x and 
t— >£ decays of the PF-boson. The combined selection efficiency in the kinematical region of interest, 



including the double 6-tagging, is assumed to be 50%. Figure 38 demonstrates that it will be essential to 
measure both processes pp — > bbW and pp — » bbW + jet at the LHC. The allowed region for each single 
process is a rather large annuli, but the overlapping region is much smaller and allows an improvement 
of the sensitivity on anomalous couplings of an order of magnitude with respect to the Tevatron. Since 
the production rate is large, even after strong cuts, expected statistical errors are rather small, and the 
systematic uncertainties (from luminosity measurements, parton distribution functions, QCD scales, m t , 
. . . ) will play an important role. As it is not possible to predict them accurately before the LHC startup, 



we show here how the results depend on the assumed combined systematic uncertainty. Figure 38 (lower 
part) shows how the exclusion contours deteriorate when systematic errors of 1% and 5% are included. 
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Fig. 38: Limits on anomalous couplings after optimised cuts from two processes pp — > bbW and pp 
plot). Dependence of the combined limits on the values of systematic uncertainties (lower plot). 



bbW + jet (upper 



Note that a systematic error of 10% at the LHC will diminish the sensitivity significantly and the allowed 
regions will be comparable to those expected at the upgraded Tevatron. 

The rate of single top production at LHC is different from the rate of single anti-top production. 
This asymmetry provides an additional observable at LHC that is not available at the Tevatron and which 
allows to reduce systematic uncertainties. 

The potential of the hadron colliders can be compared to the potential of a next generation e + e~ 



linear collider (LC) where the best sensitivity could be obtained in high energy 7e-collisions Q208| , |215| ]. 
The results of this comparison are shown in Table |l^. From the table we see that the upgraded Tevatron 
will be able to perform the first direct measurements of the structure of the Wtb coupling. The LHC 
with 5% systematic uncertainties will improve the Tevatron limits considerably, rivalling with the reach 
of a high-luminosity (500 fb^ 1 ) 500 GeV LC option. The very high energy LC with 500 fb^ 1 luminosity 
will eventually improve the LHC limits by a factor of three to eight, depending on the coupling under 
consideration. 
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Table 15: Uncorrected limits on anomalous couplings from measurements at different machines. 





F L 2 


F R2 


Tevatron (A sys 


. « 10%) 


-0.18 


... +0.55 


-0.24 


. . . +0.25 


LHC (A sys . ~ 


5%) 


-0.052 


. . . +0.097 


-0.12 


... +0.13 


7e (yJSe+e- = 


0.5 TeV) 


-0.1 


... +0.1 


-0.1 


... +0.1 


Qs e +e- = 


2.0 TeV) 


-0.008 


. . . +0.035 


-0.016 


... +0.016 



7.3 FCNC in top quark physics 

In the previous subsections, we analysed top quark anomalous couplings as small deviations from the 
ordinary SM interactions (gti and tWb vertices). Here we consider new processes which are absent at 
tree-level and highly suppressed in the SM, namely the FCNC couplings tVc and tVu (V = g, 7, Z). 
The SM predicts very small rates for such processes [216] (see Table ^). The top quark plays therefore 
a unique role compared to the other quarks, for which the expected FCNC transitions are much larger: 
the observation of a top quark FCNC interaction would signal the existence of new physics. As an 
illustration, Table 16 shows predictions for the top quark decay branching ratios evaluated in the two- 
Higgs doublet model [217], the SUSY models [218], and the SM extension with exotic (vector-like) 
quarks [ll98t|. 



Table 16: Branching ratios for FCNC top quark decays as predicted within the SM and in three SM extensions. 





SM 


two-Higgs Q217P 


SUSYQ218H 


Exotic quarks [198] 


B 

B(t^qZ) 


5 x lO" 11 
5 x 10~ 13 

~ 10- 13 


~ io- & 

~ 10~ 7 
~ 10" 6 


~ 10~ 3 
~ 10~ 5 

~ 10- 4 


~ 5 x 10" 4 
~ 10~ 5 

~ 10- 2 



In the effective Lagrangian description of (pq), ( p9| ) it is straightforward to calculate the top quark 
decay rates as a function of the top quark FCNC couplings: 

r(t- 97 ) = 



v(t- 
r(t- 

T(t- 



■qZ\ 




-a s m 3 



I Z\2 , I Z\2 
\ u tq\ I" 



Hq\ 



1 




a 



1 



4M| sin 2 26 w 



a m+ 



sin 29w 



Ml 
m? 




(41) 



-f 1 + , (42) 



(43) 



For comparison, Table In] collects the rare top decay rates normalised to k 9 



tq 



nf q = 1, and for the SM. We assume m t = 175 GeV, A = 1 TeV, a = j^, a s = 0.1 and sum the 



"tq 



,,Z\2 
hq\ 



+ a 



Z |2 
tq\ 



decays into q = u,c. In this 'extreme' case with the anomalous couplings equal to one the top can decay 
into a gluon or a Z boson plus a light quark q = u,c and into the SM mode bW at similar rates. 



7.31 Current Constraints on FCNC in top quark physics 

Present constraints on top anomalous couplings are derived from low-energy data, direct searches of top 
rare decays, deviations from the SM prediction for it production and searches for single top production 
at LEP2. 

Indirect constraints: The top anomalous couplings are constrained by the experimental upper 
bounds on the induced FCNC couplings between light fermions. For example, the 7^ term in the Ztq 
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Table 17: Top quark decay widths and corresponding branching ratios for the anomalous couplings equal to one and for the 
SM. In the fourth line we gather the values of the corresponding anomalous couplings giving the same decay rates as in the SM. 





Top decay mode 

W + b (c + u)g (c + u)7 (c + u)Z 1 (c + u)Z a 


FCNC coupling 




1 


1 


1 


1 


r(GeV) 


1.56 


2.86 


0.17 


2.91 


0.14 


B 


0.20 


0.37 


0.022 


0.38 


0.018 


FCNC coupling 




8 x 10" 5 


3 x 10" a 


4 x 10~ y 




r SM (GeV) 


1.56 


8 x 10 -11 


8 x 10~ 12 


2.2 x 10~ 13 




E>SM 


1 


5 x 10 -11 


5 x 10~ 13 


1.5 x 10~ 13 





vertex generates an effective interaction of the form [219] 



C 



eff 



fj 



cos 9w 



aijfil' 



1 



75 



fjZn + h.c, 



where are two different light down-type quarks. The one-loop estimate of the vertex gives: 



1 m 2 r 
16-7T 2 v 2 



VuK + <)V Q3 + V*(vf q * + af q *)V tj 



In 



A 2 



m 



2 - 



(44) 



(45) 



where V{j are the CKM matrix elements. Then, using the results of Q219j ] and the experimental constraints 
from [ |192[ ] on Kl—*[j, + [j,~, the K^-Kg mass difference, B° — B° mixing, B^l + l~X and 6^57, one 
obtains: 



a sd < 2 x 10 5 , a bd < 4 x 10 4 , a bs < 1.4 x 10 3 , 
and, taking v = 250 GeV, m t = 175 GeV and A = 1 TeV: 

\v? u + a?u\ < 0-04, \vg + ag\ < 0.11. 



(46) 



(47) 



J tq 



af q do not contribute to aij for massless external fermions. However, both chiralities of the Ztq 



vertex contribute, for instance, to the vacuum polarization tensor Ii iJ,u {q 2 ). Thus, using the recent value 
for the p parameter, p = 0.9998 ± 0.0008 (+0.0014) [ ]i92T l, the following 2a limit is obtained: 



vt q \ 2 + K\ 2 < 0.15. 



(48) 



CDF results: The CDF collaboration has searched for the decays t— >^c{u) and t— >Zc(u) in the 
reaction pp — > itX at y/s = 1.8 TeV, obtaining the following 95% CL limits JT3|]: 



BR(t^c7) + BR(t^« 7 ) < 3.2%, BR(t->cZ) + BR(t^uZ) < 33% . 
These translate into the bounds on the top anomalous couplings 



(49) 



kL < 0-78, 



v t9 ^w., u , y,^| 2 + |af ? | 2 < 0.73. (50) 

ti production via FCNC: Constraints on the vertex gtq can be derived form the study of the tt- 
pair production cross-section. Imposing that the tt-pair production cross-section, including the possible 
effect of anomalous couplings, should not differ from the observed one (assumed in this study to be 



a 



exp 



6.7 db 1.3 pb [pp) by more than 2 pb, leads to the constraint [220]: 



K 



■P- < 0.47 TeV -1 . 
A ~ 



(51) 
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Table 18: Short summary of the LEP 2 results for e e — ► tq. The theoretical value a t h is evaluated assuming the limit on 
the corresponding anomalous coupling in (p0[). 



Collab. 


^(GeV) 




a exp (95%CL) 




DELPHI 

ALEPH 

DELPHI 


183 GeV [|222 
189 GeV [Em 
189 GeV [|223 


] 
] 
] 


47.7 
174 
158 


< 0.55 pb 

< 0.60 pb 

< 0.22 pb 


< 0.15 pb 

< 0.30 pb 

< 0.30 pb 



=5 

N 

t 
-(-» 

+ 

N 

T 
■»-> 

$-1 



0.1 







DELPHI upper limit (preliminary) 



0.3 



0.2t 







CDF i 


lpper limit 
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Fig. 39: Upper limit on branching fraction of t^Zq resulted from LEP 2 data. Dashed curve corresponds to kJ. = 0, while 



solid one corresponds to < 0.78. 



FCNC at LEP 2: Since 1997, LEP2 has run at cms energies in excess of 180 GeV, making the 
production of single top quark kinematically possible through the reaction: 



7 *(Z*) 



(52) 



Two LEP experiments p21| , [222], [223J have presented the results of their search for this process. A short 
summary of these data is given in Table Il8l The production cross section is very sensitive to the top 
quark mass, at q ~ (1 — ^-) 2 (see [224] for details). Therefore, the upper limit on the corresponding 
branching ratio depends from the exact value of m t as well, as shown in Fig.B^. The current constraints 



on the top quark FCNC processes are summarised in Table 19. Note that the LEP2 limit is slightly 
better then that given by CDF (|f 
highest-energy runs are analysed. 



These constraints should further improve once the data from the 



7.4 Search for FCNC in top quark production processes 

FCNC interactions of top quarks will be probed through anomalous top decays (as discussed in Sec- 
tion |8]), and through anomalous production rates or channels, as discussed in the remainder of this 
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Table 19: Current constraints on top quark FCNC interactions. 



t- 


^91 


BR< 17% 


g 

K tq 


< 0.47 


(other FCNC coupling 


is zero) 


t- 




BR< 3.2% 


7 
K tq 


< 0.78 


(other FCNC coupling 


>s zero) 


t- 


+Zq 


BR< 22% 


M\ 2+ 


\a? q \ 2 < 0.55 


(other FCNC coupling 


>s zero) 



Table 20: Upper bounds on the anomalous couplings tif u and nf c from single top production processes. The symbols 2— >1 



and 2^2 correspond to the reactions quark-gluon fusion, and single top production, respectively [225, 142]. 





Tevatron 
Run 1 Run 2 Run 3 


LHC 


V~s (TeV) 
^(fb- 1 ) 


1.8 2.0 2.0 
0.1 2 30 


14.0 
10 


4(2-1) 
4(2-2) 

4(2-1) 
4(2-2) 


0.058 0.019 0.0092 
0.082 0.026 0.013 
0.22 0.062 0.030 
0.31 0.092 0.046 


0.0033 
0.0061 
0.0084 
0.013 



section. 



7.41 Deviations from SM expectations for tt production 

As shown in the previous subsection, the FCNC tgg-vertex contributes to gg—>tt transitions, and to a 



possible enhancement of the top quark production at large E t and M t p A recent study [ |220| ] shows that 
at the LHC the sensitivity to these couplings is equivalent to that found with the data of Run 1 at the 
Tevatron: 




LHC 




0.5 TeV" 



(53) 



FNAL 



7.42 'Direct' top quark production (2— >1) 



The 'quark-gluon' fusion process [ ]225| ] g + u(c)—>t is characterised by the largest cross-section for top 
quark production through FCNC -interactions assuming equal anomalous couplings. At the LHC, using 

a 

1 TeV -1 are equal to: 

<r(cg-*t) ~ 6 x 10 3 pb. (54) 



the CTEQ2L structure functions [ |115| ], these cross sections for 

a{ug-+t) ~ 4 x 10 4 pb , cr(ug^t) ~ 1 x 10 4 pb 



Note that a(ug—>t) is about 50 times larger than the SM cross section. The major source of background 
to this is the W+ jet production. The additional background due to single top production, when the 
associated jets are not observed, should not exceed 20% of the total background and was therefore 
ignored. To reproduce the experimental conditions, a Gaussian smearing of the energy of the final leptons 



and quarks was applied (see [225] for details). Cuts on the transverse momentum (pr > 25 GeV), 
pseudo-rapidity (\r/j\ < 2.0, |%| < 3.0), and lepton-jet separation (AR > 0.4) were applied. A b- 
tagging efficiency of 60% and a mistagging probability of 1 % were assumed. 

The criterion S/ y/S + B > 3 was used to determine the minimum values of anomalous couplings. 
The couplings tgu and tgc have been considered separately. The resulting constraints on k\ u 
given in Table EO, which also contains the results of an analysis done for the Tevatron. 



and k\ c are 
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qq — > tc 



eg -> tg 






tc 



Z( 7 ) c/ -HZ( 7 ) 
Fig. 41: s-channel diagrams for (V = Z, 7) production 



Fig. 40: 2^2 single top quark production. 
7.43 Single top quark production (2^2) 



Single top quark production in 2— >2 processes has been studied as well [ |142| ]. There are four different 
subprocesses, which lead to one top quark in the final state together with one associated jet (see Fig. ^ 
and []142|] for detailed considerations): 



qq^tq, gg^tq, qq^tq, qg^tg 



(55) 



The major background comes from W + 2 jets and W + bb production, as well as from single top 
production. In addition to the cuts and tagging rates used in the above analysis of 'direct' top production, 
additional cuts on the reconstructed top mass (145 GeV < M\jw < 205 GeV), on p^j > 35 GeV, 
and on jet-jet and lepton-jet separation (ARjj > 1.5, ARu > 1.0) were applied here to improve 
the signal/background separation. The corresponding limits on anomalous couplings in the top-gluon 
interaction with corn quarks are given in Table 



7.44 tZ and production 
All the anomalous couplings may contribute to the processes q g —>t Z(-y), and were considered in [226 



227 ]. The left diagram in Fig. ^T] corresponds to the Z(-j)tq coupling, while the right one shows the top- 
gluon anomalous coupling (the corresponding f-channel diagrams are not shown). For all the calculations 
presented here, the MRSA PDF set [|228|] with Q 2 = s was used. The resulting total cross sections for 



Hq 



"tql 



1 are [227] 



a(ug^Z t) 



73 pb, 
746 pb, 



a{cg^jt) 
v(cg->Zt) 



10 pb, 
114pb. 



Different background sources (W+ jets, Z+ jets, ZW+ jets, Wbb+ jets, ti, and Wt production) were 
considered. The experimental conditions were simulated by a Gaussian smearing of the lepton, photon 
and jet energies (see [ ]227| ] for details). Cuts on the transverse momenta, pt{£, j, 7) > (15, 20, 40) GeV, 
on pseudo-rapidities, |r/^ 7 | < 2.5, and on lep ton-jet-photon separation (AR > 0.4) were applied. A b- 
tagging efficiency of 60% and a mistagging probability of 1% were assumed. It was found that 6-tagging 
plays an essential role in tracing the top quark and reducing backgrounds. 

It has been shown that the best limits on the top quark FCNC couplings can be obtained from the 



decay channels Zt—>£ + £ £vb and -yt— >-y £ub (see [226] and [227] for details). Upper bounds at 95% 
CL are derived using the FC statistics [200]. Table 21 collects the corresponding limits on eight top 



anomalous couplings. Like in previous cases the bounds on u and c couplings were obtained under the 
assumption that only one anomalous coupling at a time is non-zero. The analysis was done for both 
Tevatron and LHC but with different optimized cuts. 
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different normalizations used in [ 226 , 227 ] 



Table 21: Upper bounds on top anomalous couplings (see (B8p9()) from Zt and jt production. We have corrected for the 





Tevatron 
Run 1 Run 2 


LHC 


(TeV) 
C (fbr 1 ) 


1.8 2.0 
0.1 2 


14.0 14.0 
10 100 


K 9 
K 9 

K tc 

7 
K tu 

7 
K tc 

\J\vtu\ 2 + \afu\ 2 
\/Wcl 2 + |af c l 2 

K tc 


0.31 0.057 
0.86 0.18 

0.49 0.13 
1.71 0.43 


0.0097 0.0052 
0.020 0.011 
0.013 0.0060 
0.037 0.018 

0.016 0.0078 

0.032 0.016 
0.040 0.018 
0.097 0.046 




Fig. 42: Diagram describing like-sign top quark pair production 
7.45 Like -sign tt (ii) pair production 

Additional evidence for a FCNC gtq coupling can be sought through the production of like-sign top pairs 
(see Fig. 42). 



pp->ttX, pp^ttX (56) 
The ATLAS collaboration performed a detailed investigation of this reaction for the case of high lumi- 



nosity, d n t = 100 fb 1 (see [ f30| ] and [220] for details). All the three anomalous couplings contribute to 
this process and the kinematics of the it-pair is almost the same as for the conventional it-pair production. 

An experimentally clean signature of tt (it) production is the production of like-sign high pr 
leptons plus two hard 6-jets. The main sources of background are qq'—>-tiW and qq—i-W^q'W^q'. The 
expected cross sections for the signal (with K 9 g = K~J q = \v z c \ 2 + \af c \ 2 = 1) and background processes 
are equal to: 

cr(it) = 1920 pb, a(tt) = 64 pb, 

a(W + ti) = 0.5 pb, a(W~ti) = 0.24 pb, 

a(W+W + qq) = 0.5 pb, a(W~W~qq) = 0.23 pb. 



CTEQ2L structure functions [ |115| ] were used with the evolution parameter Q = mf for the signal 
and Q 2 = for the background calculations. PYTHIA 5.7 [Q| was used for the fragmentation and 
all events were passed through the ATLFAST detector simulation. An additional reducible like-sign 
di-lepton background is due to ti events with a b semi-leptonic decay. The initial selection required 
therefore two like-sign isolated leptons with pt > 15 GeV and \r/\ < 2.5 as well as at least two jets with 
Pt > 20 GeV and \r]\ < 2.5. In order to get a better signal/background separation jets with pt > 40 GeV 



(with at least one tagged as a fc-jet) were required (see [ ]30| , [22CJ| ] for other cuts). The potential reach of 



this study, using the S/ y/ S + B > 3 criterion, is given in Table I 
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Table 22: The limits on anomalous couplings from an improved ATLAS analysis [|3C|, [22C| ] of like-sign top-pair production at 
the LHC for the case of high luminosity, Ant = 100 fb -1 . The contribution from the term in the Ztq vertex is ignored. 



g 
K tu 


g 

K tc 


7 
K tu 


7 
K tc 


\f\vtu\ 2 + \atu\ 2 


\f\v? c \ 2 + \ag\ 2 


0.078 


0.25 


0.14 


0.32 


0.27 


0.85 



Table 23: Summary of the LHC sensitivity to the top quark anomalous couplings K^ q , n] q and yj\vf u \ 2 + \af u \ 2 . The resulting 
constraints are presented in terms of 'branching ratio', T(t—>qV)/TsM(= 1-56 GeV). The results for the Tevatron option are 
also given (see text for explanation). 2— >1, 2^2, tV, and tt stand for quark-gluon fusion, single top production, t + ^(Z) 
production, and like-sign top-pair final states, respectively. The 'decay', 'ATLAS', and 'CMS' labels denote the results obtained 
from the study of top decay channels, documented in Section |i] 





Tevatron 


LHC 


V^(TeV) 


1 8 
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30 


10 


1zt 
100 
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3.1 x 10~ 3 

6.2 x 1(T 3 
1.8 x 1CT 1 


3.3 x 1CT 4 
6.2 x 1CT 4 
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1.5 x 10~ 4 

2.7 x 10" 3 


1.0 x 10~ b 
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1.7 x 10~ 4 

1.5 x 10" 2 


3.2 x 10~ b 
1.1 x 10" 5 
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2^1 [ 
2^2 [ 

tv | 

decay [ 
tt | 
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142 
226 
229 
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227 
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1.9 x 10~ 2 
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2.0 x 10~ 3 

2.7 x 10~ 3 


6.5 x 10" 5 
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1.6 x 10" 1 
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decay [ 
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1 
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1.9 x 10~ 4 
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tv | 

tt [ 
ATLAS [ 
CMS [ 
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220 
30] 
230 


in 

3j 


I 
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1.5 x 10" 4 
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1.9 x 10" 4 
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tt [ 
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30] 
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3fJ] 


tuZ 


4.5 x 1CT 1 


3.2 x 10" 2 
1.1 x 10~ 2 


5.2 x 10~ 3 


4.8 x 10" 4 

5.8 x 10~ 4 

1.9 x lO" 1 
6.5 x 10~ 4 
1.4 x 10" 3 


1.1 x 10" 4 
1.9 x 10~ 4 
6.8 x 10~ 2 
1.0 x 10~ 4 
1.4 x 10" 4 
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tt | 
ATLAS [ 
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226 
219 
22C 
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7.5 Conclusion on tqV anomalous couplings 

Table ^ presents a short summary of LHC sensitivities to anomalous FCNC couplings of the top quark. 
For comparison, we present also the estimates of the corresponding sensitivities at Tevatron. For com- 
pleteness we anticipate and include here the results from rare decays discussed in the next section (see 
also [219, 229| |). To unify the description of the LHC potential to detect top anomalous couplings from 
production and decay processes, all results in Table 23 are expressed in terms of limits on top decay 
branching ratios: T(t^qV) /Fsm (= 1-56 GeV). The results were obtained using m t = 175 GeV, 
a s = 0.1, and a = 1/128. When needed the limits quoted in the table have been rescaled to the different 
luminosities and to the S/ V 'S + B > 3 criterion by using a simple linear extrapolation of the available 
bounds (see [30, [23C| ] and Section |8]). The limits on the top anomalous couplings from tV production 
in Table ^l] were obtained using the FC prescription [ 200 ] and have been multiplied by a factor of y/2, 
which roughly relates this prescription with the statistical criterion adopted in Table ^3| [ {226 > 227| ]. 

At present, only few cases (like-sign top-pair production, t—>qZ and t— >qj decays, see [|^, 23C]) 
were investigated with a more or less realistic detector simulation (ATLFAST and CMSJET). Other in- 
vestigations were done at the parton level (the final quarks were considered as jets and a simple smearing 
of lepton, jet and photon energies was applied). Of course, more detailed investigations with a more 
realistic simulation of the detector response may change these results. 

The most promising way to measure the anomalous FCNC top-gluon coupling seems to be the 
investigation of single top production processes, as the search for t—>gq decays would be overwhelmed 
by background from QCD multi-jet events. At the same time, both top quark production and decay would 
provide comparable limits on top quark anomalous FCNC interactions with a photon or a Z-boson. In 
general, the studies shown above indicate that the LHC will improve by a factor of at least 10 the Tevatron 
sensitivity to top quark FCNC couplings. Of course, the results presented here are not complete, since 
other new kinds of interactions may lead to the appearance of unusual properties of the top quark. For 
example, recently proposed theories with large extra-dimensions predict a significant modification of tt 
pair production (see, for example, [231] and references therein). It was found that the exchange of spin-2 
Kaluza-Klein gravitons leads to a modification of the total tt production rate as well as to a noticeable 
deviation in the and M t i distributions with respect to the SM predictions. Naturally, we may expect 
also the modifications of spin-spin correlations due to graviton exchange. 

It has to be stressed that different types of new interactions may affect the same observable quan- 
tity. Only a careful investigation of different aspects of top quark physics may provide a partial separation 
of these interactions. 



8. RARE DECAYS OF THE TOP QUARK0 

The production of 10 — 10 8 top quark pairs per year at LHC will allow to probe the top couplings to 
both known and new particles involved in possible top decay channels different from the main t^bW. 
Thanks to the large top mass, there are several decays that can be considered, even involving the presence 
of on-shell heavy vector bosons or heavy new particles in the final states. On a purely statistical basis, 
one should be able to detect a particular decay channel whenever its branching ratio (BR) is larger than 
about 10~ 6 — 10~ 7 . In practice, we will see that background problems and systematics will lower this 
potential by a few orders of magnitude, the precise reduction being dependent of course on the particular 
signature considered. We will see, that the final detection threshold for each channel will not allow the 
study of many possible final states predicted in the SM, unless new stronger couplings come into play. 



8.1 Standard Model top decays 

In this section, we give an overview of the decay channels of the top quark in the framework of the SM. 
In the SM the decay t—>bW is by far the dominant one. The corresponding width has been discussed in 

23 Section coordinators: B. Mele, J. Dodd (ATLAS), N. Stepanov (CMS). 
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Section 2A . The rates for other decay channels are predicted to be smaller by several orders of magnitude 
in the SM. The second most likely decays are the Cabibbo-Kobajashi-Maskawa (CKM) non-diagonal 
decays t— >sWa.nd t—>dW. Assuming \Vt a \ — 0.04 and \ Vtd\ — 0.01, respectively [192], one gets 



BR(t^sW) ~ 1.6 x 10~ 3 and BR(t-*dW) ~ 1 x 10~ 4 
in the SM with three families. From now on, for a generic decay channel X, we define 

r(t-x) 



BR{t^X) 



T(t^bW) SM ' 



(57) 



(58) 



The two-body tree-level decay channels are the only ones that the LHC could detect in the framework of 
the SM. With the exception of higher-order QED and QCD radiative decays, the next less rare processes 
have rates no larger than 10~ 6 . 




At tree level, the decay t^bWZ (Fig. 43) has some peculiar features, since the process occurs 
near the kinematical threshold (mt ~ Mw + Mz + m&) [232, [233|, [234 W~ 



This fact makes the W 



and Z finite- width effects crucial in the theoretical prediction of the corresponding width [233]. Because 
the W and Z are unstable and not observed directly, more than one definition of the t—>bW Z branching 
ratio is possible. If defined according to 



T(t->bWZ) 



T{t^b^iv^u e u e ) 



BR(W^/iJ/^)BR(Z^i/ e P c ) 



(59) 



including a consistent treatment of W and Z width effects, the branching ratio is to a very good approx- 



imation given by the double resonant set of diagrams (shown in Fig. 43), since the background to the 
neutrino decay of the Z is negligible. One obtains | |235[ ], for mt = 175 GeV, 



BR(t^bWZ) = BR res {t^bWZ) = 2.1 x 10" 



(60) 



However, the signature bpLV^,v e v e is not practical from an experimental point of view. In [ |233[ ], a first 
estimate of BR(t—>bW Z) was given on the basis of the definition 



r(t->bwz) 



T(t^bfih i ^e + e 



BR(W^fj,u^)BR(Z^e+e-) ' 



(61) 
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which involves experimentally well-observable decays, but includes contributions to the numerator from 
t—tbW'y decays (with 7— >e + e~) and other "background" diagrams. The estimate for the corresponding 
branching ratio is 

BR cut (t^bWZ) ~ 6 x 1(T 7 , (62) 

for mt = 175 GeV, assuming a minimum cut of 0.8Mz on the e + e _ -pair invariant mass. This cut tries 
to cope with the contribution of background graphs where the e + e~ pair comes not from a Z boson but 
from a photon. 



If the Higgs boson is light enough, one could also have the decay t^bWH (Fig. 43), although the 



present limits on mn strongly suppress its rate. For ran ^ 100 GeV, one gets [233]: 

BR(t^bWH) < 7 x 10~ 8 . (63) 

in the amplitude. One then 



Finally, the decay t^cWW is very much suppressed by a GIM factor 

M w 



gets [g3J: 

BR{t^cWW) ~ 10~ 13 . (64) 

One can also consider the radiative three-body decays t—>bWg and t—^bW'y. These channels suffer 
from infrared divergences and the evaluation of their rate requires a full detector simulation, including 
for instance the effects of the detector resolution and the jet isolation algorithm. In an idealised situation 
where the rate is computed in the t rest frame with a minimum cut of 10 GeV on the gluon or photon 



energies, one finds [236]: 

BR{t^bWg) 



0.3 , BR(t^bW-f) ~ 3.5 x 10" 



(65) 



The FCNC decays t—>cg, t—*cry and t— >cZ occur at one loop, and are also GIM suppressed by a 



factor jjfi- in the amplitude. Hence, the corresponding rates are very small [ [2490 : 



BR(t^cg) ~ 5 x 10~ n , BR(t^c7) ~ 5 x 10~ 13 , BR(t->cZ) ~ 1.3 x 10~ 13 (66) 
For a light Higgs boson, one can consider also the FCNC decay t—>cH . A previous evaluation of its 



rates [249] has now been corrected. For mjj ~ 100 (160) GeV, one gets [237] 

BR(t^cH) ~ 0.9 x 10~ 13 (4 x 10~ 15 ). 



(67) 



To conclude the discussion of rare SM decays of the top quark, we point out here the existence 
of some studies on semi-exclusive £-quark decays where the interaction of quarks among the t decay 



products may lead to final states with one hadron (meson) recoiling against a jet. In [238] decays with 
an T meson in the final state and decays of the top through an off-shell W with virtual mass Myy* near 
to some resonance M, like 7r + , p + , K + , Df, were considered. An estimate for the latter case is 



T(t^bM) 



144vr 



fll\ V qq'\ 



The typical values of the corresponding branching ratios are too small to be measured: 

BR(t^6vr) ~ 4 • 10~ 8 , BR(t^bD s ) ~ 2 • 10~ 7 . 



(68) 



(69) 



In Table 24 we summarize the expected decay rates for the main top decay channels in the SM. 



8.2 Beyond the Standard Model decays 

The fact that a measurement of the top width is not available and that the branching ratio BR(t—>bW) 
is a model dependent quantity makes the present experimental constraints on the top decays beyond the 
SM quite weak. Hence, the possibility of t decays into new massive states with branching fraction of 
order BR(t— >bW) is not excluded. Apart from the production of new final states with large branching 
fractions, we will see that new physics could also give rise to a considerable increase in the rates of many 
decay channels that in the SM framework are below the threshold of observability at the LHC. 



75 



Table 24: Branching ratios for the main SM top decay channels. 



channel 


BRsm 


channel 


BRsm 


bW 


1 


sW 


1.6 • 10~ 3 


dW 


~ 10~ 4 


bWg 


0.3 (E g > 10 GeV) 


bW-f 


3.5 • 10" 3 (£ 7 > 10 GeV) 


bWZ 


2 • KT S 


cW + W~ 


~ io- 13 


bW+H 


< io- y 


qg 


5 • lO" 11 


<n 


5 • 10" 13 


qZ 


1.3- 10" i3 


cH 


< 10~ 13 



8.21 4 th fermion family 

Extending the SM with a 4 th fermion family can alter considerably a few t decay channels. First of 
all, when adding a 4 th family to the CKM matrix the present constraints on the \ Vt q \ elements are con- 



siderably relaxed. In particular, \Vt s \ and \V t d\ can grow up to about 0.5 and 0.1, respectively Q192| ]. 
Correspondingly, assuming \Vtb\ ~ 1 for the sake of normalisation, one can have up to : 

BR 4 (t^sW) ~ 0.25 and BR 4 (t^dW) ~ 0.01, (70) 



to be confronted with the SM expectations in (57). 



The presence of a 4 th fermion family could also show up in the t direct decay into a heavy b 1 



quark with a relatively small mass (my ~ 100 GeV) [239]. This channel would contribute to the 
t^cWW decay, with a rate: 

BR(t^W + b'(-+W~c)) ~ 10" 3 (10^ 7 ) at my = 100 (300) GeV, (71) 

to be confronted with the SM prediction in 



8.22 Two Higgs Doublet models (2HDM's) 

The possibility that the EW symmetry breaking involves more than one Higgs doublet is well motivated 
theoretically. In particular, three classes of two Higgs doublet models have been examined in connection 
with rare top decays, called model I, II and III. The first two are characterised by an ad hoc discrete 
symmetry which forbids tree-level FCNC [ |240| ], that are strongly constrained in the lightest quark sector. 
In model I and model II, the up-type quarks and down-type quarks couple to the same scalar doublet 
and to two different doublets, respectively (the Higgs sector of the MSSM is an example of model II). 



In model III [241, 242], the above discrete symmetry is dropped and tree-level FCNC are allowed. In 
particular, a tree-level coupling tcH is predicted with a coupling constant < ^/m t m c /v (where v is the 
Higgs vacuum expectation value). 

Since enlarging the Higgs sector automatically implies the presence of charged Higgs bosons in 
the spectrum, one major prediction of these new frameworks is the decay t—>bH + , possibly with rates 
competitive with BR(t^bW) for m H + < 170 GeV. In the MSSM, one expects BR(t^bH + ) ~ 1, both 
at small and large values of tan f3. The interaction Lagrangian describing the H + t 6-vertex in the MSSM 



is [243]: 



C Htb = —J- H + [t (m t cot (3 P L + m b tan [3 P R )b + v (mi t&nl3P R )£] + h.c. , (72) 

v 2 M w 

where Pl,r = 1/2(1 =F 75) are the chiral projector operators. 

At tree level the corresponding decay widths of t—>bH + , H + ^tu, and H + ^tb (or, analogously, 
of H + ^cs) are equal to [ p43| ] 

r(^ + ) = a L i^i 2 Wi,4,4) >< 

647T M^mt \ mf mf I 
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T{H + ^tb) 



where A (a, b, c) = a 2 + b 2 + c 2 



(m 2 + m 2 — mj])(m 2 cot 2 f3 + m 2 tan 2 /3) + Am 2 ml 
^ 11 to 2 tan 2 /3, 



327TM 2 , 



(73) 
(74) 



3 5 2 



\ rrijj m H 



(m 2 H — m 2 t — ml)(m 2 cot 2 (3 + ml tan 2 f3) — 4rn 2 m 2 



(75) 



- 2(ab + ac + be), and m# = mn+. 

mt — rrib, one expects H + —>t + v (favoured for large tan j3) and/or 



Consequently, if mn < 

H + ^cs (favoured for small tan/3) to be the dominant decays. Hence, for tan/3 > 1 and H + ^t + v 
dominant, one can look for the channel t—>bH + by studying a possible excess in the r lepton signature 



from the t pair production [244]. On the other hand, if tan/3 < 2 and m# > 130 GeV, the large mass 
(or coupling) of the t-quark causes ER(H + ^>t*b^W + bb) to exceed BR(H + ^cs) (Fig. |44], see Q245Q 
for details). 




As a consequence, new interesting signatures at LHC such as leptons plus multi-jet channels with 
four 6-tags, coming from the gluon-gluon fusion process gg^tbH^ , followed by the H~^ib decay, 
have been studied [ 246 ]. These processes could provide a viable signature over a limited but interesting 
range of the parameter space. 

One should recall however that both BR(t^bH + ) and BR(H + —>W + bb) are very sensitive to 
higher-order corrections, which are highly model dependent [ ]247j ]. 

In model III, the tree-level FCNC decay t—>ch can occur with branching ratios up to 10~ 2 [ |242| ]. 
In [ |248| ], the rate for the channel t^ch^cWW (cZ Z) has been studied. Accordingly, BR(t^cWW) 
can be enhanced by several orders of magnitude with respect to its SM value. In particular, for an on- 
shell decay with 2A% < m^ < wi[, one can have up to BR(t— >cWW) ~ 10~ 4 from this source. The 
same process was considered in a wider range of models, where the decay t^cWW can occur not only 
through a scalar exchange but also through a fermion or vector exchange [239]. In this framework, the 
fermion exchange too could lead to detectable rates for t^cWW, as in (71 ). 
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In 2HDM's, the prediction for the FCNC decays t— >cg, t—*cy and t— >cZ can also be altered. 
While in models I and II the corresponding branching fractions cannot approach the detectability thresh- 
old | p49|], in model III predicts values up to BR(t^cg) ~ 1(T 5 , BR(t^c7) ~ KT 7 and BR(t^cZ) ~ 
10~^W7[j. 

By further extending the 2HDM's Higgs sector and including Higgs triplets, one can give rise to a 
vertex HWZ at tree level in a consistent way [250]. Accordingly, the t^bWZ decay can be mediated 
by a charged Higgs (coupled with m t ) that can enhance the corresponding branching fraction up to 
BR(i— >bWZ) ~ 10~ 2 . Large enhancements can also be expected in similar models for the channels 
t->sWZ and t^dWZ. 



8.23 Minimal Super symmetric Standard Model (MSSM) 

Supersymmetry could affect the t decays in different ways. (Here, we assume the MSSM frame- 
work [Eg], with (or without, when specified) R parity conservation.) 

First of all, two-body decays into squarks and gauginos, such as i— >tig, t^bixf, i— ►iiX?, could 



have branching ratios of order BR(i— >bW), if allowed by the phase space (see, i.e. [251] for references) 



QCD corrections to the channel t^t\g have been computed in [252] and were found to increase its width 



up to values even larger than T(t—*bW). Three-body t decays in supersymmetric particles were surveyed 



in 02510 . 

The presence of light top and bottom squarks, charginos and neutralinos in the MSSM spectrum 
could also give rise to a CP asymmetry of the order 10~ 3 in the partial widths for the decays t—>bW + 



and t->bW~ [14C, 253]. 



Explicit i?-parity violating interactions [254] could provide new flavour-changing t decays, both 



at tree-level (as in the channels t— >fb and t^-rbxi t p55p ) and at one loop (as in t— >cD Q256| ]), with 
observable rates. For instance, BR(t—>ci>) ~ 10" 4 — 10~ 3 in particularly favourable cases. 

Another sector where supersymmetric particles could produce crucial changes concerns the one- 
loop FCNC decays t-^cg, t—>cy, t^cZ and t—>cH , which in the SM are unobservably small. In the 
MSSM with universal soft breaking the situation is not much affected, while, by relaxing the universality 



with a large flavour mixing between the 2 nd and 3 rd family only, one can reach values such as J257| , [25 8| ] : 
BRmssm (t-c 5 ) ~ HT 6 , BR M ssm(*-c7) ~ KT 8 , BR MS sm(*-cZ) ~ KT 8 , (76) 
which, however, are still not observable. The introduction of baryon number violating couplings in 



broken i?-parity models could on the other hand give large enhancements [ |218| ], and make some of these 
channels observable: 

BR n (t^cg) ~ KT 3 , BR^(t^c7) ~ 10~ 5 , BR fi (t^cZ) ~ 10" 4 . (77) 

A particularly promising channel is the FCNC decay t—>ch in the framework of MSSM, where 
h = h°,H°,A° is any of the supersymmetric neutral Higgs bosons [ ]259| ]. By including the leading 
MSSM contributions to these decays (including gluino-mediated FCNC couplings), one could approach 
the detectability threshold, especially in the case of the light CP-even Higgs boson, for which one can 
get up to: 

BR MSSM {t^ch°) ~ KT 4 . (78) 

8.24 Anomalous couplings 

In the framework of the top anomalous couplings described in Section 7., one can predict large enhance- 
ments in different FCNC top decay channels. While the t-^cg, t— >cj, and t— >cZ processes are analysed 
in section 7., here we concentrate on the possible FCNC contributions to the top decays into two gauge 
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bosons, t—>qVV, where V is either a W or a Z and q = c,u: 

t^qW + Z , t^qW + W~ , t-> 9 ZZ. (79) 

In the SM, the first two decays occur at tree level, while t^qZZ proceeds only through loop contribu- 
tions. We will see that within the present experimental limits on the top anomalous couplings, the rates 
for these processes can be large with respect to the SM prediction, but are still below the detectability 
threshold at the LHC. 

The FCNC contribution to the first channel in (f79|), for the anomalous coupling k z ps 0.3, has a 
rate of the the same order of magnitude as the SM BR(t—>-bWZ) J260| ] : 

BR FC Nc(t^cW Z) ~ 10~ 6 ps BRsm (t—*bWZ) . (80) 



Top anomalous FCNC interactions with both a photon and a Z-boson contribute to the second 
Contrary to the SM case this amplitude has no GIM suppression. As a result, the 



process in eq.(79) 

corresponding branching ratio can have almost the same value as that of the t 



>qWZ decay [ g60J : 



BRfcncO^cW+W - ) ~ 10- ' > BRsm(*->cWW) 



(81) 



For the t—>qZZ decay mode, a coupling k z ~ 0.3 gives a branching ratio much greater than the 
corresponding SM one ( < 10" 13 p6C|]): 

BRfcnc (t^qZZ) ~ 10- 8 » BR SM (t^qZZ) , (82) 

but still too small to be detected at LHC. 

In summary, the observation of any of these decays at LHC would indicate new physics not con- 
nected with the top FCNC interactions (see, for example, [ |248| ]). 



8.3 ATLAS studies of (rare) top quark decays and couplings 

In ATLAS various analyses have been performed on top decays, using the PYTHIA Monte-Carlo inter- 
faced to a fast detector simulation (ATLFAST). In the following, the most relevant results are reported. 



8.31 BR(t— * bX) and measurement of \ Vtb\ 

The SM prediction BR(i^ W + b) ps 1 can be checked by comparing the number of observed (1 or 2) 
6-tags in a ti sample. The first &-tag is used to identify the event as a ti event, and the second fo-tag (if 
seen) is used to determine the fraction of top decays producing a b quark. Within the three-generation 
SM, and assuming unitarity of the CKM matrix, the ratio of double 6-tag to single 6-tag events is given 
by: 

R 2b/lb = BR(i- Wb)/BR(t^ Wq) = \V tb \ 2 /{\Vtb\ 2 + \V ts ? + \V td ?) = \V tb \ 2 (83) 
The CDF collaboration has used the tagging method in leptonic ti events to obtain the result 



R-26/16 = 0.99 ± 0.29 Q261D , which translates to a limit of \V tb \ > 0.76 at the 95% CL assuming three- 
generation unitarity. If this constraint is relaxed, a lower bound of \Vtb\ > 0.048 at the 95% CL is 
obtained, implying only that |T4&| is much larger than either \Vt s \ or \Vu\- 

The LHC will yield a much more precise measurement of R26/1&. For example, ti events in the 
single lepton plus jets mode can be selected by requiring an isolated electron or muon with pr > 20 
GeV, E™ iss > 20 GeV, and at least four jets with p T > 20 GeV. Requiring that at least one of the 
jets be tagged as a 6-jet produces a clean sample of ti events, with S/B = 18.6, with the remaining 
background coming mostly from M^+jet events [0O|]. Assuming a 6-tagging efficiency of 60%, a sample 
of 820 000 single 6-tagged events would be selected for an integrated luminosity of 10 fb _1 . Of these, 
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276 000 would be expected to have a second 6-tag, assuming the SM top quark branching ratios. This 
ATLAS study indicates that the statistical precision achievable would correspond to a relative error of 
dR-2b/ib/R-2b/ib (stat.) = 0.2% for an integrated luminosity of 10 fb _1 . The final uncertainty will be 
dominated by systematic errors due to the uncertainty in the 6-tagging efficiency and fake 6-tag rates, as 
well as correlations affecting the efficiency for 6-tagging two different jets in the same event. Further 
study is needed to estimate the size of these systematic uncertainties. 



8.32 BR(i— > WX) 

The measurement of the ratio of di-lepton to single lepton events in a ti sample can be used to determine 
BR(i— > WX). In this case, the first lepton tags the ti event, and the presence of a second lepton is used 
to determine the fraction of top quark decays producing an isolated lepton, which can be then be related 
to the presence of a TV (or other leptonically decaying states) in the decay. The SM predicts that R21/11 = 
BR(W — > Iv) « 2/9 where £ = (e, y). Deviations from this prediction could be caused by new physics, 
for example, the existence of a charged Higgs boson. The dominant H + decays in such instances are 
usually considered to be H + —> tv or H + —> cs. In either case, the number of isolated electrons and 
muons produced in top decay would be reduced, and R2//1/ would be less than the SM prediction. 

A study performed by ATLAS [[}(]] shows that with an integrated luminosity of 10 fb _1 , a clean 
sample of about 443 000 ti events in the single lepton plus jets mode could be selected by requiring an 
isolated electron or muon with > 20 GeV, E™ isa > 20 GeV, and at least two 6-tagged jets with p? 
> 20 GeV. To determine R2Z/1Z, one then measures how many of these events have a second isolated 
electron or muon, again with px > 20 GeV, and of the opposite sign to the first lepton. Assuming the 
SM, one would expect a selected sample of about 46 000 di-lepton events with these cuts. Given these 
numbers, the statistical precision achievable would correspond to a relative error of 5R2z/iz/R2//iz (stat.) 
= 0.5% for an integrated luminosity of 10 fb _1 . Further study is required to estimate the systematic 
uncertainty on R21/11 due to the lepton identification and fake rates. 



8.33 Radiative Decays: t— > WbZ, t— > WbH 



The 'radiative' top decay t— > WbZ has been suggested D233Q as a sensitive probe of the top quark mass, 
since the measured value of m t is close to the threshold for this decay. For the top mass of (173 ± 5.2) 



GeV [ ]192| ], the SM prediction, based on the Z— >ee signature and a cut m ee > 0.8Mz (see Section |8.1| ), 
is BR cut (i^ WbZ) = (5A±£l) x 10~ 7 [ g33| ]. Thus, within the current uncertainty 5m t « 5 GeV, the 
predicted branching ratio varies by approximately a factor of three. A measurement of BR(t— ► WbZ) 
could therefore provide a strong constraint on the value of mt- Similar arguments have been made for 
the decay t— ► WbH, assuming a relatively light SM Higgs boson. 

ATLAS has studied the experimental sensitivity to the decay t— > WbZ [30, 262], with the Z 
being reconstructed via the leptonic decay Z— >ll {£. = e,y), and the W through the hadronic decay 
W— > jj. The efficiency for exclusively reconstructing t— > WbZ is very low, due to the soft px spectrum 
of the 6-jet in the t— * WbZ decay. Instead, a semi-inclusive technique was used, where a WZ pair 
close to threshold was searched for as evidence of the t— > WbZ decay. Since the t— ► WbZ decay is 
so close to threshold, the resolution on mwz is not significantly degraded with respect to the exclusive 
measurement. The selection of Z—> 11 candidates required an opposite-sign, same-flavor lepton pair, 
each lepton having pr > 30 GeV and \rj\ < 2.5. The clean Z^ 11 signal allows a wide di-lepton mass 
window to be taken (60 GeV < m« < 100 GeV) in order to have very high efficiency. Candidates for 
W— ► jj decay were formed by requiring at least two jets, each having pr > 30 GeV and \rj\ < 2.5, and 
satisfying 70 GeV < mjj < 90 GeV. The 11 jj invariant mass resolution was a[mwz] = 7.2 ± 0.4 GeV, 
and the signal efficiency was 4.3%. 

The dominant backgrounds come from processes with a Z boson in the final state, primarily 
Z+jet production, and to a much lesser extent from WZ and ti production. In order to reduce the 
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Z+jet background, an additional cut requiring a third lepton with px > 30 GeV was made. For the 
signal process tt — ► (WbZ)(Wb), this cut selects events in which the W from the other top decays 
leptonically. After this selection, and with a cut on mwz of ±10 GeV around the top mass, the total 
expected background was reduced to 1.5 events (mostly from WZ production) per 10 fb _1 . Requiring 
at least five events for signal observation leads to a branching ratio sensitivity of order 10~ 3 . Since the 
background has been reduced essentially to zero, the sensitivity should improve approximately linearly 
with integrated luminosity. However, even with a factor of ten improvement for an integrated luminosity 
of 100 fb _1 , the sensitivity would still lie far above the SM expectation of order 10~ 7 — 10~ 6 . 

Given this result, observation of the decay t— > WbH does not look possible. The current LEP 
limit on mu implies that the Higgs is sufficiently heavy that, in the most optimistic scenario that the 
Higgs mass is just above the current limit, BR(i— > WbH) Ss BR(t— ► WbZ). As m# increases further, 
BR(t^ WbH) drops quickly. Assuming mjj ~ mz, one would have to search for t— ► WbH using 
the dominant decay H— > bb. The final state suffers much more from background than in the case of 
t—> WbZ, where the clean Z — ► £ + £~ signature is a key element in suppressing background. Although 
~BR(H—> bb) in this mu range is much larger than BR(Z — > £ + l~ ), the large increase in background 
will more than compensate for the increased signal acceptance, and so one expects the sensitivity to 
BR(i— >• WbH) to be worse than for BR(i— >• WbZ). The decay t— ► WbH has therefore not been studied 
in further detail. 



8.34 i-> H+b 

Limits on the mass of the charged Higgs have been obtained from a number of experiments. An indirect 
limit obtained from world averages of the r branching ratios excludes at 90% CL any charged Higgs with 
mjj+ < 1.5 tan (3 GeV Q263| ] , where tan f3 is the ratio of the vacuum expectation values of the two Higgs 
doublets. CLEO indirectly excludes m H + < 244 GeV for tan > 50 at 95% CL, assuming a two-Higgs- 
doublet extension to the SM [212], while the LEP experiments directly exclude m H + < 59.5 GeV/c 2 at 
95% CL [264]. Searches at the Tevatron have extended the region of excluded [m H ± , tan 0] parameter 
space, particularly at small and large tan (3, and set a limit on the branching ratio BR(t^ H + b) < 
0.45 at 95% CL [§65|]. Run 2 at the Tevatron will be sensitive to branching fractions BR(i-> H + b) > 



11% [266] 



ATLAS has performed an analysis of the experimental sensitivity to the t— > H + b decay, followed 
by H + ^ tv, in the context of the MSSM |f(J gg7| |. Since the relevant t— ► H + b branching ratio is 
proportional to (mf cot 2 (3 + ml tan 2 (3) (see (^3])), for a given value of m H ± the branching ratio for 
such decays is large at small and at large tan [3, but has a pronounced minimum at tan f3 ~ \fm~tfm~b ~ 
7.5. The exact position of this minimum and its depth is sensitive to QCD corrections to the running 
6-quark mass. 

In the ATLAS analysis, an isolated high-p^ lepton with | rj | < 2.5 is required to trigger the ex- 
periment, which in signal events originates from the semi-leptonic decay of the second top quark. One 
identified hadronic tau is then required, and at least three jets with pt > 20 GeV and | r\ \ < 2.5, of which 
two are required to be tagged as 6-jets. This reduces the potentially large backgrounds from l^+jet and 
bb production to a level well below the tt signal itself. These cuts enhance the r-lepton signal from H 
decays with respect to that from W decay, and select mostly single-prong r-decays. After the selection 
cuts and the r identification criteria are applied, t— * H + b decays appear as final states with an excess of 
events with one isolated T-lepton compared to those with an additional isolated electron or muon. 

A signal from charged Higgs-boson production in tt decays would be observed for all values of 
m H ± below m t — 20 GeV over most of the tan f3 range. For moderate values of tan (3, for which the 
expected signal rates are lowest, the accessible values of m H ± are lower than this value by 20 GeV. The 
limit on the sensitivity to BR(t^ H + b) is dominated by systematic uncertainties, arising mainly from 
imperfect knowledge of the r-lepton efficiency and of the number of fake r-leptons present in the final 
sample. These uncertainties are estimated to limit the achievable sensitivity to BR(t^ H + b) = 3%. 
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For charged Higgs masses below 150 GeV and for low values of tan (3, the H — ► cs and if =t — > cb 
decay modes are not negligible. In the same mass range, the three-body off-shell decays if 1 * 1 — > hW*, 
H^ 1 ^ AW* and if^— » bt*— > bbW also have sizeable branching ratios. When the phase-space increases, 
for 150 GeV < m H ± < 180 GeV, both the bbW and the hW* mode could be enhanced with respect 
to the tv mode. Decays into the lightest chargino Xi an d neutralino x.i or decays into sleptons would 
dominate whenever kinematically allowed. For large values of tan [3 the importance of these SUSY 
decay modes would be reduced. However, for values as large as tan (3 = 50, the decay H —* fv would 
be enhanced, provided it is kinematically allowed and would lead to r's in the final state. Their transverse 
momentum spectrum is, however, expected to be softer than that of r's from the direct if ± — > tv decays. 

The if ± — > cs decay mode has been considered as a complementary one to the H^—t tv channel 
by ATLAS for low values of tan f3. In the ATLAS analysis, one isolated high px lepton with | rj \ < 
2.5 is required to trigger the experiment, which in signal events originates from the semi-leptonic decay 
of the second top quark. Two 6-tagged jets with px > 15 GeV and | rj \ < 2.5 are also required, with 
no additional 6-jet. Finally, at least two non-6 central jets with | r\ | < 2.0 are required for the H^— > cs 
reconstruction, and no additional jets above 15 GeV in this central region. Evidence for is searched 
for in the two-jet mass distribution. The mass peak from an decay can be reconstructed with a 
resolution of a = (5 — 8)GeV if the mass of the if ± is in the ranee between 110 and 130 GeV. In this 
mass range, the peak sits on the tail of the reconstructed W— > jj distribution from ti background events 
which decay via a Wb instead of a H ± b. In the mass range 110 < H± < 130 GeV, the ff± peak can 
be separated from the dominant W^> jj background, with S/B m 4-5% and S/\ /r B ~ 5. This channel 
is complementary to the if ± — > tv channel for low tan j3 values. Whereas the if^— » tv channel allows 
only the observation of an excess of events, it is possible to reconstruct a mass peak in the if 1 * 1 — > cs 
decay mode. 

The if±-> hW*, if ^ AW* and fl^-c bt*^ bbW have not been studied so far by ATLAS. With 
the expected ^-tagging efficiency, these multi-jet decay modes are very interesting for a more detailed 
investigation. 

8.35 t — > Zq decay 

The sensitivity to the FCNC decay t — > Zq (with q = u,c) has been analyzed [ ]268| ] by searching for a 
signal in the channel ti — ► (Wb)(Zq), with the boson being reconstructed via the leptonic decay Z — > 11. 
The selection cuts required a pair of isolated, opposite sign, same flavor leptons (electrons or muons), 
each with px > 20 GeV and \rj\ < 2.5 and with \mu — mz\ < 6 GeV. The dominant backgrounds come 
from Z + jet and WZ production. Not only cuts were applied on the Zq final state, but also on the Wb 
decay of the other top quark in the event, to further reduce the background. Two different possible decay 
chains have been considered: the first ("leptonic mode") where the W decays leptonically W — > Iv, 
and the second ("hadronic mode") with W — > jj. The hadronic W decay signature has a much larger 
branching fraction, but suffers from larger backgrounds. The search in the leptonic mode required, in 
addition to the leptons from the Z boson decay, a further lepton with px > 20 GeV and \r]\ < 2.5, E™ lss 
> 30 GeV, and at least two jets with px > 50 GeV and \rj\ < 2.5. Exactly one of the high px jets was 
required to be tagged as a 6-jet. The invariant mass spectrum of each Zq combination was then formed 
from the Z —> 11 candidates taken with each of the non 6-tagged jets. The Zq invariant mass resolution 
was 10.1 GeV. Combinations were accepted if mz q agreed with the known top mass within ± 24 GeV. 
Assuming an integrated luminosity of 100 fb _1 , 6.1 signal events survive the cuts with 7 background 
events. A value of BR(t^Zg) as low as 2 • 10~ 4 could be discovered at the 5a level. 

The search in the hadronic mode required, in addition to the Z —> 11 candidate, at least four jets 
with px > 50 GeV and \r]\ < 2.5. One of the jets was required to be tagged as a 6-jet. To further reduce 
the background, the decay t — > jjb was first reconstructed. A pair of jets, among those not tagged as a 
6-jet, was considered a W candidate if \rrijj — My/\ < 16 GeV. W candidates were then combined with 
the 6-jet, and considered as a top candidate if \mjjt, — m t \ < 8 GeV. For those events with an accepted 
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t — > jjb candidate, the invariant mass of the Z candidate with the remaining unassigned high px jets was 
calculated to look for a signal from t — > Zq decays. Combinations were accepted in case \mz q — m t \ 
< 24 GeV. Assuming an integrated luminosity of 100 fb _1 , one would get 0.4 signal events, with 2 
background events. 



8.36 t — > 7^ decay 

The FCNC decay t — > (with q = u,c) can be searched for as a peak in the Mjj spectrum in the 
region of m t . The requirement of a high px isolated photon candidate in tt — ► {Wb)(jq) events is 
not sufficient to reduce the QCD multi-jet background to a manageable level. Therefore, the t — > Wb 
decay of the other top (anti-) quark in the event was reconstructed using the leptonic W — > iv decay 
mode, and looking for the tt — ► (Wb)(^q) —>■ (lvb)(^q) final state. For the event selection, the ATLAS 
collaboration [ 30|, 262] required the presence of an isolated photon with pt > 40 GeV and < 2.5, an 



isolated electron or muon with px > 20 GeV and \rj\ < 2.5, and E™ ss > 20 GeV. Exactly 2 jets with px 
> 20 GeV were required, in order to reduce the tt background. At least one of the jets was required to 
be tagged as a 6-jet with px > 30 GeV and \rj\ < 2.5. The t — > ivb candidate was first reconstructed. 
The combination was accepted as a top quark candidate if m^ agreed with mt within ±20 GeV. For 
these events the t — > jq decay was sought by combining the isolated photon with an additional hard jet 
with px > 40 GeV and \rj\ < 2.5. The invariant mass of the jj system was required to agree with the 
known value of m t within ±20 GeV. The m 7 j resolution with the cuts described above was 7.7 GeV, and 
the signal efficiency (not counting branching ratios) was 3.3%, including a 6-tagging efficiency of 60%. 
The background (155 events for an integrated luminosity of 100 fb -1 ) is dominated by events with a real 
W — ► iv decay and either a real or a fake photon. These processes include tt, single top production, 
W + jets and Wbb production. The corresponding 5a discovery limit is 

BR(t -> -yq) = 1.0 x 10 -4 . (84) 



8.37 t — > gq decay 

The search for a FCNC tgq coupling (with q 



u, c) through the decay t — > gq was analyzed in [229] 



for the Tevatron. However, as can be seen from Table |23| in Section 7., the sensitivity for such a coupling 
turns our to be much larger in the t production processes than in the decay t — > gq, whose signal will be 
overwhelmed by the QCD background. We refer the reader to Section 7. for a detailed discussion of this 
point. 



8.4 CMS studies of FCNC top quark decays and t H+ b 



The CMS sensitivity to t — * j(Z)(u,c) decays was studied recently (see [23C] for details). The 
PYTHIA 5.7 \$2\ \ generator was used for the signal and background simulations and the detector re- 
sponse was simulated at the fast MC level (CMSJET [ |177[ ]). For the t — ► 7(1*, c) signal the exact 2 — ► 5 
matrix elements gg(qq) — * tt — * ju(c) + W*b{-^ Ivb) were calculated and included in PYTHIA. The 
t — * j(Z)(u, c) decays would be seen as peaks in the Mjfz)jet spectrum in the region of m t . To separate 
the signal from the background one has to exploit the presence of the additional top decaying to the Ivb 
in the same event. The signature with the hadronic decay of the additional top was found to be hopeless. 



8.41 t -> 7(n, c) 

In order to separate the (75) (ivb) final state from the backgrounds several selection criteria were 
found to be effective. First, the presence of one isolated photon with E t > 75 GeV and |r/| > 2.5, one 
isolated lepton (/i, e) with E t > 15 GeV and |?7| > 2.5, and at least two jets with E t > 30 GeV and 
\r]\ > 2.4 is required. One top quark has to be reconstructed from the photon and jet (Mjj et C mt ± 15 
GeV), the corresponding jet is not allowed to be 6-tagged. On the contrary, the jet with maximal E t , 
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which is not involved in the (7, jet) system has to be 6-tagged, should have E t > 50 GeV and contribute 
to another reconstructed top quark (M C m t ± 25GeV). There must be no additional jets with E t > 50 
GeV. The 6-tagging efficiency was assumed to be 60% for the purity 1%(10%) with respect to the gluon 
and light quark jets (c-quark jets). After this selection, approximately 270 background events dominated 
by the ti and W + jets, including Wbb, survive for the integrated luminosity of 100 fb , while the signal 
efficiency is 9.1%. The S/B ratio is about 1 for BR(i — > j(u, c)) = 10 -4 and the 5 a discovery limit is 
as low as 3.4 x 10~ 5 for 100 fb -1 . 



8.42 t -> Zq 

The t — > Zq signal was searched for in the ti — > {££q)£ub final state. Three isolated leptons with 
Et > 15 GeV and \rj\ < 2.5, and exactly two jets with E t > 30 GeV and \q\ < 2.5 are required. The 
pair of the opposite-sign same-flavour leptons has to be constrained to the Z mass (Mgi C Mz ± 6 GeV) 
and one jet, combined with the reconstructed Z, has to form the top system (M^- C m t ± 15 GeV). 
This jet is not allowed to be the 6-jet, but the last "free" jet in the event has to be 6-tagged. For the 
integrated luminosity of 100 fb -1 just ~ 9 background events coming from the WZ, tiZ and Z + jets 
processes survive. The signal efficiency is about 6.8% which corresponds, however, only to ~ 12 events 
for BR(i — ► Z[u, c)) = 10~ 4 . The indication is that one can reduce the background rate to the nearly 
zero level tightening the selection criteria. In particular, requiring in addition E™ %ss > 30 GeV and a 
harder jet involved in the top (££j) system (E t > 50 GeV) one can reduce the background to the level 
of ~ 0.6 events still keeping ~ 3.7% of the signal (6.6 events for BR(t -> Z(u,c)) = 10~ 4 and 100 
fb -1 ). One can conclude that the t — > Z (u, c) signal should be very clean but, due to the low signal event 
rate, only ~ 3 x 100fb _1 of integrated luminosity would allow one to probe BR(t — > Z(u, c)) as low as 
10~ 4 , provided the present background understanding is correct and the detector performance will not 
be deteriorated during the long run. The 5a reach for 100 fb -1 is ~ 1.9 x 10 -4 . 



8.43 t-> H+b 

CMS has investigated the production of the light charged Higgs, m H ± < m t ,m it events using the decay 



chain fii-+H ± bWb-+(T ± v T b) + (£ub) [ Et»|| . The -^tv branching ratio is large ~98% in this mass 
range for tan (3 >2 and only slightly dependent on tan /?. The t—tH^b branching ratio is large both at 
high and at low tan/3 values and has a minimum of ~0.8% around tan/3 ~6. Since the Higgs mass 
cannot be reconstructed in this process the signal can be only inferred from the excess of r production 
over what is expected from the SM t—>Wb, W — ► decay. 

An isolated lepton with p t > 20 GeV is required to identify the top decay and to trigger the 
event. The r's are searched starting from calorimeter jets with E t >40 GeV within \rj\ <2.4. For the 
r identification the tracker information is used, requiring one hard isolated charged hadron with p t >30 
GeV within the cone of AR <0.1 inside the calorimeter jet. The algorithm thus selects the one prong r 
decays. 

The main backgrounds are due to the it events with ti^>WbWb^(r v T b) + (£ub) and W + 
jet events with W—>tv. The ti background is irreducible, but can be suppressed by exploiting the 



r polarisation effects [27C]. Due to the r polarisation the charged pion from T—t^v decay has a 



harder p t spectrum when coming from H^^tu than from W^^tv. The decay matrix elements with 



polarisation [271] were implemented in PYTHIA [52]. Due to the polarisation, the efficiency of the above 



r selection is significantly better for H^—tTv (~19%) than for W^^tv (~6%). 

The events were required to have at least one 6-jet with E t > 30 GeV tagged with an impact 



parameter method [272]. This 6-tagging suppresses efficiently, by a factor of ~70, the background from 
VF+jet events. The efficiency for ti events is ~35%. The expected 5c discovery range for 10 fb -1 
in the MSSM (uia, tan /3) parameter space was found to be: uia < 110 GeV for all tan /3 values and 
somewhat extended (rriA ~ 140) for tan [3 ^ 2. 
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8.5 Conclusions on rare top decays 

In the framework of the SM, the top rare decays (that is any channel different from t-^qW) are definitely 
below the threshold for an experimental analysis at LHC. On the other hand, LHC experiments will be 
able to probe quite a few predictions of possible extensions of the SM. 

An extended Higgs sector will be looked for through the tree-level decay t^bH + . ATLAS esti- 
mates its sensitivity to this channel in the MSSM, through an excess in the tau lepton signal, to be around 
BR(t—>H + b) = 3% (that is almost 4 times better than what expected from Run 2 at the Tevatron). This 
would allow to probe all values of m H ± below mt — 20 GeV over most of the tan (3 range. For low 
tan/3, the complementary decay mode H^^-cs has been considered. In the mass range 110 < < 
130 GeV, the if* peak can be reconstructed and separated from the dominant W^jj background. 

For CMS, using the r excess signature, the expected 5a discovery range for 10 ftr 1 in the MSSM 
(m a, tan (3) parameter space is uia < 110 GeV, for all tan/3 values, and somewhat extended (wia ^ 
140), for tan < 2. 

Other interesting signatures like H^ — >hW*, H^—i-AW* and H —*bt*—*bbW are very promising 
in particular parameter ranges, but have not yet been thoroughly investigated. 

ATLAS has studied its sensitivity to the radiative decay t—>WbZ. This has been found to be at 
most of the order 10~ 4 , that is insufficient for the study of a SM signal (~ 10~ 6 ), but possibly useful for 
exploring the predictions of some extended Higgs-sector model, for which BR(t— >WqZ) 10~ 2 . On 
the other hand, the radiative Higgs decay t—>WbH seems out of the reach of LHC in any realistic model. 

The LHC reach for the FCNC decays t—>qZ, t—^qj and t—>qg has also been thoroughly in- 
vestigated. Apart from the t—>qg, which is completely overwhelmed by the hadronic background, both 
ATLAS and CMS have a sensitivity of about 2 x to the t^qZ channel, while the CMS reach for 
the t^qj channel is about 3.4 x 10 -5 , that is slightly better than the ATLAS sensitivity (1.0 x 10~ 4 ), 
assuming an integrated luminosity of 100 fb~ 4 . These thresholds could be largely sufficient to detect 
some manifestation of possible FCNC anomalous couplings in the top sector. 

ATLAS has also investigated its sensitivity to a measurement of \Vtb\ through a determination 
of the rate BR(t—>bX), by comparing the number of observed (1 or 2) fe-tags in a ti sample. Within 
the three-generation SM, the ratio of double 6-tag to single 6-tag events is R,2b/ib = \Vtb\ 2 - LHC will 
allow a much more precise determination of R2&/1& with respect to the Tevatron (where, presently, one 
gets \Vtb\ > 0.76 at the 95% CL). On a purely statistical basis, the expected relative error on R-26/ifo ^ 
^F^b/ib/R^ft/ifc (stat.) = 0.2% for an integrated luminosity of 10 fb _1 , that would imply a relative error 
on \ Vtb\ of about l%o . On the other hand, the final uncertainty will be dominated by systematic errors 
related to the 6-tagging. Further study is needed to estimate the size of these systematic uncertainties. 



9. ASSOCIATED TOP PRODUCTION^ 

The associated production of a Higgs boson (both SM-like and MSSM) with a top-antitop pair, is one of 
the most promising reactions to study both top quark and Higgs boson physics at the LHC. 

The pp — > tiH channel can be used in the difficult search for an intermediate mass Higgs (mjj ~ 
100— 130 GeV), as first proposed in [273]. In this mass region, the associated top production cross 
section is quite high but still smaller than the leading gg — > H and qq — > Hqq cross sections by two 
orders and one order of magnitude, respectively. However, since the final state tiH signature is extremely 
distinctive, even such a small signal production rate can become relevant, especially if identifying the 
Higgs through its dominant H — > bb decay becomes realistic, as will be discussed in the following. 

Associated tiH production will furthermore provide the first direct determination of the top quark 
Yukawa coupling, allowing to discriminate, for instance, a SM-like Higgs from a more general MSSM 
Higgs. Processes like gg — > H or H — > 77 are also sensitive to the top Yukawa coupling, but only 

24 Section coordinators: A. Belyaev, L. Reina, M. Sapinski (ATLAS), V. Drollinger (CMS). 
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Fig. 45: Cross section for ttH production at the LHC as a 
function of the Higgs mass, for /i = rriH- 
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Fig. 46: Diagrams for gg—>ttH, the leading parton 
level process for pp — * ttH. 



through large top loop corrections. Therefore loop contributions from other sources of new physics can 
pollute the interpretation of the signal as a measurement of the top Yukawa coupling. 

In the following we will concentrate on the case of a SM-like Higgs boson, whose top Yukawa 
coupling (y t = 2 3 / 4 G]/ 2 m t ) is enhanced with respect to the corresponding MSSM (scalar Higgs) cou- 
pling for tan > 2, the region allowed by LEP data. Predictions for the MSSM case can be easily 
obtained by rescaling both the ttH coupling and any other coupling that appears in the decay of the 
Higgs boson. 

The cross section for pp — ► ttH at LO in QCD has been known for a long time [ 274 ] and has 
been confirmed independently by many authors. We have recalculated it and found agreement with the 
literature. Of the two parton level processes (qq — » ttH and gg — > ttH), gg — ► ttH dominates at 
the LHC due to the enhanced gluon structure function. The complete gauge invariant set of Feynman 
diagrams for gg — > ttH is presented in Fig. 46. The corresponding analytical results are too involved to 
be presented here. The numerical results for yfs = 14 TeV and a few values of the QCD scale \i are given 

-vriH- For consistency, we have used the 



in Table 25, and illustrated in Fig. 
leading order CTEQ4L PDFs [fiT: 



<°0 = M z ) = 0.1317 for A, 



J) 

QCD 



3 as functions of m#, for fj, ■- 

as well as the leading order strong coupling constant (for reference, 
0.181). The cross section, as expected from a LO calculation, 
shows a strong scale dependence, as can be see in Table 25 , where results for \i = mn , nit, m jj + 2mt 
and y/§ are presented. In comparison with fj, = 2mt + mg, for fi = ma we have 80-50% higher cross 
sections, when 100 GeV < m# < 200 GeV. Since the choice of the QCD scale at LO is pretty arbitrary, 
and since we expect NLO QCD corrections to enhance the LO cross section, we decide to use \x = uih in 



Fig. |45| and in the following presentation. These calculations have been performed independently using 
the CompHEP software package p75ll and MADGRAPH ||7|]+HELAS $27% . 

The NLO QCD corrections are expected to enhance the cross section, but their complete evaluation 
is still missing at the moment. Associated top production is in fact the only Higgs production mode for 
which the exact NLO QCD corrections have not been calculated yet. The task is very demanding, since 
it requires the evaluation of several one loop five-point functions for the virtual corrections and the 
integration over a four-particle final state (three of which massive) for the real corrections. 

For large m#, the cross section for ttH has been calculated including a complete resummation 



of potentially large logarithms, of order ln(m.ff/mt), to all orders in the strong coupling [ [27 8Q . These 
effects can almost double the cross section for mjj = 1 TeV. 
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Table 25: Leading order cross sections for ttH production at the LHC. The individual parton level channels {qq — * ttH and 



gg — * ttH) as well as their sum are given for a few values of the renormalization scale fi. 





qq [fb] 


99 [fb] 


qq+gg [fb] 


qq+gg [fb] 


qq+gg [fb] 


qq+gg [fb] 




pi = m H 


fi = m H 


fi = m H 


(j, = m t 


/j, = 2m t + m H 


fi = y/s 


100 


348. 


990. 


1340. 


1070. 


765. 


685. 


110 


279. 


740. 


1020. 


840. 


596. 


534. 


120 


227. 


558. 


785. 


674. 


473. 


422. 


130 


186. 


428. 


613. 


542. 


379. 


338. 


140 


153. 


334. 


487. 


445. 


308. 


273. 


150 


128. 


263. 


391. 


367. 


251. 


224. 


160 


107. 


210. 


317. 


306. 


207. 


184. 


170 


90.5 


169. 


260. 


257. 


173. 


152. 


180 


76.8 


139. 


216. 


218. 


145. 


128. 


190 


65.7 


115. 


181. 


187. 


124. 


108. 


200 


56.4 


97.1 


153. 


162. 


106. 


92.4 


300 


15.0 


29.5 


44.5 


55.7 


33.2 


28.4 


400 


5.11 


15.6 


20.7 


29.6 


16.2 


13.8 


500 


2.04 


9.51 


11.5 


18.4 


9.32 


7.98 


600 


0.909 


6.00 


6.91 


12.1 


5.73 


4.93 


700 


0.439 


3.86 


4.29 


8.20 


3.63 


3.14 


800 


0.226 


2.50 


2.72 


5.62 


2.34 


2.04 


900 


0.122 


1.65 


1.76 


3.90 


1.54 


1.35 


1000 


0.0684 


1.10 


1.16 


2.73 


1.02 


0.900 



For an intermediate mass Higgs, the K factor ((Tnlo/plo) has been estimated in the Effective 
Higgs Approximation (EHA) [ |279j ]. The EHA neglects terms of 0(rriH / \fs) and higher and works 
extremely well for e + e~ — ► ttH already at yfs = 1 TeV. However, it is a much poorer approximation in 
the pp — > ttH case, since it does not include the t -channel emission of a Higgs boson for gg — > ttH. 
Indicatively, at y/s = 14 TeV, for a SM-like Higgs boson with m# ~ 100 — 130 GeV, the EHA gives 
K ~ 1.2 — 1.5, with some uncertainty due to scale and PDF dependence. Only the complete knowledge 
of the NLO level of QCD corrections will allow to reduce the strong scale and PDF dependence of the 
LO and EHA cross sections. For the following analysis we choose to use the pure LO cross section 
with no K-factor, both due to the uncertainty of the result and for consistency with the corresponding 
background cross sections. However, one should point out that, due to the choice of a quite low QCD 
scale (/i = mu), a sort of effective K-factor has been automatically included in our analysis. 

In the following subsection we present the analysis and results from the ATLAS collaboration as 
well as a discussion of the main backgrounds. The analysis mainly focuses on the search and study 
of an intermediate mass Higgs boson. To introduce the study, it is useful to discuss and qualitatively 
understand the size of the possible irreducible backgrounds in the 100 < mu < 140 GeV mass region. 

Given the relatively small number of events that will be available, one should try to consider all 
possible decay channels of the Higgs boson in the intermediate mass region: H — > bb , tt , 77 , WW 
and ZZ. The corresponding irreducible backgrounds are: 1) ttbb, 2) tirf, 3) ££77, 4) ttWW, and 5) 
tiZZ. The number of events expected from signal and background signatures for l)-5) are presented in 
Fig. 47. This figure shows the number of signal and background events in each bin of the corresponding 
invariant mass: M^, M TT , M 77 , Myyw or Mzz- They are obtained multiplying the ttH cross section 
by the respective Higgs boson branching ratios. In order to take into account finite mass resolution 
effects, we have chosen 10 GeV bins for the M 77 distribution and 50 GeV for the others. The presented 
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Table 26: Leading-order cross sections for various ttXX backgrounds. 





ttbb 


ttrf 


tt77 


ttww 


ttzz 


cuts 


\Vb\ < 3 
E b T > 15 GeV 
m bb > 90 GeV 


\Vt\ < 3 
E T T > 15 GeV 


Kl < 2.5 
W T > 15 GeV 
M 77 > 90 GeV 






a [fb] qq 
99 
QQ+99 


41.2 
846. 
887. 


2.9 
15.7 
18.6 


2.73 
1.82 
4.55 


0.50 
1.52 
2.53 


1.11 
0.567 
1.68 



numbers correspond to 30 fb 1 of integrated luminosity. The corresponding total cross sections are given 



in Table 26 



Cross sections for backgrounds l)-3) were calculated with the kinematic cuts shown in Table [26 
while for processes 4) and 5) no cuts were applied. We have used CTEQ4L PDF and fi 2 = Mxx, where 
XX is bb, t + t~, 77, W + W~ or ZZ depending on the channel. One can see that the ttbb signature 
has the highest signal (and background) event rate. It has been the object of the study of the ATLAS 
collaboration and will be discussed in the next section. The W77 channel has also been the subject 
of [ |280| ] where signal as well as reducible and irreducible backgrounds have been studied in details at the 
parton level. However, one can see that other signatures could also be interesting and helpful in searching 
for the Higgs boson and measuring the ttH Yukawa coupling, and should be taken into account in future 
studies. 



9.1 ttH : Analyses and Results 

The ATLAS collaboration has studied several channels in which the discovery of a SM-like Higgs boson 
would be possible and obtained a quite complete Higgs discovery potential [|30|]. One of the most impor- 
tant channels for discovery of a low mass Higgs boson (100 — 130 GeV) is the ttH, H —> bb channel, 
in which it is possible to obtain quite large signal significance [ 281 ] and also to measure the top-Higgs 
Yukawa coupling. 

The final state of this channel consists of two W bosons and four b— jets: two from the decay of 
the top quarks, and two from the decay of the Higgs boson. In order to trigger signal events, one W 
boson is required to decay leptonically. The second W boson is reconstructed from the decay to a q'q 
pair. This channel could be also investigated with both W bosons decaying leptonically. However, for 
this signature the total branching ratio is much smaller and, in addition, it is more difficult to reconstruct 
two neutrino momenta from the measured missing energy. 

In the analysis both top quarks are fully reconstructed, and this reduces most of the VF+jets back- 
ground. The reconstruction is done using strategies similar to those discussed in Section ^ for the 
kinematic studies of ti production. The main backgrounds for this process are: 

• the irreducible continuum ttbb background; 

• the irreducible resonant ttZ background, which is not very important for this channel as it has a 
very small cross section; 

• the reducible backgrounds which contain jets misidentified as 6-jets, such as ttjj, Wjjjjjj, 
WWbbjj, etc. 

After the reconstruction of the two top quarks, it has been found that the most dangerous back- 
ground is ttbb (56% of all tt+jets background). In Table 27 we give ax BR, where BR represents the 
product of the branching ratios for t — > Wb,Wi — > £v,W2 — ► qiq~2, and H — > bb. We also give 
the number of events expected after the reconstruction procedure for 3 years of low luminosity oper- 
ation (6-tagging efficiency e b = 60%; probability to mistag c-jet as 6-jet e c = 10%; probability to 



88 



100 



10' 



a> 
> 



o 10 



CD 



10 




100 



200 



300 



400 500 
M„[GeV 1 



100 



200 



300 



400 



500 100 



200 



300 



"V GeV ] 



400 
M ww [ GeV 



500 



c 
CD 

> 1 
<" 10 



CD p 

-Q 10 2 

Z3 



10 




100 



200 



300 



400 
M zz [ GeV 



500 



Fig. 47: Number of events for tiH signal (solid line) and background tibb, ttyy, ttWW, tiZZ signatures (histogram), as a 
function of the corresponding invariant masses Mxx, assuming 30 fb _1 of integrated luminosity at \/s= 14 TeV. The bin size 
is 10 GeV for the Af 77 distribution and 50 GeV for the others. 

mistag any other jet as tb-jet ej = 1%; pj? > 15 GeV; lepton identification efficiency q = 90%; 
p e j? > 20 GeV), and after one year of high luminosity operation (for high luminosity the 6-tagging effi- 
ciency is degraded to = 50% (e c , ej and eg remain unchanged), the threshold on jet reconstruction is 
raised to pt > 30 GeV and the electron pt threshold is raised to pj. > 30 GeV). Combined results are 
also shown. 

Figure [18| shows the signal and background shapes for uih = 120 GeV and 100 fb _1 of integrated 
luminosity obtained with combined detector performance (30 fb _1 with low luminosity and 70 fb _1 with 
high luminosity). On the other hand, Fig. 49 illustrates the signal shape for mjj = 100 GeV, as obtained 
by using the full (GEANT) simulation of the detector. In this figure, the shaded area represents the true 
signal where both fe-jets come from the Higgs boson, and the solid line stands for the signal obtained 
through the method that we described above. The combinatorial background, which comes from taking 
at least one 6-jet from a top instead the one from the Higgs, is quite large and the signal purity is at the 
level of 60% for low luminosity. 

For the fast simulation the m b i peak mass resolution is c mbi = 19.0 GeV, while for the full 
simulation, including the influence of electronic noise and the threshold on cell energy, a resolution 

0",, 



: 20.0 GeV has been obtained. 
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Table 27: Cross sections multiplied by branching ratios and numbers of events after all cuts, including the ±30 m b i mass win- 
dow cut, for 30 fb _1 (low luminosity detector performance), 100 fb _1 (high luminosity detector performance) and combined 
100 fb _1 (30 fb _1 with low luminosity and 70 fb _1 with high luminosity detector performance) of integrated luminosity. 





U s\ L> 1\ 




nr. of 






(r>b) 


reconstructed events 






low lumi 


high lumi 


combined 


ttH, m H = 120 GeV 


0.16 


40 


62 


83 


tt + jets 


87 


120 


242 


289 


Wjjjjjj 


65200 


5 


10 


12 


ttz 


0.02 


2 


5 


6 


total background 




127 


257 


307 


S/B 




0.32 


0.24 


0.27 


S/V(B) 




3.6 


3.9 


4.8 


^H^bb/ $ total 




59% 


50% 




Syt/yt (stat.) 




16.2% 


14.4% 


11.9% 



Similar analyses have been performed for the ttH, H — > 77 channel. Since the signal rate for this 
channel is very small, it will not be useful during the low luminosity period. However, thanks to the high 
purity of the signal, it will be possible to obtain between 4 or 5 signal events above 1 event from ££77 
background per one year of high luminosity operation [282]. To increase the signal rate, WH and ZH 
with H —>■ 77 channels have been included into the analysis and 14 signal events above 5 background 
events (VK77, Z77, ££77 and 6677) are expected for one year of high luminosity operation [30]. 

The statistical uncertainty in the determination of the top-Higgs Yukawa coupling y t is given in 
the last row of Table 27 These results assume that the theoretical uncertainty is small, as we expect to 
be the case by the time the LHC turns on. Many statistical uncertainties of the direct measurement of 
yt, such as those associated with uncertainties in the integrated luminosity and in the ti reconstruction 
efficiency, could be controlled by comparing the ttH rate with the ti rate. 

To conclude, the ttH, H — > bb and H — > 77 channels are very useful for Higgs boson discovery 
as well as for the measurement the of top-Higgs Yukawa coupling. 



9.11 A closer look at the ttbb background: CompHEP versus PYTHIA 

It is necessary to stress that the correct understanding of the tibb background is one of the main points 
of this study. One can simulate this background using PYTHIA, by generating events of top pair produc- 
tion and emitting bb pairs from the gluon splitting after the initial and final state radiation. In order to 
understand how good or bad this approximation is, one needs to calculate and simulate the complete tibb 
process. We have done this using the CompHEP package [275], 

In order to compare CompHEP and PYTHIA on the same footing, one should take into account 
the effects of the initial and final state radiation in CompHEP. This has been done through a CompHEP- 
PYTHIA interface [283]. We use parton level events generated by CompHEP and link them to PYTHIA in 
order to include initial and final state radiation effects as well as hadronization effects. 

Table ^8| presents parton level CompHEP and PYTHIA cross sections including branching ratios 
of the M^-boson decay, for the same choice of structure function (CTEQ4L Ql 1 5p ) and QCD scale (/i 2 = 
m t +Pr(average)). We can see a good agreement for the total cross sections between the exact calculation 
and the gluon splitting approximation. 

In Fig. ^ we present the distribution of 6-jet separation in tibb events. One can see a quite 
good agreement between CompHEP and PYTHIA. Figures |T| and 52 compare the transverse momentum 
distributions of the most energetic 6-jet and of the least energetic 6-jet in tibb production, as reproduced 
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Fig. 48: Invariant mass distribution of tagged 6-jet 
pairs in fully reconstructed tiH signal events and back- 
ground events, obtained using the fast simulation of the 
ATLAS detector, for ran = 120 GeV and integrated lu- 
minosity of 100 ftT 1 (30 ftT 1 at low plus 70 ftT 1 at 
high luminosity). The points with error bars represent 
the result of a single experiment and the dashed line rep- 
resents the background distribution. 
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Fig. 49: Invariant mass distribution of tagged 6-jet pairs in 
fully reconstructed tiH signal events, obtained using a full 
(GEANT) simulation of the ATLAS detector, for uih = 
100 GeV and low luminosity performance. The shaded area 
denotes those events for which the jet assignment in the Higgs 
boson reconstruction is correct. 



Table 28: Results for the ttbb background, assuming an integrated luminosity d n t = 30 fb 1 : CompHEP (ISR and FSR 



included) versus PYTHIA (default). 



Selection 


CompHEP 


PYTHIA 


CompHEP / PYTHIA 


4 6-quarks with 

p T > 15 GeV/c ; \rj\ < 2.5 


92000 events 
a = 3.1 pb 


87600 events 
a = 2.9 pb 


1.05 


AR(b,b) > 0.5 

6-quarks not from top decay 


54000 events 
59% of prev. Step 


48900 events 
56% of prev. Step 


1.10 



using PYTHIA and CompHEP respectively. These distributions also confirm that PYTHIA describes well 
the ttbb background. 

9.2 Summary and conclusions for ttH production 

The associated production of a Higgs boson with a top-antitop pair is important for the discovery of an 
intermediate mass Higgs boson (mjj ~ 100— 130 GeV) and provides a direct determination of the top- 



Higgs Yukawa coupling. From studies of the couplings and of the CP-parity of the Higgs boson [284] it 
will be possible to discriminate, for instance, a SM-like Higgs boson from a generic MSSM one. 

The ATLAS analysis has focused on the tiH , H —> bb channel for the low luminosity run of the 



LHC (30 fb -1 of integrated luminosity). The results presented in Section 9T are very encouraging and 
indicate that a signal significance of 3.6 as well as a precision of 16% in the determination of the Yukawa 
coupling can be reached (for mjj = 120 GeV). Better results can be obtained from the high luminosity 
run of the LHC (100 fb -1 of integrated luminosity), when also the high purity ttH, H — > 77 channel is 
available. 
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Fig. 50: Results for tibb background, assuming an integrated luminosity Li 
versus PYTHIA (default). 
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Fig. 5 1 : Transverse momentum distribution of the most 
energetic 6-jet from tibb processes at the LHC: CompHEP 
(dark-grey histogram) versus PYTHIA (light-grey) 



Fig. 52: Transverse momentum distribution of the least 
energetic 6-jet from tibb processes at the LHC: Com- 
pHEP (dark-grey histogram) versus PYTHIA (light- 
grey) 

A APPENDIX: b-TAGGING AND JET E-SCALE CALIBRATION IN TOP EVENTS0 

For the reconstruction of the top events and in particular for the precision measurement of the top mass 
two important aspects in the detector performance have to be considered: 

• the 6-quark jet tagging capabilities and efficiency in top events, and 

• the jet energy scale calibration for the light quark jets but in particular for the 6-jets. 

In both experiments ATLAS and CMS several studies have been made on these, highlights of which are 
presented here. From the preliminary results available so far, there is confidence that the numbers used 
or implied in the analyses presented in this report are realistic. Needless to say that these are preliminary 
results and several detailed studies need to be performed with the final detector simulations and the first 
LHC data. 



25 Section coordinator: I. Efthymiopoulos 
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Fig. 53: Left: jet rejection factors for the vertex 6-tagging method, with high luminosity pile-up. Open symbols: mjj = 
100 GeV, full symbols: vtih = 400 GeV. Right: background rejections as a function of the p T for eb = 50% with high 
luminosity pile-up included. From [boj] , Chapter 10. 



Al b-jet tagging in the top events 

ATLAS has done extensive studies for the 6-tagging performance using jets from the decay of 100 and 



400 GeV Higgs bosons ([pO|], Chapter 10). In Fig. |53| the rejection factors for the light quark jets versus 
the b tagging efficiency and the jet p T are shown. 

Typically in the ATLAS analyses discussed here, and in particular for the fast simulation studies, 
an overall 6-jet tagging efficiency of 60% (50%) for low (high) luminosity of LHC is used. The mis- 
tagging inefficiencies for the c-jets (or other light quark jets) were 10% (1%) for the p T range interesting 
for the top physics. Although most of the studies were done with events from the Higgs decays, the 
results were verified with the top events themselves and no significant differences were found. 

A2 Absolute jet energy scale calibration 

Determining the absolute jet energy scale at LHC will be a rather complex issue because it is subject 
to both physics (initial-final state radiation, fragmentation, underlying event, jet algorithm etc.) and 
detector (calorimeter response over a wide range of energies and over the full acceptance of the detector, 
non-linearities at high energies, e/h ratio etc.) effects. All these have to be understood at the level of a 
fraction of a percent in terms of systematic uncertainties as required for the precision measurements of 
the top mass. 

ATLAS has done an extensive study of the possible in situ jet scale calibration methods using 



specific data samples available at LHC ([p0[], Chapter 12). In general, good candidate event classes at 
LHC will be: 



• reconstruction of W — > jj decays within the top events themselves Q12| ] to obtain the light quark 
jet calibration and, 

• events containing a Z boson decaying into leptons balanced with one high-p T jet to cross-check 
the light quark jet calibration but in addition to calibrate the 6-jets and extend the energy reach to 
the TeV range. 



In Fig. 54 the results obtained are shown. As can be seen (left plot) for the case of W — ► jj events, 
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Fig. 54: Left: Ratio of the original parton p T to the p T of the reconstructed jet as a function of the p T of the jet for the 
W — > jj decays reconstructed in include ti events. The jets were reconstruced using a fixed cone jet algorithm with cone size 
DR= 0.4, (optimized for high luminosity operation of LHC). Right: Average fractional imbalance between the p T of the Z 
boson and the p T of the leading jet as a function of the p T of the jet for the sample of Z + jets events. A cone of DR= 0.7 is 
used to collect the jet energy. 



once the jet 4-vectors are reseated using the Mw constraint the required 1% uncertainty is reached for 
jets with p T > 70 GeV up to several hundred GeV. The lower and upper end of this range will depend on 
how well residual systematic effects can be controlled in the data and the Monte Carlo simulation [ J285| ]. 

The use of the Z + jets sample in LHC is a bit less straightforward than at the Tevatron [ 286 ] 
due to the ISR radiation which produces an additional high-p T jet which degrades the quality of the 
p T -balance between the Z boson and the leading jet. In Fig. 54 (right) the variation of the average 
fractional imbalance between the pt of the leading jet and the Z boson as a function of the pt of the jet. 
Rescaling the jet pr to satisfy pr balance with the Z boson and applying tight selection criteria (jet veto 
and difference in azimuth 8(f) between the reconstructed Z and the leading jet) the desired goal of ±1% 
systematic uncertainty on the absolute jet energy scale can be achieved for jets with px > 50 GeV and 
up to the TeV range [287]. 



However, as shown in Fig. q4| (right), it is possible, taking advantage of the large rate and requiring 
tight event selection criteria, to obtain the required precision for jets with p T > 40 GeV and up to the TeV 
range. 

Clearly more studies are needed, and will be done in the years to come, to understand the limita- 
tions of the proposed methods and to devise possible improvements. 



B APPENDIX: DIRECT MEASUREMENT OF TOP QUANTUM NUMBERS^ 
Bl Top spin and experimental tests 

Evidence to date is circumstantial that the top events analysed in Tevatron experiments are attributable 
to a spin- 1/2 parent. The evidence comes primarily from consistency of the distribution in momentum 
of the decay products with the pattern expected for the weak decay t — > b + W, with W —> £ + v or 
W — ► jets, where the top t is assumed to have spin- 1/2. 

26 Section coordinators: E.L. Berger, U. Baur. 
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Fig. 55: (a) Partonic cross sections a(M) as functions of partonic sub-energy M for the gg channel. 
da/dM in proton-proton collisions at 14 TeV as functions of pair mass. The top quark mass mt = 
(stop) mass m t - = 165 GeV. 
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It is valuable to ask whether more definitive evidence for spin- 1/2 might be obtained in future 
experiments at the Tevatron and LHC. We take one look at this question by studying the differential 
cross section da/dM t i in the region near production threshold [ |288| ]. Here M t j is the invariant mass of 
the tt pair. We contrast the behaviour of tt production with that expected for production of a pair of 
spin-0 objects. We are motivated by the fact that in electron-positron annihilation, e + + e~ — > q + q, 
there is a dramatic difference in energy dependence of the cross section in the near-threshold region for 
quark spin assignments of and 1/2. 

For definiteness, we compare top quark t and top squark t production since a consistent phe- 
nomenology exists for top squark pair production, obviating the need to invent a model of scalar quark 
production. Moreover, top squark decay may well mimic top quark decay. Indeed, if the chargino x + 
is lighter than the light top squark, as is true in many models of supersymmetry breaking, the dominant 
decay of the top squark is t — ► b + \ + . If there are no sfermions lighter than the chargino, the chargino 
decays to a W and the lightest neutralino \°- m another interesting possible decay mode, the chargino 
decays into a lepton and slepton, \ + — ► £ + v. The upshot is that decays of the top squark may be very 
similar to those of the top quark, but have larger values of missing energy and softer momenta of the 



visible decay products. A recent study for Run II of the Tevatron [289] concluded that even with 4 fb _1 
of data at the Tevatron, and including the LEP limits on chargino masses, these decay modes remain 
open (though constrained) for top squarks with mass close to the top quark mass. 

At the energy of the CERN LHC, production of it pairs and of pairs is dominated by gg subpro- 
cess, and the threshold behaviours in the two cases do not differ as much as they do for the qq incident 
channel. In Fig. 55(a), we show the partonic cross sections a(ys) as functions of the partonic sub- 
energy vi for the gg channel. In Fig. 55(b), we display the hadronic cross sections for pp — > tiX and 
pp — > ttX at proton-proton center-of-mass energy 14 TeV as a function of pair mass. We include the 
relatively small contributions from the qq initial state. After convolution with parton densities, the shape 
of the tt pair mass distribution is remarkably similar to that of the tt case. 

Based on shapes and the normalisation of cross sections, it is difficult to exclude the possibility that 
some fraction (on the order of 10%) of top squarks with mass close to 165 GeV is present in the current 
Tevatron tt sample. The invariant mass distribution of the produced objects, M t i, is quite different at the 
partonic level for the qq initial state (dominant at the Tevatron), but much less so for the gg initial state 
(dominant at the LHC). However, after one folds with the parton distribution functions, the difference in 
the qq channel at the Tevatron is reduced to such an extent that the M t i distribution is not an effective 
means to isolate top squarks from top quarks. 

Ironically, the good agreement of the absolute rate for tt production with theoretical expecta- 
tions I 



, 47] would seem to be the best evidence now for the spin- 1/2 assignment in the current Tevatron 



sample. 
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Fig. 56: Distribution in the invariant mass of a bottom quark and charged lepton (X) for a top quark or top squark decay, with 
relative size normalized to the cross sections at the LHC. The top squark decay and sparticle masses are discussed in the text. 

A promising technique to isolate a top squark with mass close to m t would be a detailed study 
of the momentum distribution of the top quark decay products (presumably in the top quark rest frame). 
One could look for evidence of a chargino resonance in the missing transverse energy and charged lepton 
momentum, or for unusual energy or angular distributions of the decay products owing to the different 
decay chains. One could also look for deviations from the expected correlation between angular distri- 
butions of decay products and the top spin [|167|]. 



As a concrete example of an analysis of this type, in Fig. 56 we present the distribution in the 
invariant mass X of the bottom quark and charged lepton, with X 2 = (p b + Pe+) 2 , where the bottom 
quark and lepton are decay products of either a top quark with m t = 175 GeV or a top squark t — ► 
X + b -» W + x°b -> £ + veX°b, with m~ t = 165 GeV, = 130 GeV, m^o = 40 GeV, and m b = 5 GeV. 
The X distribution is a measure of the degree of polarisation of the W boson in top quark decay [ ]290[ ], 
and the figure shows that the different dynamics responsible for top squark decay result in a very different 
distribution, peaked at much lower X. The areas under the curves are normalised to the inclusive ti and 
Urates at the LHC. 

In this simple demonstration potentially important effects are ignored such as cuts to extract the ti 
signal from its backgrounds, detector resolution and efficiency, and ambiguities in identifying the correct 
b with the corresponding charged lepton from a single decay. Detailed simulations would be required 
to determine explicitly how effective this variable would be in extracting a top squark sample from top 
quark events. Nevertheless, such techniques, combined with the large ti samples at the Tevatron Run II 
and LHC, should prove fruitful in ruling out the possibility of a top squark with mass close to the top 
quark mass, or alternatively, in discovering a top squark hidden in the top sample. 



B2 Direct Measurement of the Top Quark Electric Charge 

In order to confirm that the electric charge of the top quark is indeed Qtop = 2/3, one can either measure 
the charge of the fe-jet and W boson, or attempt to directly measure the top quark electro-magnetic 
coupling through photon radiation in 

pp^ti-f, pp->ti, t^Wbj. (85) 

Since the process pp—>ti~/ is dominated by gg fusion at the LHC, one expects that the tij cross section is 
approximately proportional to Q 2 op . For radiative top decays the situation is more complicated because 
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the photon can also be radiated off the 6-quark or the W line. 

The charge of the 6-jet can most easily be measured by selecting events where the 6-quarks are 
identified through their semi-leptonic decays, b^ivc with £ 1 = e, \i. The small semi-leptonic branching 
ratio of the 6-quark (Br(6— >£vc) as 10%) and wrong sign leptons originating from B — B mixing are the 
main problems associated with this method. For a quantitative estimate realistic simulations are needed. 
Nevertheless, we believe that the enormous number of top quark events produced at the LHC should 
make it possible to use semi-leptonic 6-tagging to determine the electric charge of the top quark. 

In our analysis, we focus on top charge measurement through the photon radiation processes listed 
in (p5|), concentrating on the lepton+jets mode, 

pp^jivjjbb. (86) 

We assume that both 6-quarks are tagged with a combined efficiency of 40%. Top quark and W decays 
are treated in the narrow width approximation. Decay correlations are ignored. To simulate detector 
response, the following transverse momentum, rapidity and separation cuts are imposed: 

Pr (6)>15GeV, \y(b)\<2, (87) 

p T (£) > 20 GeV, \r){£)\ < 2.5, (88) 

p T (j)>20GeV, \ V (j)\<2.5, (89) 

p T (7)>30GeV, 1^(7)1 < 2.5, (90) 

i> T > 20 GeV, all AR's > 0.4. (91) 

In addition, to suppress contributions from radiative W decays, we require that 

m (jjl) > 90 GeV and m T {£ T ^ T ) > 90 GeV, (92) 
where rriT is the cluster transverse mass of the £7 system. 



The events passing the cuts listed in ( ]8 8[ ) - ( J92D can then be split into three different subsamples: 

1. By selecting events which satisfy 

mibjj-y) > 190 GeV and m T (b£r,0T) > 190 GeV, (93) 

radiative top quark decays can be suppressed and an almost pure sample of ttj events is obtained 
("«7 cuts"). 

2. For 

m T (b h2 £r,M < 190 GeV and m(b 2jl jjj) > 190 GeV, (94) 

the process pp—>ti, t—>Wlr/, W^lv dominates ("t—>Wb~/, W^£v cuts"). 

3. Requiring 

m T (b lt2 £r^T) > 190 GeV and 150 GeV < m(b 2 ,ijn) < 190 GeV, (95) 

one obtains an event sample where the main contribution originates from the process pp—>ti, 
t^Wbj, W^jj ("t^Wlrf, W^jj cuts"). 

For m t = 175 GeV, Q tcyp = 2/3, and J Cdt = 100 fb~\ one expects about 2400, 11000 and 9400 events 
in the regions of phase space corresponding to the three sets of cuts. We have not studied any potential 
background processes. The main background should originate from jets production and should be 
manageable in a way similar to the W+ jets background for regular ti production. 

The differential cross section for the photon transverse momentum at the LHC is shown in Fig. |57[ 
Results are shown for m t = 175 GeV and three "top" quark charges: Qtop = 2/3, Qtop = —4/3, and 
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Fig. 57: The differential cross section for the photon transverse momentum in the reaction pp~* r yivjjbb at the LHC for three 
different "top" quark charges. 

Qtop = 1/3- For Qtop = —4/3, the "top" quark decays into a W~ and a 6-quark instead of t— >W + b. 
If Qtop = 1/3, the "6"-quark originating from the "top" decay is a (exotic) charge —2/3 quark In the 
tij region (Eq. ( |93| ) and Fig. |57|a), the pp^^lvjjbb cross section for a charge —4/3 (1/3) "top" quark 
is uniformly a factor « 3.3 larger 2.3 smaller) than that for Qtop = 2/3, reflecting the dominance 
of the gg—>ti / y process for which the cross section scales with Q\ op . On the other hand, for the pp—>ti, 
t^-Wbj, W~^lv selection cuts (Eq. ( |94| ) and Fig. |57|b), the cross section for Qtop = —4/3 is a factor 3 
to 5 smaller than that for a charge 2/3 top quark, due to destructive interference effects in the t— >Wb^/ 
matrix element. If Qtop = 1/3, the interference is positive, and the cross section is about a factor 2 to 2.5 
larger than for Qtop = 2/3. The results for the t—>Wbj, W—>jj selection cuts ((95)) are similar to those 
shown in Fig. |57|b, and are therefore not shown here. Note that the photon p-p distribution for radiative 
top decays is much softer than that for tij production. 

From our (simplified) calculation we conclude that the large number of double-tagged ^Ivjjbb 
events, together with the significant changes in the ttj and the tt, t—>Wlry cross sections should make it 
possible to accurately determine Qtop at the LHC. 

C APPENDIX: 4 th GENERATION QUARKS[] 

For completeness, we present here results for the total cross section of possible heavy quarks above the 
top quark mass. The scale and PDF dependences are shown in Fig. |58| The uncertainty due to the choice 
of scale is comparable to that of the tt cross section, although the effects of the higher order corrections 



are more and more important at large masses (see Fig. 59). The uncertainty induced by PDF changes 
becomes very large at large masses, in particular if one considers sets such as CTEQ5HJ which have 
harder gluons. Notice however that the relative effect due to the resummation corrections depends only 



very weakly upon the choice of PDF's (cf. Section 3.2) 



Section coordinator: M.L. Mangano 
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Fig. 58: Heavy quark total production rates. Left figure: scale dependence at fixed NLO (dashed lines in the lower inset), and 
at NLO+NLL (solid lines). Right figure: PDF dependence. See the Section 



3.2 



for details. 
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Fig. 59: Heavy quark total production rates. Left figure: fractional contribution induced by resummation contributions of order 
0(ag 4 ). Right figure: initial state composition. 

D APPENDIX: MONTE CARLO TOOLS0 

Dl Parton shower Monte Carlos 

General purpose Monte Carlo event generators like HERWIG, PYTHIA and ISAJET are essential tools 
for measuring the top quark cross section, mass and other production and decay properties. They are 



complementary to the QCD tools described in Section since, although they are less reliable for 
inclusive quantities like the total cross section, they provide a fully exclusive description of individual 
events at the hadron level. These can be analysed in exactly the same way as experimental data and can be 
put through full or fast detector simulations to estimate experimental systematics. In certain kinematic 
regions, particularly the quasi-elastic limit in which accompanying radiation is suppressed, they give 
more reliable QCD predictions than the available calculations. They include approximate treatments of 
higher order perturbative effects, hadronisation, secondary decays and underlying events. 

The three programs we discuss have the same basic structure, although the precise details vary 
enormously. Events are generated by starting with the hardest (highest momentum scale) interaction, 
described by exact QCD (or EW) matrix elements. This is usually only done to leading order so describes 

28 Section coordinators: M.L. Mangano, M. Seymour. 
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a 2^2 scattering process. The production of multi-parton final states is described as the emission of 
additional partons from the incoming and outgoing partons of the hard process. This is simulated by a 
parton shower algorithm in which the partons evolve downwards in some energy-like scale according to 
perturbatively-calculable probabilistic distributions. When the evolution scale becomes small the running 
coupling grows, phase space fills with (mostly soft) partons and perturbation theory breaks down. At this 
point a model of the non-perturbative physics is needed: the perturbative emission is cutoff by a fixed 
infrared cutoff and the system of partons is confined into hadrons. Having treated all outgoing partons 
we are left with the remnants of the incoming protons, stripped of the partons that participated in the hard 
process. These remnants can interact with each other, to produce additional soft hadrons in the event, 
known as the underlying event. 

Parton shower algorithms are developed by studying the amplitude to emit an additional parton 
into a given process. This is enhanced in two kinematic regions: collinear, where two massless partons 
are much closer to each other than any others or where a massless parton is close to the incoming proton 
direction; and soft, where a gluon is much softer than any other parton. In both cases the enhanced terms 
are universal, allowing a factorisation of emission by a system of partons from the process that produced 
them. In the collinear case, this factorisation works at the level of cross sections, so it is not surprising 
that a probabilistic approach can be set up. In the soft case however, the factorisation theorem is valid 
at amplitude level and it turns out that in any given configuration, many different amplitudes contribute 
equally. It therefore seems impossible to avoid quantum mechanical interference and so to set up the 
evolution in a probabilistic way. The remarkable result though is that, due to coherence between all the 
coloured partons in an event, the interference is entirely destructive outside angular-ordered regions of 
phase space. This means that the soft emission can be incorporated into a collinear algorithm, simply 
by choosing the emission angle as its evolution variable, as is done in HERWIG. The most important 
effects of coherence can be approximately incorporated by using some other evolution variable, like 
virtuality, and vetoing non-angular-ordered emission, as is done in PYTHIA. If the colour-coherence is 
not treated at all, one obtains the wrong energy-dependence of jet properties. Such models, like ISAJET, 
are completely ruled out by e + e~ annihilation data. Colour coherence effects are also important in 
determining the initial conditions for the parton evolution, resulting in physically-measurable inter-jet 
effects Q292T J, which are also in disagreement with ISAJET. 

Since the top quark decays faster than the typical hadronisation time, its width cuts off the parton 
shower before the infrared cutoff. Its decay then acts as an additional hard process and the resulting 
bottom quark (and two more partons if the W decays hadronically) continue to evolve. Additional coher- 
ence effects mean that radiation from the top quark is suppressed in the forward direction (the dead cone 
effect), as is radiation in the W direction in the top decay. These effects are again included in HERWIG, 
partially included in PYTHIA and not included in ISAJET. Since the top quark is coloured, the b quark in 
its decay is colour-connected to the rest of the event, meaning that its properties are not necessarily the 
same as in a 'standard' b production event. As mentioned in Section 4^6 and as discussed in more detail 
in j64|], such non-universal effects are small. 

Although parton showers are reliable for the bulk of emission, which is soft and/or collinear, it is 
sometimes the rare hard emissions that are most important in determining experimental systematics and 
biases. Such non-soft non-collinear emission should be well described by NLO perturbation theory, since 
it is far from all divergences. However, it is not straightforward to combine the advantages of the parton 
shower and NLO calculation, so it has only been done for a few specific cases. Most notable for hadron- 
hadron collisions are the Drell-Yan process, for which matrix element corrections are included in both 
HERWIG and PYTHIA, and top decay, which is included in HERWIG and discussed earlier in Section 4.62| 
in this report. The corrections to Drell-Yan events are particularly important at high transverse momenta, 
where the uncorrected algorithms predict far too few events. It is likely that implementing corrections to 
tt pair production would cure the analogous deficit at high p" seen in Fig. 

Hadronisation models describe the confinement of partons into hadrons. Although this process 
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is not well understood from first principles, it is severely constrained by the excellent data from LEP, 
SLD and HERA. The string model, used by PYTHIA, and the cluster model, used by HERWIG, both 
take account of the colour structure of the perturbative phase of evolution, with colour-connected pairs 
producing non-perturbative singlet structures that decay to hadrons. The biggest difference between these 
models is in how local these colour-singlet structures are. In the string model they stretch from a quark 
(or anti-di-quark) via a series of colour-connected gluons to an antiquark (or di-quark). In the cluster 
model each gluon decays non-perturbatively to a quark-antiquark pair and each resulting quark-antiquark 
singlet (coming one from each of two colour-connected gluons) decays to hadrons. The independent 
fragmentation model, used by ISAJET, on the other hand, treats each partem as an independent source of 
hadrons and is strongly ruled out by e + e~ data, for example on inter-jet effects in three-jet events, the 
so-called string effect. Of the other two models, PYTHIA gives the better description of e + e _ data, but 
HERWIG also gives an adequate description, despite having a lot fewer adjustable parameters. 

Models of the underlying event are not strongly constrained by either theoretical understanding or 
experimental data. Two extreme models are available and the truth is likely to lie between them. In the 
soft model, used in HERWIG, the collision of the two proton remnants is assumed to be like a minimum 
bias hadron-hadron collision at the same energy. A simple parametrisation of minimum bias data (from 
UA5 [293]) is used with little additional physics input. In the mini-jet model, used in PYTHIA and 
available as an additional package for HERWIG, on the other hand, the remnant-remnant collisions act as 
a new source of perturbative scattering, which ultimately produce the hadrons of the underlying event. 
To avoid regions of unstable perturbative predictions and problems with unitarity, a cutoff must be used, 
Pt,min ~ 1 GeV. Presumably for a complete description, some soft model should describe the physics 
below pt,min sucn that the results do not depend critically on its value. Unfortunately no such model 
exists at present. Although the two models give rather similar predictions for average properties of the 
underlying event, they give very different probabilities for the rare fluctuations that can be most important 
in determining jet uncertainties. This is an area that needs to be improved before LHC running begins. 



D2 Parton-level Monte Carlos 

With few exceptions (e.g. 3 or 4-jet final states in e + e~ collisions) multi-jet final states are not accurately 
described by the shower MC's described above. This is because emission of several hard and widely sep- 
arated partons is poorly approximated by the shower evolution algorithms, and exact (although perhaps 
limited to the tree level) matrix elements need to be used to properly evaluate quantum correlations. 
Parton-level Monte Carlos are event generators for multi-parton final states, which incorporate the exact 
tree-level matrix elements. They can be used for parton-level simulations of multi-jet processes, under 
the assumption that each hard parton will be identified with a final-state physical jet with momentum 
equal to the momentum of the parent parton. Selection and analysis cuts can be applied directly to the 
partons. In some cases, the partonic final states can be used as a starting point for the shower evolution 
performed using a shower MC such as HERWIG, PYTHIA, or ISAJET. For a discussion of the problems 
involved in ensuring the colour-coherence of the shower evolution when dealing with multi-parton final 
states, see [294]. 

In the following, we collect some information on the most frequently used parton-level MCs used 
in connection with top quark studies. 



D21 VECBOrf^ 

VECBOS [ ]150| ] is a Monte Carlo for inclusive production of a VF-boson plus up to 4 jets or a Z-boson 
plus up to 3 jets. VECBOS is therefore a standard tool used in the simulation of backgrounds to ti 
production. The matrix elements are calculated exactly at the tree level, and include the spin correlations 
of the vector boson decay fermions with the rest of the event. Various parton density functions are 

29 VECBOS authors: F.A. Berends, H. Kuijf, B. Tausk and W.T Giele. Contacts: giele@fnal.gov 
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available and distributions can be obtained by using the kinematics of the final state, available on an 
event-by-event basis together with the corresponding event weight. The code and its documentation can 
be obtained from: 



nttp : //www- theory . f nal . gov/people/giele/vecbos . html 



Documentation on the use of VECBOS within ATLAS can be found in [295] 



D22 CompHElf] 

CompHEP is a package for the calculation of elementary particle decay and collision properties in the 
lowest order of perturbation theory (the tree approximation). The main purpose of CompHEP is to gener- 
ate automatically transition probabilities from a given Lagrangian, followed by the automatic evaluation 
of the phase-space integrals and of arbitrary distributions. The present version has 4 built-in physical 
models. Two of them are the versions of the Standard Model (SU(3)xSU(2)xU(l)) in the unitary and 
t'Hooft-Feynman gauges. The user can change the models or even create new ones. 

The symbolic part of CompHEP allows the user to perform the following operations: 

1 . to select a process by specifying incoming and outgoing particles for the decays of 1 — > 2, . . . , 1 — > 
5 types and the collisions of 2— ► 2, . . . , 2 — > 4 types, 

2. to generate Feynman diagrams, calculating the analytical expressions for the squared matrix ele- 
ments, 

3. to save the algebraic symbolic results and to generate the optimized Fortran and C codes for the 
squared matrix elements for further numerical calculations. 

The numerical part of CompHEP allows to convolute the squared matrix element with structure functions 
and beam spectra, to introduce various kinematic cuts, to introduce a phase space mapping in order to 
smooth sharp peaks of a squared matrix element, to perform a Monte Carlo phase space integration by 
VEGAS, to generate events and to display distributions for various kinematic variables. Recently, an 



interface with PYTHIA has been created [283]. This allows to perform realistic simulations of the process 



including hadronisation effects as well as the effects of the initial and final state radiation. 

The CompHEP codes and manual are available from the following Web sites: 
|http : / /theory, npi.msu. su/~comphep 



http: //www. ifh. de/~pukhov 



D23 ALPHAfJ\ 

ALPHA is an algorithm introduced in [ ]296| ] for the evaluation of arbitrary multi-parton EW matrix el- 
ements. This algorithm determines the matrix elements from a (numerical) Legendre transform of the 
effective action, using a recursive procedure which does not make explicit use of Feynman diagrams. The 
algorithm has a complexity growing like a power in the number of particles, compared to the factorial- 
like growth that one expects from naive diagram counting. This is a necessary feature of any attempt 
to evaluate matrix elements for processes with large numbers of external particles, since the number of 
Feynman diagrams grows very quickly beyond any reasonable value. 



An implementation of ALPHA for hadronic collisions was introduced in Q294[ ], where the algorithm 
was extended to the case of QCD amplitudes (see also Q297| ]). The main aim of the hadronic version of 
ALPHA is to allow the QCD parton-shower evolution of the multi-parton final state, in a way consistent 
with the colour-coherence properties of the soft gluon emission dynamics. This is achieved by evaluating 



the QCD amplitudes in an appropriate colour basis [294], such that the assignement of a specific colour 



flow configuration on an event-by-event basis. The pattern of colour flow defines the colour currents 



30 CompHEP authors: A. Pukhov, E. Boos, M. Dubinin, V. Edneral, V. Ilyin, D. Kovalenko, A. Kryukov, V. Savrin, 
S. Shichanin, A. Semenov. Contacts: pukhov@theory.npi.msu.su, ilyin@theory.npi.msu.su 

31 ALPHA authors: F. Caravaglios, M. Moretti. The version for hadronic collisions received additional contributions from 
M.L. Mangano and R. Pittau. Contact: moretti@fe.infn.it 
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required to implement the angular ordering prescription which embodies, at the leading order in the 1/N C 
expansion, the quantum coherence properties of soft-gluon radiation, as discussed in Appendix pl[ A 



version of the code is being completed [ 298 ] , which incorporates the evaluation of Wbb + n jets (n < 4), 
with all 6-mass effects included. This program will allow a complete evaluation of the W+ multijet 
backgrounds to single top and ti production. The code contains 3 modules: the first for the generation 
of parton-level events, with the assignement of partonic flavours, helicities and colour flows. The second 
for the unweigthing of the events, and the third for the parton-shower evolution of the initial and final 
states, done using the HERWIG MC. The code will soon be available from the URL: 



http : / /home . cern . ch/~mlm/ alpha 
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